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ABSTRACT 


This  research  dealt  with  the  measurement  of 

subnanosecond  nuclear  lifetimes  using  the  pulsed  beam 

delayed-coincidence  technique.  Measurements  were 

7 

performed  on  isotopes  in  the  f  shell  and  specifically 

the  isotopes  of  Titanium  and  Vanadium.  Experimental 

investigations  were  also  pursued  in  two  other  isotopes, 
59  .  65 

Ni  and  Zn.  Several  new  lifetimes  were  determined 
and  confirmation  was  obtained  for  some  previous  values 
which  were  measured  with  different  techniques.  More 
information  was  also  obtained  on  certain  levels  where 
previous  results  are  in  disagreement. 

A  new  technique  has  been  developed  to  provide 
simultaneous  data  acquisition  of  both  delayed  and  prompt 
time  distributions.  The  centroids  of  the  prompt  time 
distributions  were  used  to  define  the  so-called  prompt 
time  centroid  curve.  The  Centroid  Shift  technique  was 
employed  to  calculate  the  lifetimes  from  the  centroid 
shifts  of  the  delayed  time  distributions  with  respect 
to  the  prompt  time  centroid  curve  at  the  appropriate 
energies.  Additional  'prompt'  reference  y-rays  were 
obtained  by  using  a  second  target.  For  this  purpose, 
a  two-target  rotating  system  was  developed.  The  system 
has  the  electronics  capable  of  routing  into  separate 
regions  of  computer  memory,  information  obtained  from 
each  target.  Thus  errors  introduced  by  background 
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subtractions  are  minimised  since  the  peak  to  background 
ratio  is  not  adversely  affected  by  the  presence  of  the 
second  target. 

A  major  part  of  the  work  involved  the  improve¬ 
ment  of  the  time  resolution  and  study  of  the  timing 
properties  of  coaxial  Ge(Li)  detectors.  The  significant 
improvement  in  the  time  resolution  was  achieved  through 
the  development  of  the  Triple  Constant  Fraction  Discrimi¬ 
nator,  an  electronic  circuit  designed  to  improve  the 
time  resolution  by  a  pulse  shape  discrimination  technique. 
The  improvement  in  the  time  resolution  leads  to  greater 
accuracy  in  the  determinations  of  the  centroids  of  the 
time  distributions,  which  is  the  main  limitation  of  the 
Centroid  Shift  technique  for  measuring  lifetimes  in  the 
tens  of  picosecond  region. 

Finally  a  theoretical  description  of  the  energy 
levels  and  electromagnetic  transition  rates  of  the 
isotopes  studied  was  performed  using  the  Strong  Coupling 
Symmetric  Rotator  Model,  which  included  band  mixing,  and 
the  Intermediate  Coupling  Model.  The  predicted  level 
schemes  are  in  good  agreement  with  experimental  results 
and  fairly  good  agreement  is  obtained  between  experi¬ 
mental  and  theoretical  electromagnetic  transition  rates. 
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CHAPTER  1 


INTRODUCTION 


The  knowledge  of  the  lifetime  of  a  nuclear  state 
is  important  in  the  determination  of  the  electromagnetic 
transition  matrix  elements  for  comparison  with  predic¬ 
tions  based  on  theoretical  nuclear  models.  These  matrix 
elements  which  involve  an  overlap  of  the  wavef unctions 
of  the  states  between  which  the  transition  occurs  are 
more  sensitive  to  the  details  of  the  nuclear  model  than 
many  other  quantities  which  could  be  measured.  To  date, 
there  is  only  a  limited  knowledge  of  the  wavef unctions 
of  nuclear  levels  and  this  essentially  in  the  low  mass 
region.  Thus  there  is  interest  in  the  measurement  of 
lifetimes  of  nuclear  states  to  obtain  more  knowledge 
about  the  form  of  the  wavef unctions . 

For  lifetime  determinations  in  the  range 
-3  -12 

10  >  x  >  10  secs,  the  most  common  method  is  the 

so-called  delayed-coincidence  technique  which  involves 
the  detection  of  the  times  of  population  and  decay  of 
the  state  whose  lifetime  is  to  be  measured.  In  the 
most  widely  used  arrangement,  scintillator  photo¬ 
multiplier  combinations  determine  both  the  time  of 
emission  and  the  energy  of  the  radiations  announcing 
the  birth  and  decay  of  the  state  as  illustrated  in 
fig.  1.  The  current  literature  contains  many  examples 
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Fig  .  1 


Schematic  diagram  of  the  various  delayed  coincidence 
techniques ,  data  acquisition  and  display  of  the 
prompt ,  P(t),  and  delayed,  D(t)  time  distributions. 
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of  y-y  and  3_Y  delayed-coincidence  measurements.  How¬ 
ever,  this  method  limits  the  number  of  nuclear  states 
that  can  be  investigated  due  to  the  poor  energy  resolu¬ 
tion  of  scintillator  detectors  and  the  number  of  states 
that  are  populated  by  3  and  y-ray  transitions. 

Alternatively,  the  population  of  a  state  can  be 
achieved  by  using  a  pulsed  beam  to  excite  the  state  via 
a  nuclear  reaction,  the  timing  signals  coming  from  a 
capacitive  beam  pick-off  before  the  target  and  from 
the  detector  corresponding  to  the  detection  of  the 
radiation  from  the  decay  of  the  state.  The  University 
of  Alberta  has  a  7 . 5MV  van  de  Graaff  which  has  a  Mobley 
Compression  System  capable  of  producing  every  microsecond 
a  beam  burst  with  a  width  of  0.4ns.  This  system  is 
ideal  for  the  measurement  of  lifetimes  by  direct  timing 
using  the  delayed-coincidence  technique. 

With  the  advent  of  Lithium  drifted  Germanium 
[Ge(Li)]  detectors  with  their  superior  energy  resolution 
capability,  one  of  the  major  difficulties  associated  with 
the  commonly  employed  scintillator-photomultiplier 
combinations  can  be  eliminated.  However,  very  little 
work  has  been  reported  utilizing  coaxial  Ge(Li)  detec¬ 
tors  for  lifetime  measurements  in  delayed-coincidence 
systems,  especially  pulsed  beam  experiments  (Ya  69, 

Wa  70,  Ki  71),  due  to  the  poor  time  resolution  charac¬ 
teristics  of  Ge(Li)  detectors.  It  is  clear,  however, 
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that  their  extremely  good  energy  resolution  would  give 
a  great  advantage  over  scintillator  detectors  for 
lifetime  measurements  of  excited  states  in  nuclei  with 
complicated  decay  schemes. 

In  general,  two  different  time  distribution 
curves  are  obtained  in  delayed-coincidence  experiments 
as  illustrated  in  fig.  1.  One,  the  prompt  distribution, 
P(t),  is  obtained  from  events  in  which  a  state  is  produced 
and  decays  during  an  interval  short  compared  to  the 
time  resolution  of  the  apparatus.  The  other,  called  the 
delayed  distribution  ,  D  (t)  characterises  the  lifetime  to 
be  measured.  For  lifetimes  which  are  much  greater  than 
the  time  resolution  of  the  apparatus,  the  delayed  dis¬ 
tribution  exhibits  an  exponential  tail.  However,  for 
lifetimes  shorter  than  the  time  resolution  of  the 
measuring  apparatus  the  delayed-coincidence  method  is 
restricted  to  experimental  observations  of  the  centroid 
shift  between  the  delayed  and  prompt  time  distribution 
curves  for  a  given  energy. 

The  accuracy  of  the  centroid  shift  method,  how¬ 
ever,  is  severely  limited  by  the  time  resolution  of 
Ge(Li)  detectors.  If  electronic  noise  were  the  only 
limitation,  timing  with  Ge(Li)  detectors  could  be  very 
precise,  particularly  at  high  energies  since  the  full 
width  at  half  maximum  (FWHM)  of  the  time  distribution 
curves  would  scale  inversely  as  the  detected  energy. 


However,  the  time  resolution  is  seriously  degraded  by 
the  pulse  shapes  and  the  pulse  rise-time  variations 
which  depend  on  the  y-ray  interaction  position  in  the 
detector  intrinsic  volume,  the  detector  geometry  and 
the  electric  field  as  a  function  of  the  position  inside 
the  detector. 

The  purpose  of  this  research  was  the  significant 
improvement  of  the  time  resolution  of  coaxial  Ge(Li)  de¬ 
tectors  so  that  lifetimes  as  short  as  10  secs  could 
be  measured  using  the  pulsed  beam  facility.  The  Doppler 
Shift  Attenuation  Method  (DSAM)  (Li  63,  Wa  63)  and  the 
Recoil  Distance  Method  (RDM)  (A1  65,  A1  70)  are  other 

techniques  used  to  measure  lifetimes  in  this  region.  How 

-12  -11 

ever,  in  the  range  10  <t<5x10  secs,  the  DSAM  is  res 

tricted  to  gas  targets  and  the  RDM  requires  high  energy 
heavy  ion  beams  to  provide  sufficiently  large  recoil  vel 
ocities . 

The  improvement  in  the  time  resolution  was  accom¬ 
plished  through  the  development  of  the  so-called  Triple 
Constant  Fraction  Discriminator,  an  electronic  circuit 
designed  to  improve  the  time  resolution  by  a  pulse  shape 
discrimination  technique  (Wh  74,  Me  74). 

The  accuracy  of  the  centroid  shift  technique  is 
also  determined  by  the  definition  of  the  so-called 
prompt  time  distribution  centroid  curve,  a  curve  through 
the  points  representing  the  centroids  of  the  prompt  time 
distributions  plotted  against  the  corresponding  energy 
of  the  'prompt'  y-ray.  Since  a  shift  from  this  curve 
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corresponds  to  a  y-ray  with  a  measurable  lifetime,  it 
is  essential  that  the  curve  be  defined  accurately  in 
the  region  of  interest.  However,  in  most  cases  the 
isotope  under  study  does  not  have  a  sufficient  number 
of  'prompt'  deexcitation  y-rays  in  the  required  energy 
region.  To  overcome  this  problem,  a  rotating  target 
system  was  developed.  To  the  rotator  are  attached  two 
targets,  each  rotating  alternately  in  and  out  of  the 
pulsed  beam,  the  one  target  producing  the  energy  level 
whose  lifetime  is  to  be  measured  and  the  other,  supple¬ 
mentary  reference  'prompt'  y-rays  in  the  energy  region 
of  interest. 

This  simultaneous  acquisition  of  data  for  the 
prompt  and  delayed  time  distributions  is  a  new  develop¬ 
ment  of  the  technique.  In  most  previous  delayed-coincidence 
techniques  alternate  experimental  runs  were  done  to 
obtain  the  prompt  and  delayed  time  distributions.  Thus 
errors  due  to  electronic  drifts  and  shifts,  and  pulsed 
beam  conditions  between  separate  runs,  are  overcome  by 
our  technique.  Furthermore,  by  the  definition  of  the 
prompt  time  distribution  centroid  curve,  centroid  shifts 
are  calculated  with  respect  to  this  curve  at  the  appro¬ 
priate  energy  whereas  previously  this  was  done  with 
respect  to  a  prompt  time  distribution  from  a  y-ray  of 
approximately  the  same  energy. 
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The  collection  of  data  by  our  technique  is  made 
possible  by  a  two  dimensional  sort  routine,  a  software 
component  of  the  on-line  Honeywell  516  computer. 

Using  this  routine  the  total  time  distribution  spectrum 
can  be  sorted  by  setting  up  to  32  digital  windows  on  the 
energy  spectrum  and  on  that  region  of  the  'pulse  shape' 
spectrum  which  produces  improved  time  resolution. 

A  further  refinement  of  the  technique  was  the 
development  of  an  electronic  Target  Discriminator  module 
which  can  rout  the  information  from  each  target  into 
separate  regions  of  memory  of  the  computer.  In  this 
way  the  peak  to  background  ratio  in  the  energy  spectrum 
is  not  affected  by  the  presence  of  the  second  target,  a 
consideration  which  is  very  important  for  the  reduction 
of  errors. 

Measurements  were  performed  on  isotopes  in  the 
f  7/2  shell  especially  the  isotopes  of  Titanium  and 

Vanadium.  Investigations  in  two  other  isotopes, 

65  59  . 

Zn  and  Ni,  also  proved  fruitful  with  several  life¬ 
time  determinations  being  made.  Level  schemes  and 
electromagnetic  transition  probabilities  for  the  iso¬ 
topes  studied  were  compared  with  theoretical  results 
obtained  from  the  Intermediate  Coupling  Model  and  the 
Strong  Coupling  Symmetric  Rotator  Model,  using  the 
Nilsson  Hamiltonian  to  represent  the  interaction  between 
the  last  unpaired  nucleon  and  the  core. 


CHAPTER  2 


THE  ANALYSIS  OF  DELAYED-CO INCIDENCE  EXPERIMENTS 

The  so-called  delayed-coincidence  technique  is 
frequently  used  to  determine  the  lifetimes  of  excited 
nuclear  states  by  employing  a  Time  to  Amplitude  Converter 
(TAC) ,  which  translates  time  differences  into  voltage 
signals,  and  an  Analogue  to  Digital  Converter  (ADC)  which 
digitises  these  signals.  There  are  several  different 
methods  of  determining  the  unknown  lifetime  from  the 
two  time  distributions  experimentally  obtained.  These 
are  commonly  referred  to  as  the  prompt  and  delayed  time 
distributions . 

A  frequently  used  method  is  the  so-called  slope 
method  (Ne  50) ,  in  which  a  fit  of  a  pure  exponential 
function  to  that  part  of  the  delayed  curve  well  outside 
the  region  of  the  prompt  time  curve  yields  the  desired 
decay  constant,  A.  This  method  can  be  applied  if  the 
lifetime  is  sufficiently  long  to  give  a  significant 
exponential  tail,  but  otherwise  is  statistically  poor 
in  that  only  part  of  the  data  is  used  in  the  evaluation 
of  the  lifetime.  Statistically,  the  methods  utilizing 
the  moments  of  the  measured  time  distributions  (Ba  50, 

Su  62)  provide  a  more  accurate  determination  of  the  life¬ 
time  since  the  entire  data  are  used  in  the  calculations. 
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The  most  important  of  these  methods  is  the  first 
moment  or  centroid  shift  method;  however,  higher  moments 
can  also  yield  information  about  the  lifetimes.  Another 
method  which  has  statistical  errors  the  same  magnitude 
as  those  of  the  centroid  shift  method  is  the  integral 
transform  method  (Sp  66A) . 

A  theoretical  description  of  the  above  methods 
of  analysis  is  presented  below  with  a  discussion  of  the 
advantages  and  disadvantages  of  each.  The  prompt  and 
delayed  distributions  are  denoted  by  P(x)  and  D(x),  and 
the  corresponding  histograms  obtained  from  the  ADC  are 
denoted  by  and  D^.  P(x)  and  D(x)  are  assumed  iden¬ 
tically  zero  for  x  <  0  and  the  number  of  counts  under  the 
two  curves  (their  areas)  are  Np  and  N  ,  respectively. 


2 . 1  The  Slope  Method 

The  mathematical  relation  between  the  two  time 
distributions  is  the  convolution  integral 


D(x)  = 


N 

N, 


D 


w  (t )  P  (x  -  t )  dt 


2-1 


0 


where  if  the  decay  is  a  radioactive  decay  with  a  mean- 
life  x  (=  1/A),  the  decay  function  is 


w (t)  =  -  exp (-  t/T ) 


2-2 


Substitution  of  eg.  2-2  in  2-1  gives 
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D(x)  = 


N 

N 


D 
P 


A  e  At  P(x  -  t)  dt 


0 


N 

N, 


D 


x 
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P  (y)  A  e  A  (x  Y)  dy 
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where  y  =  x  -  t.  Differentiating  eq.  2-3  with  respect 
to  x  gives 


dD  (x) 
dx 


ND 

A  [r=—  P  (x)  -  D  (x)  ] 
P 


From  this  relation  follows 


dx 


[ In  D  (x)  ]  =  -  A  [  1  - 


P(x)  ND 
D(x)  N 


] 
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Consequently,  the  decay  constant  can  be  found  from  the 
slope  of  the  logarithmic  plot  of  the  delayed  curve  in 
the  region  where  P(x)/Np  is  negligible  compared  to 
D (x ) /Np .  Thus  only  that  part  of  the  delayed  curve  which 
is  outside  the  prompt  curve  can  be  used  in  the  evaluation 
of  the  lifetime.  This  restricts  lifetime  measurements  by 
this  method  to  those  which  are  greater  than  the  time 
resolution  of  the  apparatus  used.  The  method,  however, 
is  insensitive  to  electronic  drifts  and  the  prompt 


admixtures. 
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2 . 2  First  Moment  or  Centroid  Shift  Method 
th. 

The  n  moment  of  the  delayed  time  distribution 
is  defined  through  the  formula 


M  (D)  = 
n 


n 


x  D  ( x )  dx 


2-5 
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Using  n  =  1  and  substituting  in  eq.  2-1 


N 

M1(D>  =  N 


D 


x  dx 


w  ( t )  P  (x  -  t )  dt  . 


0 
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Substituting  y  =  x  -  t  and  rearranging  one  obtains 
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M1<D)  =  N 
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y  P  (y)  dy 


w  (t )  dt  + 


N 
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P (y)  dy 


t  w(t)  dt 
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—  M  (P ) 
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Since  M  (w)  =  n! 
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Mx  (D) 


N 


D 


+  Nd  (w)  . 

xn  it  follows 
MX(P) 
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Thus  the  lifetime  to  be  measured  is  the  centroid  shift 
between  the  delayed  and  prompt  time  distributions.  In 
practice  the  formula  takes  the  form 


T 


Z  D 


2-7 


and  the  statistical  error  can  be  calculated  using 
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where  A  =  f(a,b,c, 
variables  a,  b,  c, 


)  is  a  function  of  the  independent 
.  Therefore  the  statistical  error 


m  t  is 


s  = 


sx±  Di 


(IX.  D.)2  p± 
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Equations  2-7  and  2-9  are  independent  of  the  origin  of 
the  coordinate  system  and  if  the  particular  choice  of 
Zxi  P.  =  0  is  made  expression  2-9  can  be  transformed 
using  eq.  2-7  and  2-11  to 


s  M2  (P )  2.  i  i 

2  lN^  N^j  N 


_T 

2 


2-10 


T 


D 


D 


where  M^CP)  is  the  normalised  second  moment  of  the  prompt 

curve  about  its  own  centroid  and  is  a  measure  of  the 

width  of  the  curve  (normalised  means  adjusted  to  unit  area) 
c  2 

If  M ^  (P ) /t  is  much  larger  than  unity  i.e.  the 
lifetime  is  small  compared  to  the  time  resolution  of 
the  apparatus,  then  the  statistical  error  is  approximately 
equal  to 


M2<P>^  +  t ^ 


This  stresses  the  importance  of  good  time  resolution  when 
short  lifetimes  are  to  be  measured.  It  is  for  this 
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reason  that  a  large  section  of  this  thesis  is  devoted 
to  the  improvement  of  the  time  resolution  of  coaxial 
Ge(Li)  detectors. 

On  the  other  hand,  when  x  is  large  the  relative 
error  is  approximately  (N  ) *  2  and  nearly  independent  of 

the  shape  of  the  prompt  curve,  the  slope  method  is  then 
better  in  these  cases.  The  centroid  shift  method  is 
simple  to  use  and  has  small  statistical  errors,  however, 
electronic  shifts  and  drifts  play  an  important  role 
since  they  give  rise  to  systematic  errors.  These  pro¬ 
blems  can  be  overcome  if  the  prompt  and  delayed  time 
distributions  are  obtained  simultaneously. 


2 . 3  Methods  of  Higher  Moments 

Using  the  definition  of  the  moments  eq.  2-5  one 

obtains 
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00 
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Substituting  y  =  x  -  t  and  rearranging 
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Therefore  using  relation  2-6 


m2  (D) 


N 


M  (P)  M-.  (P) 

+  2x  [- - - -  +  x] 
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N, 


N, 


m2(p)  m1(d) 
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N, 
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In  an  analogous  way  the  relation  for  the  third 
moment  can  be  obtained  and  the  meanlife  from  the  result¬ 


ing  formulae  are: 


2nc*  moment: 
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m2  (D) 


2Mx(D) 


N 


D 


m2  (P) 

n" 


2-12 


,rd 


N 


moment :  x  = 


D 


3M2 (D) 


M  (D)  M  (P) 

- 4—] 


N 


D 


N, 


2-13 


Again  the  origin  of  the  coordinate  system  is 
arbitrary  and  by  choosing  the  origin  in  such  a  way  that 
M-^(P)  =  0  the  formulae  for  the  statistical  errors  reduce 
to  simpler  forms  (Sp  66)  as  shown  below.  These  are  for 
the  second  moment 
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m4  (P) 

-  — 4 
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N, 


4M2  (P) 
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2 


2-14 


and  for  the  third  moment 
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Another  version  of  the  higher  moments  was  proposed 
by  Sundstrom  (Su  62)  in  which  the  normalised  moments  of 
the  delayed  and  prompt  curves  (both  curves  assumed 
normalised  to  unit  area)  are  computed  with  respect  to 
their  centroids  i.e. 


M°(g> 
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Using  this  relationship  formula  2-16  holds 


M?(D)  =  M? (P )  +  M?(w)  i  =  1,2,3  2-16 

and  the  meanlife  is  evaluated  from  eq.  2-16  to  give 


2nd  moment:  t  =  [M^CD)  -  M^CP)]  2  2-17 

3rd  moment:  x  =  ^-[M^  (D)  -  M^(P)]lf  2-18 


As  pointed  out  by  Weaver  and  Bell  (We  60)  the 
quantity  M^(P) ,  the  third  moment  of  the  prompt  curve 
P  (x)  about  its  own  centroid,  is  usually  small  and 
frequently  negligible  so  that  for  many  cases  it  is 


sufficient  to  write 
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Thus  it  is  possible  to  evaluate  the  required  lifetime 
from  the  shape  of  the  delayed  curve  with  only  rough 
knowledge  of  the  shape  of  the  prompt  curve. 


The  great  advantage  of  the  higher  moment  methods 


is  that  if  the  delayed  and  the  prompt  time  distribution 
curves  are  recorded  under  different  experimental  con¬ 
ditions  and  a  pure  shift  (translation) results ,  these 
methods  eliminate  such  effects.  Electronic  drifts, 
however,  can  have  serious  effects  especially  for  the 
second  moment  method  as  they  result  in  broadening  of 
the  two  time  distribution  curves.  Also  since  these 
methods  depend  more  on  the  outer  less  well-defined 
parts  of  the  curves,  these  methods  suffer  from  inaccuracy 
in  the  determination  of  accidentals. 

2 . 4  Integral  Transform  Method 

In  this  method,  developed  by  Sparrman  (Sp  66A) , 
the  Laplace  integral  transforms  of  P(x)  and  D(x)  are 
defined  as 
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Using  the  measured  distributions  P.  and  D.  and 

1  1 

a  selected  value  of  r,  numerical  values  of  L(P)  and 
L (D )  are  obtained  by  numerical  integration.  The 
choice  of  r  is  arbitrary  except  for  the  condition, 
r  <  Xr  the  decay  constant.  Using  eq .  2-3  the  Laplace 
integral  can  be  obtained. 
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Changing  the  order  of  integration  and  evaluating  the 


inner  most  integral  yields 
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and  hence 
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Using  eq.  2-8  the  statistical  error  in  x  is  given  by 
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As  c  =  Np/Np ,  {sc/c}  is  easily  obtained: 
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For  the  purpose  of  calculating  the  errors  in  L(P)  and 
L (D )  the  values  of  L(P)  and  L(D)  are  approximated  by  the 


sums 
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Using  eq.  2-8  and  the  fact  sD 
expression  is  obtained 
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Dn ,  the  following 
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error,  s  ,  therefore  depends  on  r.  Sparrman  (Sp  66A) 
has  shown  that  for  a  minimum  statistical  error  an  r- 
value  roughly  equal  to  A/2  is  the  optimum  choice  in  most 
cases . 

To  anticipate  the  statistical  errors  in  the  above 

methods  Sparrman  and  Falk  (Sp  66)  assumed  that  the  prompt 

2 

curve  was  of  the  form  cos  xn/2a  (FWHM  =  a)  and  plotted 
the  relative  errors  against  a/x  as  shown  in  fig.  2. 

The  prompt  curve  was  chosen  to  have  this  form  since  the 
experimental  prompt  curves  have  approximately  this  form 
and  partly  because  the  convolution  integral  could  be 
found  analytically. 


In  conclusion  the  first,  second,  and  third  moments 


can  be  interpreted  in  terms  of  apparent  features  of  the 
delayed  curve  as  compared  to  the  prompt  curve.  The 
first  moment  is  sensitive  to  the  position  of  the  centroids 
on  the  time  axis,  the  second  moment  to  the  width  of  the 
curve,  and  the  third  to  the  asymmetry  of  the  curve. 
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Fig.  2.  The  relative  statistical  error  in  the  mean- 
life,  t,  for  the  moment  methods  and  the 
integral  transform  method  as  a  function  of 
the  FWHM  of  the  prompt  curve  over  t.  The 
prompt  curve  is  assumed  to  be  of  the  form 
cos27Tx/2a.  (FWHM  =  a,  ND=Np)  (Sp  66). 
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Even  higher  moments  can  be  employed  but  the  statistical 
errors  of  the  moments  increase  with  the  order  and  thus 
only  a  limited  number  of  moments  are  useful  for  practical 
purposes . 


CHAPTER  3 


THE  IMPROVEMENT  IN  THE  TIME  RESOLUTION  OF 
COAXIAL  Ge(Li)  DETECTORS 

The  improvement  in  the  time  resolution  of 
Ge(Li)  detectors  is  crucial  to  the  measurement  of  life¬ 
times  shorter  than  10  ^  secs  by  the  Centroid  Shift 
technique.  As  shown  in  Chapter  2,  when  the  lifetime 
to  be  measured  is  small  compared  to  the  time  resolution 
of  the  equipment,  the  statistical  error  is  roughly  equal 
to  M^tPHl/Np  +  1/Np]  .  M^(P)  is  the  normalised  second 

moment  of  the  prompt  curve  about  its  own  centroid  and 
Np  and  Np  are  the  total  number  of  counts  in  the  delayed 

and  prompt  time  distribution  curves,  respectively. 

c 

Since  IV^  (P)  is  a  measure  of  the  width  of  the  curve,  this 
stresses  the  importance  of  good  time  resolution  when 
short  lifetimes  are  to  be  measured. 

Ge(Li)  detectors,  with  their  extremely  good  energy 
resolution,  are  essential  for  lifetime  measurements  in 
nuclei  with  complicated  decay  schemes.  However,  the 
time  resolution  of  Ge(Li)  detectors  is  seriously 
degraded  by  the  pulse  shape  and  pulse  rise-time  varia¬ 
tions  which  depend  on  the  y-ray  interaction  position  in 
the  detector  intrinsic  volume,  the  detector  geometry  and 
the  electric  field  as  a  function  of  the  position  inside 
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the  detector.  As  a  result  one  observes  a  wide  variation 
of  pulse  shapes  and  pulse  rise-times  as  indicated  in 
fig.  3, which  is  a  photograph  of  the  actual  signals 
obtained  from  the  preamplifier  of  a  23c. c.  coaxial 
Ge(Li)  detector. 

Several  methods  (Ch  68,  Ja  69,  Br  69,  Mo  70, 

Be  72,  Ch  72)  have  been  introduced  to  compensate  for 
the  time  walk  due  to  rise-time  changes  in  order  to 
improve  the  time  resolution  of  Ge(Li)  detectors.  The 
ARC  zero-crossing  method  (Amplitude  and  Rise-time 
Compensation)  introduced  by  Chase  (Ch  68)  is  one  of  the 
most  successful  methods  to  date  for  improving  the  time 
resolution  in  Ge(Li)  detectors.  The  method,  which  uses 
a  Constant  Fraction  of  Pulse  Height  Trigger  (Ge  67) ,  was 
originally  intended  for  use  with  scintillator  detector 
pulses  which  have  uniform  shapes.  Chase  has  shown  that 
the  same  electronic  technique  applied  to  the  Ge(Li) 
detector  pulses  does  a  first  order  correction  to  the 
rise-time  variations  in  the  pulses.  In  fact,  it  corrects 
exactly  for  pulses  which  have  an  initial  rising  portion 
which  is  linear.  However,  the  early  slope  of  the  time 
distribution  curves  are  often  distorted  significantly 
due  to  slope  changes  that  occur  in  the  pulses  before  the 
zero-crossing  of  the  ARC  shaped  pulses.  Thus  a  second 
order  correction  would  be  useful  to  further  improve  the 
time  resolution  of  Ge(Li)  detectors. 
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A  theoretical  base  to  our  technique  for  improving 
the  time  resolution  is  formed  by  an  analysis  of  the 
pulse  shapes  from  coaxial  Ge(Li)  detectors  as  a  function 
of  the  position  of  irradiation  in  the  detector,  and  the 
synthesis  of  these  pulses  in  a  Constant  Fraction  of  Pulse 
Height  Trigger.  The  results  obtained  in  our  investiga¬ 
tions,  which  involve  a  pulse  shape  discrimination  tech¬ 
nique,  are  considerably  better  than  those  of  previous 
methods  and  are  such  that  lifetimes  as  short  as  10  secs 
have  been  measured  by  the  Centroid  Shift  technique. 

In  this  work  some  interesting  problems  developed 
which  induced  a  more  detailed  investigation  of  the  tim¬ 
ing  properties  of  coaxial  Ge(Li)  detectors.  The  results 
of  this  investigation  will  be  discussed  later  in  this 
chapter  and  Chapter  5. 

3 . 1  Pulse  Shape  Analysis 

The  timing  properties  of  semiconductor  detectors 
are  intimately  related  to  the  pulse  shapes  (charge  or 
current  pulses)  which  the  detector  delivers  to  the 
amplifier  as  a  consequence  of  the  detector-radiation 
interaction.  An  interaction  in  the  sensitive  region 
of  a  detector  releases  a  large  number  of  electron-hole 
pairs  which  are  separated  and  collected  by  the  electric 
field  present  in  the  intrinsic  region.  The  pulse  shapes 
depend  on  the  charge  transit  time  which  is  determined  by 
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the  following  major  factors  (Be  68): 

a)  The  electric  field  as  a  function  of  the  posi¬ 
tion  inside  the  detector  sensitive  region, 

b)  the  electron  and  hole  mobilities,  y^  and  y^, 
which  in  a  given  material  are  functions  of  the  electric 
field  and  temperature, 

c)  the  distribution  of  charge  created  by  the 
detected  radiation  which  is  dependent  on  the  energy 
and  the  type  of  incident  radiation. 

Other  factors  such  as  the  noise  or  charge  trapp¬ 
ing  effects,  due  to  impurities  and  defects  in  the 
semiconductor  material,  lead  to  additional  requirements 
for  the  total  collection  of  the  charges.  At  higher 
energies  the  plasma  effect,  which  is  the  production  of 
secondary  electron-hole  pairs  by  the  moving  charges 
from  the  initial  interaction,  also  introduces  further 
complications  to  the  charge  collection  time. 

In  the  following  analysis  we  consider  an  ideal 
coaxial  semiconductor  detector  in  which  no  charge 
trapping  or  recombination  effects  occur.  The  equation 
of  motion  of  a  carrier  is  given  by  the  empirical  for¬ 
mula  stated  by  Moszyriski  and  Bengtson  (Mo  72A) 

v  =  vs/[l  +  vs/ye]  3-1 

where  v  is  the  drift  velocity,  v  is  the  carrier  sat- 

s 

uration  velocity,  e  is  the  electric  field  and  y  is  the 
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zero  field  carrier  mobility.  However,  this  formula  leads 
to  complicated  numerical  integration  in  order  to  deter¬ 
mine  the  pulse  shapes  from  coaxial  detectors.  To 
simplify  the  mathematics  we  neglect  the  carrier  velocity 
saturation  effect  and  assuming  that  y  is  independent  of 
the  electric  field,  we  obtain  as  a  first  approximation 
of  the  empirical  formula  the  equation  of  motion  of  a 
carrier  given  by  (Be  68) 

v  =  y  e  .  3-2 

Sakai  et  al  (Sa  69)  have  shown  that  the  saturation  velo¬ 
city  for  both  holes  and  electrons  is  approximately 

7 

1.1 x  10  cm/sec  although  this  limiting  velocity  is  not 
reached  as  abruptly  for  the  holes.  Electrons  reach  this 
saturation  velocity  at  an  electric  field  of  about 
100  V/mm  in  Ge(Li)  semiconductor  material  at  77°K. 
However,  much  higher  fields  are  required  for  the  holes 
to  reach  saturation  velocity  (Sa  69) .  Cho  and  Llacer 
(Ch  72A)  did  not  observe  hole  saturation  velocity  in 
high  purity  Ge  coaxial  detectors  at  100-500  V/mm  fields. 
Hence  for  typical  fields  in  Ge(Li)  coaxial  detectors 
eq.  3-2  is  a  reasonable  approximation  for  the  holes  but 
not  for  electrons.  However,  since  we  wish  only  a 
qualitative  understanding  of  the  problem,  eq.  3-2  is 
used  for  both  the  electron  and  hole  motion  to  obtain 
conformity  on  the  final  pulse  shape  equation.  As  will 
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be  shown  later  the  form  of  the  pulse  shape  equation  is 
not  affected  by  considering  saturation  velocities. 

The  simplest  way  to  calculate  the  signal  is  to 
consider  the  energy  balance  in  the  system  (Go  66). 

The  collection  of  the  carriers  corresponds  to  an  energy 
dissipation  in  the  electrostatic  field  inside  the 
detector.  To  maintain  the  energy  stored  in  the  di¬ 
electric,  a  corresponding  charge  must  flow  in  the 
external  circuit.  The  collection  of  a  single  electron- 
hole  pair  is  depicted  in  fig.  4. 

Now  we  have  that  the  energy  W  stored  in  the 

1  2 

dielectric  is  equal  to  j  Q  /C  where  Q  is  the  charge 
in  the  detector  of  capacity  C.  Therefore, 


dW  _  Q  dQ 
dt  C  dt 
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The  field  strength  in  the  case  of  the  double-ended  co¬ 
axial  detector  is 


e  (x)  =  - 


V 


x  ln(b/a) 
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where  a  and  b  are  the  inner  and  outer  radii  of  the 
interaction  zone  of  the  coaxial  detector,  and  V  is  the 
applied  voltage.  Thus  the  drift  velocity  for  the 
electrons  is 

,  Vu 

_  dx  _  _  e 

dt  ~  x  ln(b/a) 
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Fig.  4. 


Collection  of  a  single  electron-hole  pair 
in  a  double-open-end  coaxial  Ge(Li) 
detector . 
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Integrating  this  relation,  taking  note  of  the  limits  of 
integration,  yields  for  the  electrons  and  holes, 
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respectively.  The  work  done  on  the  electron  in  time  dt 
is : 


force  •  distance  =  e(x)*e*vdt 


Therefore , 


dW  2  r  1  , 

dt  e  ye  V  {x  In  (b/a) } 


Equating  eqs.  3-3  and  3-8  we  have 
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where  I  (t)  is  the  current  due  to  the  motion  of  the 
e 

electrons.  Similarly  for  the  motion  of  the  holes 
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Substituting  A  =  y  V/ln(b/a)  and  A,  =  y,  V/ln(b/a) 

00  n  n 

into  eqs.  3-9  and  3-10,  and  using  eqs.  3-6  and  3-7,  the 
current  contributions  for  each  carrier  are: 
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By  integrating  these  expressions  with  respect  to 
charge  Q  induced  by  each  carrier  is  given  by 
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o 

At  this  point  another  approximation  is  made  in  that  the 
mobilities  for  the  electrons  and  holes  are  equal  (Ae=A^=A) 
which  is  a  good  approximation  for  Ge(Li)  semiconductor 
material  at  77°K.  Further,  if  one  assumes  that  the 
interaction  of  a  y-ray  with  the  detector  is  such  that 
the  charge  produce d) n^  e  ,  is  localised  to  a  point  which 
is  small  compared  to  the  detector  dimensions,  (the 
range  of  1  Mev  electrons  in  Germanium  is  ^0.8  mm  (St  67)), 
then  the  total  charge  induced  is  given  by 


n  e  2At  2At 

Q(t)=2  ln°b/a)  U°(l  +  — )-  in  (!-—)>  ■  3-15 

o  o 

If  saturation  velocities  for  both  the  electrons  and  holes 
had  been  used,  eq.  3-15  would  have  the  form 

n  e  ve  t  v*1  t 

Q  (t)  =  t — r&7"T-  { In  (1  +  — - )  -  In  ( 1  -  — - )  }  3-16 
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which  is  very  similar  to  that  of  eq.  3-15. 
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The  maximum  collection  time,  which  corresponds  to 
a  generation  point  near  one  of  the  two  electrode,  is 
given  by 


T 

max 


a2 ) /2A 


where 

A 


UV  , 

In  (b/a) 


Substituting  for  T  in  eq .  3-15,  one  obtains 
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Using  eq.  3-17,  pulse  shape  calculations  for  seven 
different  interaction  positions  within  a  coaxial  Ge(Li) 
detector,  outer  radius  1.54cm  and  inner  radius  0.34  cm, 
are  shown  in  fig.  5.  It  is  obvious  from  these  curves 
that  the  variations  in  pulse  shapes  and  pulse  rise-times 
are  a  major  factor  contributing  to  the  poor  timing  reso¬ 
lution  of  coaxial  Ge(Li)  detectors. 
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Fig.  5.  Calculated  pulse  shapes  for  local  interactions 
at  seven  different  positions  along  a  radius  of 
an  ideal  double-open-end  coaxial  Ge(Li)  detector. 
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The  above  pulse  shapes  are  only  a  good  approxima¬ 
tion  to  reality.  y-rays  can  interact  at  two  or  more 
locations  in  the  detector  to  produce  a  charge  which  is 
not  localised  at  one  point  but  at  several  points.  This 
would  occur  when  Compton  scattered  photons,  annihilation 
quanta,  or  bremsstrahlung  are  reabsorbed. 

Furthermore,  the  range  of  electrons  or  positrons 
above  1  MeV  is  not  negligible  within  a  1  cm  detector. 
Thus,  the  charge  produced  by  high  energy  y-rays  and  by 
y-rays  interacting  at  several  locations  tends  to  be 
distributed  in  the  detector  volume.  The  shape  of  the 
pulses  resulting  from  such  events  will  depend  upon  the 
the  spatial  distribution  of  the  charge,  including  the 
length  of  the  ionization  tracks  and  their  orientation 
with  respect  to  the  collecting  electrodes.  Thus  the 
resulting  pulses  are  likely  to  be  a  superposition  of  a 
number  of  the  waveforms  as  shown  in  fig.  5.  Neverthe¬ 
less,  the  actual  pulse  shapes  obtained  as  shown  in  fig. 3 
show  a  similar  variety  of  shapes  as  those  obtained  by 
the  above  theoretical  calculations. 

A  signal  can  be  synthesized  from  the  above 
theoretical  pulses  by  inverting  and  attenuating  it 
by  a  factor  f,  and  adding  the  result  to  the  original 
pulse  delayed  by  the  time  t^.  This  is  the  basis  of 
the  Constant  Fraction  of  Pulse  Height  Trigger  (Ge  67) 
(CFPHT)  and  the  Amplitude  Rise-Time  Compensation 
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technique  (ARC)  suggested  by  Chase  (Ch  68) .  Assuming 
that  one  has  a  fast  charge  amplifier,  the  timing  pulse 
from  the  CFPHT  is  initiated  at  the  zero-crossing  point 
of  the  synthesized  signal 
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Using  the  result,  eq.  3-17,  we  have 
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The  zero-crossing  point  of  this  function  was  found 
by  successive  approximations  on  a  computer  for  various 
values  of  f  and  T^.  Figs.  6  and  7  show  examples  of 
the  crossing  time  t  versus  the  slope  change  time  tg , 
which  is  a  function  of  the  interaction  position  in  the 
detector,  both  normalised  to  the  delay  time  t^.  For 
each  fraction  the  curves  are  made  up  of  two  parts,  one 
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Fig.  6.  Zero-crossing  time  vs.  slope  change  time,  both 

normalised  with  respect  to  the  delay  time,  t^, 

for  the  pulse  shapes  of  fig.  5,  when  applied 

to  a  CFPHT.  Delay  time  t,  =  0.05  T 

T  a  max 


10.5 


101  0- 
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Fig.  7.  Zero-crossing  time  vs.  slope  change  time,  both 

normalised  with  respect  to  tile  delay  time,  t^, 

for  the  pulse  shapes  of  fig.  5,  when  applied 

to  a  CFPHT.  Delay  time  t,  =  0.10  T 

q  max 


37 


corresponding  to  slope  change  due  to  the  motion  of  the 
holes  towards  the  inner  electrode  and  the  other  due  to 
the  motion  of  the  electrons  towards  the  outer  electrode 
of  the  detector. 

From  figs.  6  and  7 ,  the  walk  of  the  zero-cross¬ 
ing  time  with  slope  change  time  is  well  demonstrated. 

It  is  this  effect  that  leads  to  the  poor  time  resolu¬ 
tion  of  coaxial  Ge(Li)  detectors.  It  may  seem  at  first 
glance  that  the  smaller  the  fraction  the  better  the 
time  resolution.  However,  at  a  fraction  of  0.1  the 
attenuated  signal  is  very  small  and  in  many  cases  the 
same  magnitude  as  the  noise  in  the  region  of  zero-crossing 
point  of  the  synthesized  signal.  Hence  random  triggering 
of  the  zero-crossing  discriminator  occurs  leading  to  very 
poor  time  resolution. 

It  is  seen  from  figs.  6  and  7  that  after  a  certain 
value  of  t  /t , ,  the  walk  is  minimal  and  the  improved 

S  Cl 

timing  technique  to  be  described  effectively  takes  this 
region  of  these  zero-crossing  time  curves  to  obtain  the 
improved  time  resolution.  This  is  accomplished  by  tak¬ 
ing  the  curve  for  f  =  0.1  and  subtracting  it  from  the 
curve  for  f  =  0.5  to  obtain  time  difference  curves  as 
illustrated  in  figs.  8  and  9  for  the  cases  t^=  0.05  Tmax 
and  t^  =  0.10  Tmax*  Experimentally  this  can  be  done  by 
using  the  output  from  a  CFPHT  unit  set  at  f  =  0.1  to 
'START'  a  Time  to  Amplitude  Converter  (TAC),  the  'STOP' 
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Fig.  8 


.  Time  difference  curve  obtained  when  the  signals 

from  two  CFPHT's  set  at  fractions  0.1  and  0.5, 

respectively,  are  applied  to  a  TAC .  CFPRT  delay 

time,  t,  =  0.05  T 

d  max 


10.5 
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Fig.  9.  Time  difference  curve  obtained  when  the  signals 

from  two  CFPHT ' s  set  at  fractions  0.1  and  0.5, 

respectively,  are  applied  to  a  TAC .  CFPHT 

delay  time,  t,  =  0.10  T 

J  '  d  max 
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Theoretical  'shape'  curve  obtained  from  a  TAC 
using  the  time  difference  curve  of  fig.  8,  and 
considering  the  volume  contributions  at  different 
radii  of  a  coaxial  Ge(Li)  detector. 


Fig.  10. 
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Fig.  11.  Theoretical  'shape'  curve  obtained  from  a  TAC 
using  the  time  difference  curve  of  fig.  9,  and 
considering  the  volume  contributions  at  different 
radii  of  a  coaxial  Ge(Li)  detector. 


AT-* 
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for  which  is  obtained  from  a  CFPHT  set  at  f  =  0.5. 

Then  by  using  a  Timing  Single  Channel  Analyser  (SCA) 
that  portion  of  the  so-called  'shape'  spectrum  obtained 
from  the  TAC  which  shows  minimal  walk  can  be  selected 
by  setting  the  appropriate  window  in  the  SCA. 

When  one  considers  the  volume  contribution  at 
different  radii  of  a  coaxial  Ge(Li)  detector,  outer 
radius  1.54  cm  and  inner  radius  0.34  cm,  the  actual 
theoretical  time  difference  or  'shape'  curves  obtained 
from  the  TAC  are  shown  in  figs.  10  and  11.  These  curves 
are  not  quite  the  same  as  those  obtained  in  actual  ex¬ 
perimental  test  of  the  technique  for  reasons  that  will 
be  discussed  later. 

3 . 2  The  Triple  Constant  Fraction  Discriminator 

Initially,  as  described  in  a  recent  publication 
(Wh  74),  the  improved  time  resolution  of  coaxial  Ge(Li) 
detectors  was  obtained  by  the  above  technique  by  the 
use  of  three  separate  Constant  Fraction  Discriminators. 
Two  of  them  were  set  at  fractions  0.1  and  0.5,  and  their 
timing  outputs  were  used  to  'START'  and  'STOP'  the 
'shape'  TAC,  respectively. 

The  output  from  the  third,  the  setting  of  which 
at  f  =  0.3  gives  optimum  time  resolution  without  pulse 

shape  selection,  was  used  to  obtain  the  actual 
time  distribution  curve  from  a  second  TAC. 


43 


However,  the  method  suffered  from  a  few  drawbacks 
concerning  the  settings  of  the  three  discriminator 
units.  The  following  important  procedures  had  to  be 
observed: 

a)  The  threshold  level  discriminator, 

b)  The  external  shaping  delays,  and 

c)  The  zero-crossing  walk  adjusts  had  to  be 
set  to  precisely  the  same  value  on  two  of  the  discri¬ 
minator  units.  This  was  achieved  by  obtaining  an 
effective  'delta'  function  time  spectrum  from  the 
'shape'  TAC,  both  units  being  set  at  the  same  fraction, 
by  careful  adjustments  on  the  threshold  level  and  walk 
adjusts.  Even  with  this  procedure,  the  adjustments  might 
be  inadequate  if  the  fraction  switch  attenuators  were  not 
identical  in  each  unit. 

These  problems  were  overcome  by  the  development 
of  the  Triple  Constant  Fraction  Discriminator  (TCFD) 
which  incorporates  the  use  of  a  common  threshold  level 
discriminator,  external  50S7  shaping  delay,  and  a  common 
fraction  switch  attenuator.  The  circuit  diagram  is 
shown  in  fig.  12.  However,  the  fundamentals  of  the 
circuit  can  more  easily  be  understood  from  the  diagram¬ 
matic  block  circuit,  fig.  13. 

I.C.  1  is  the  common  threshold  level  comparator, 
a  'low'  output  being  obtained  between  the  arrows  when 
the  input  signal  goes  below  the  threshold  level.  In 
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Fig.  12. 


Circuit  diagram  of  the  Triple  Constant  Fraction 
Discriminator  (TCFD) . 


ro 


Fig.  13 


Schematic  block  diagram  of  the  electronic 
circuit  of  the  TCFD. 


Fig. 14.  A  photograph  of  the  Triple  Constant  Fraction 
Discriminator  module. 


I.C.  2,  3  or  4  a  comparison  is  made  between  the  attenua¬ 
ted  signal  and  the  original  signal  delayed  a  time  td  in 
the  common  external  50f2  delay.  The  I.C.  effectively 
works  as  a  zero-crossing  discriminator,  a  'low'  output 
being  obtained  between  the  arrows  at  the  crossing  points 
as  indicated. 

The  output  of  the  threshold  level  comparator 
primes  the  NOR  gate,  I.C.  8,  where  when  both  inputs  are 
in  the  'low'  state,  the  output  is  'high'.  This  signal 
is  then  amplified  in  I.C.  5,  6  or  7,  the  'high'  output 
of  which  is  then  sent  through  a  2.37  shaping  resistance 
R2 •  This  resistance  slows  up  the  rise-time  of  the  signal 
sufficiently  such  that  it  is  effectively  still  in  the 
'low'  state  with  respect  to  the  switching  level,  when 
the  'low'  output  of  I.C.  5,  6  or  7  is  applied  to  the  NOR 
gate  I.C.  9.  A  sharp  15  ns  pulse  is  thus  obtained,  the 
length  of  which  is  determined  by  the  resistance  R£ ,  and 
then  inverted  for  use  as  an  output  signal.  Fig.  14  shows 
a  photograph  of  the  Triple  Constant  Fraction  Discrimina¬ 
tor  module. 

3 . 3  Experimental  Tests  of  Improved  Timing  Technique 

For  bench  testing  experiments,  the  TCFD  unit  was 
placed  in  the  electronic  circuit  shown  diagrammatically 
in  fig.  15.  Timing  coincidence  was  obtained  by  placing 


ADC  Q  ADC  P  COIN.  ADC  R 

(ENERGY)  (SHAPE)  (TIME) 


Fig.  15.  Block  diagram  of  the  bench-testing 
timing  system. 
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a  Co  or  Na  source  between  a  23  c.c.  coaxial  Ge(Li) 
detector  and  a  Naton  scintillator  attached  to  a  C3/1024 
fast  multiplier  tube.  One  of  the  signals  from  the 
Ge(Li)  detector  is  fed  through  a  linear  amplifier  and 
into  ADC  Q.  The  other  is  sent  through  a  Timing  Filter 
Amplifier  (TFA)  and  into  the  TCFD  unit.  This  unit 
has  three  timing  outputs  each  of  which  can  be  individually 
set  at  any  fraction  0.5  (K) ,  0.3  (L)  and  0.1  (M)  by 

using  a  common  fraction  switch  attenuator  as  shown  in 
figs.  12  and  13.  The  output  at  fraction  0.3  is  used  to 
'STOP'  the  time  TAC .  The  'START'  signal  is  obtained 
from  the  scintillator.  The  output  of  the  time  TAC  is 
then  fed  into  ADC  R.  The  other  two  outputs  of  the  TCFD 
are  set  at  fractions  0.1  and  0.5,  'START'  and  'STOP' 
the  'shape'  TAC,  respectively,  the  output  of  which  is 
fed  into  ADC  P.  ADC's  P,  Q  and  R  require  a  coincidence 
signal  which  is  provided  by  the  three  timing  SCA's  and 
the  Universal  Coincidence  box.  The  Delay  Amplifiers 
which  are  used  to  obtain  a  coincidence  between  the  sig¬ 
nals  fed  into  ADC's  P,  Q  and  R  and  the  signal  from  the 
coincidence  box  are  not  shown  in  fig.  15. 

Fig.  16  shows  the  total  time  and  'shape'  spectra 

fi  f) 

for  the  Co  source  as  obtained  from  ADC  R  and  ADC  P, 
respectively.  The  time  spectrum  obtained  is  for  all 
energies  and  with  no  'shape'  selection.  The  FWHM  and 
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Fig.  16  a)  The  total  time  spectrum  obtained  from  ADC 
R  and 

b)  The  total  'shape'  spectrum  obtained  from 
ADC  P  for  a  C0Co  source. 
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Sorted  time  spectra  for  seven  equal  windows 
set  on  the  'shape'  spectrum  between  the 
arrows  as  indicated  in  fig.  16b),  the  system 
being  gated  by  an  energy  window  on  the  1.332 
MeV  y-rays  of  the  ^Co  source. 


Fig.  17. 


FWTM  are  2.79  ns  and  10.80  ns,  respectively.  Using  a 
two  dimensional  sort  routine,  the  total  time  spectrum 
was  sorted  into  bins  according  to  windows  set  on  the 
energy  and  'shape'  spectra.  For  the  60Co  source  an 
energy  window  was  set  on  the  1.33  MeV  peak  and  14  equal 
size  windows  were  set  on  the  'shape'  spectrum.  Fig.  17 
shows  seven  of  the  sorted  time  spectra  for  the  region 
between  the  arrows  as  indicated  in  fig.  16.  The  values 
of  the  FWHM  and  FWTM  of  the  best  timing  bin,  bin  5,  are 
1.06  ns  and  2.16  ns,  respectively. 

As  seen  in  fig.  17  there  is  a  significant  im- 

6  0 

provement  in  the  time  resolution  obtained  for  the  Co 
source,  especially  in  the  FWTM.  The  deterioration  in 
the  early  slope  of  the  time  distribution  curve  which  is 
a  feature  of  the  ARC  timing  technique  is  completely 
eliminated  by  our  pulse  shape  discrimination  technique, 
which  is  seen  to  separate  this  effect  mostly  into  bins 
1 ,  2  and  3 . 

22 

Similar  results  were  obtained  for  a  Na  source. 
Fig.  18  shows  the  total  time  and  'shape'  spectra  as 
obtained  from  ADC  R  and  ADC  P,  respectively.  In  this 
case  no  'shape'  or  energy  selection  was  imposed  on  the 
system.  Again  the  total  time  spectrum  was  sorted 
according  to  an  energy  window  set  on  the  0.511  meV  peak 
and  14  equal  size  windows  on  the  'shape'  spectrum. 

Fig.  19  shows  seven  of  the  sorted  time  spectra  for  the 
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Fig.  18  a)  The  total  time  spectrum  obtained  from  ADC 
R  and 

b)  The  total  'shape'  spectrum  obtained  from 
ADC  P  for  a  22Na  source. 
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Fig.  19.  Sorted  time  spectra  for  seven  equal  windows 
set  on  the  'shape'  spectrum  between  the 
arrows  as  indicated  in  fig.  18b) ,  the 
system  being  gated  by  an  energy  window  on 
the  0.511  MeV  annihilation  y-rays  of  the 
source. 


region  between  the  arrows  as  indicated  in  fig.  18. 

Again  there  is  a  marked  improvement  in  the  resolution 
of  the  time  spectra.  Values  for  the  FWHM  and  FWTM  of 
the  best  timing  bin,  bin  5,  are  2.29  ns  and  4.60  ns, 
respectively . 

A  more  impressive  way  of  showing  the  effectivenes 

of  this  pulse  shape  selection  technique  is  a  plot  of  the 

FWHM  and  FWTM  against  the  percentage  of  the  total  number 

of  events  as  one  adds  to  the  best  timing  bin,  the 

sorted  bins  from  windows  set  to  the  left  and  right 

across  the  shape  spectrum.  Fig.  20  shows  the  results 

for  the  1.33  MeV  y-ray  of  ^Co.  Similar  results  are 

22 

shown  in  fig.  21  for  the  0.511  MeV  peak  using  a  Na 
source . 

From  these  curves  it  is  seen  that  by  a  selection 
of  about  50  to  60%  of  the  total  number  of  events  a  sig¬ 
nificant  improvement  in  the  time  resolution  can  be 
obtained.  This  is  one  of  the  greatest  advantages  of 
the  technique  that  improved  timing  resolution  can  be 
obtained  without  a  serious  loss  of  counting  efficiency. 
In  a  recent  experimental  investigation  (Hu  74) 

the  technique  was  applied  to  the  measurement  of  life- 
35  3  7  39 

times  in  Cl,  Cl  and  K  using  a  pulsed  proton  beam 
as  a  method  of  exciting  the  states  of  interest.  The 
improvement  in  the  time  resolution  was  such  that  life¬ 
times  in  the  50  ps  region  and  below  were  measured  and 
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PERCENTAGE 

Fig.  20.  A  plot  of  the  FWHM  and  FWTM  against  the  per¬ 
centage  of  the  total  number  of  events  as 
sorted  time  bins  are  added  to  the  best  timing 
bin  from  windows  set  to  left  and  right  across 
the  shape  spectrum.  The  system  was  gated  on 
the  1.332  MeV  y-rays  of  the  60Co  source. 
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Fig.  21.  A  plot  of  the  FWHM  and  FWTM  against  the  per¬ 
centage  of  the  total  number  of  events  as 
sorted  time  bins  are  added  to  the  best  timing 
bin  from  windows  set  to  left  and  right  across 
the  shape  spectrum.  The  system  was  gated  on 
the  0.511  MeV  peak  of  a  “^Na  source. 
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time  resolutions  of  800  ps  and  650  ps  FWHM  were  obtained 
in  the  energy  range  2.20  MeV  to  3.16  MeV,  respectively. 

In  another  application  to  be  discussed  in  a  later 

chapter,  the  technique  has  been  used  to  study  lifetimes 
45  45 

in  Ti  and  Sc.  The  data  from  this  run  has  been  used 
to  show  the  change  in  FWHM  of  the  prompt  time  distribu¬ 
tion  curves  as  a  function  of  the  inverse  of  the 
associated  detected  energy  as  plotted  in  fig.  22.  The 
linearity  of  the  curve  indicates  the  way  in  which  elec¬ 
tronic  noise  is  limiting  the  time  resolution.  Also 
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shown  m  fig.  22  is  some  data  from  the  Cl,  Cl,  K 
experiment.  It  is  seen  that  this  data  falls  on  a 
separate  line.  It  is  generally  found  that  the  exact 
position  of  the  line  is  dependent  upon  the  tuning  of  the 
electronic  system  and  the  digital  window  set  on  the 
' shape '  spectrum . 

Two  interesting  electronic  effects  have  been 
found  to  occur  with  the  Triple  Constant  Fraction  Dis¬ 
criminator.  The  first  effect  is  that  which  leads  to 
the  'delta'  function  peak  in  the  'shape'  spectra  as 
shown  in  figs.  16  and  18.  This  peak  was  found  to  be 
related  to  pulses  from  the  Ge(Li)  detector  which  have 
poor  rise-time  characteristics  and  secondly  to  those 
pulses  of  very  small  amplitude  as  illustrated  in  fig.  23. 
In  these  cases  the  pulses  from  the  zero-crossing  compara¬ 
tors  prime  the  NOR  gates  I.E.  8A,  B  and  C  before  the 
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Fig.  22.  FWHM  of  the  prompt  time  distributions  plotted  as 
a  function  of  the  inverse  of  the  associated 
detected  energy.  The  circles  represent  data 
from  the  ^Sc  experimental  run  while  the  crosses 

c  o  c:  on  on 

are  some  data  for  the  Cl,  Cl,  K  experimental 
run . 
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F i g .  23.  Illustration  of  some  of  the  ways  in  which 
pulses  from  the  Ge(Li)  detector  are  mis¬ 
timed.  Pulses  of  the  form  A  and  B  trigger 
the  NOR  gate  I.C.8  before  T  ,  the  arrival 
time  for  signals  from  the  zero-crossing 
comparator,  and  are  correctly  timed.  Pulses 
of  the  form  C  and  D  occur  later  than  and 
are  thus  mistimed. 
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arrival  of  the  threshold  level  pulse.  Hence  the  timing 
is  determined  by  the  latter  pulse  and  since  this  pulse 
is  common  to  all  three  NOR  gates  a  'delta'  function 
peak  is  obtained  in  the  shape  spectrum. 

Fig.  24  shows  the  results  of  setting  a  window 

on  the  'delta'  function  peak  for  the  total  time  spec- 

22 

trum  obtained  from  a  Na  source.  It  is  seen  that  those 
pulses  that  give  rise  to  this  effect  correspond  to 
pulses  with  poor  timing  characteristics.  Hence  with 
the  pulse  shape  selection  technique  to  improve  the  time 
resolution,  these  events  are  eliminated. 

The  second  effect  is  seen  in  the  small  second 
hump  in  figs.  16  and  18.  In  general  the  zero-crossing 
comparators  are  continuously  being  triggered  by  noise 
and  real  events.  Therefore,  if  the  threshold  level 
pulse  arrives  before  the  attenuated  signal  has  attained 
a  greater  amplitude  than  the  noise,  the  NOR  gate  will 
then  fire  on  the  threshold  level  input.  As  might  be 
expected  this  effect  is  mostly  restricted  to  the  0.1 
fraction  output  and  occurs  with  slow  rise-time  pulses 
which  have  poor  timing  characteristics.  It  can  also 
occur  with  very  fast  rising  pulses.  The  effect  can  be 
diminished  by  increasing  the  size  of  the  resistance  R^ , 
thereby  slowing  down  the  rise-time  of  the  threshold 
comparator  output  and  increasing  the  possibility  of  the 
attenuated  signal  rising  above  the  noise  level. 
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Fig.  24.  a)  A  logarithmic  plot  of  the  total  time 
spectrum  obtained  for  a  ^Na  source, 
b)  The  same  but  with  the  system  gated 
by  a  window  on  the  'delta'  function 
peak  in  the  'shape'  spectrum. 

The  two  graphs  are  displaced  by  half  a 
division  to  ensure  clarity. 
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Fig.25A  shows  the  effect  on  this  second  hump  when  the 
resistance  is  set  equal  to  zero  and  when  equals 
870  ft. 

A  final  illustration  of  the  excellent  timing 
characteristics  of  the  Triple  Constant  Fraction  Dis¬ 
criminator  is  shown  in  fig.  25B.  in  this  case  the  two 
outputs  of  the  TCFD  which  are  fed  into  the  'shape'  TAC 
are  set  both  at  the  same  fractions.  The  'delta' 
functions  (right  peaks  in  each  case)  obtained,  over 
the  range  of  fraction  settings  from  0.1  to  0.5,  are 
the  desired  results  of  the  common  threshold  level  dis¬ 
criminator,  external  50  ft  delay  and  the  common  fraction 
switch  attenuator.  The  left  hand  'delta'  function  in 
each  case  is  the  same  as  that  discussed  two  paragraphs 


earlier . 
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Fig.  25.  a)  The  total  'shape'  spectra  when  resistance 
R]_  =  0  and  R]_  =  870  fi. 
b)  Spectra  obtained  from  the  ' shape 'TAC  when 
the  two  outputs  of  the  TCFD  are  set  at  the 
same  fractions. 


CHAPTER  4 


THE  BEAM  PULSING , TARGET  ROTATING  SYSTEMS  AND 
EXPERIMENTAL  METHOD 


4 . 1  The  Mobley  Compression  Pulsed  Beam  System 

In  the  University  of  Alberta's  7.5  Mv  van  de 
Graaff  the  pulsed  beam  is  obtained  by  using  both  pre¬ 
acceleration  pulsing  and  post-acceleration  compression 
as  shown  schematically  in  fig.  26.  An  outline  of  the 
system  will  be  given  here  due  to  its  importance  in  the 
experimental  technique.  However,  a  more  detailed  des¬ 
cription  can  be  obtained  in  the  Ph.D.  thesis  of  Davies 
(Da  66 ) . 

An  ion  source  located  within  the  high  voltage 
terminal  is  used  to  ionise  the  required  gas  by  a  125  MHz 
source  oscillator  and  a  d.c.  ion  beam  is  extracted  by 
means  of  a  probe  electric  field.  An  Einzel  Lens  and 
Terminal  Analysis  Magnet  are  used  to  focus  the  beam 
which  can  then  be  pulsed  by  two  sets  of  R.F.  deflection 
plates  set  at  right  angles  to  each  other  and  driven  90° 
out  of  phase  by  a  1  MHz  oscillator.  These  R.F.  deflec¬ 
tion  plates  sweep  the  beam  in  an  elliptical  path  over 
a  1/16  inch  aperture,  producing  a  beam  pulse  of  10  ns 
width.  This  beam  pulse  is  then  accelerated  and  further 


65 


. 


■ 


66 


c o 
< 
CD 


I-6S 

CL 


h- 

UJ 

2 

CD 

< 

2 

UJ 

O 

cr 

3 

O 

CO 

□ 


UJ 

o 

cr 

3 

O 

CO 


□ 


cr 

o 

H 

< 

_J 

_J 

CD 

CO 

O 


CO 

2 

UJ 


UJ 

N 

2 

UJ 


►— 

UJ 

2 

CD 

< 

2 

CO 

CO 

>■ 

_J 

< 

2 

< 


2 

cr 

uj 

K- 


<  O 


2 

O 

1- 

o 

UJ  w 
_J  UJ 
U_  J— 
UJ  < 
Q  _J 

a 

ill 

cr 


UJ 

cr 

3 

ro 

CD 

2 

H- 

CO 

Kl 

h- 

cr 

3 

>- 

UJ 

UJ 

O 

2 

CL 

o 

< 

CD 

< 

u_ 

2 

< 

< 

2 

Fig.  26.  Schematic  diagram  of  the  beam  transport 

and  Mobley  Compression  pulsed  beam  system. 


focused  by  Focus  3  and  analysed  by  a  90°  magnet  to 
select  the  desired  energy.  The  beam  then  passes 
through  a  pair  of  energy  stabilising  slits  57  inches 
away  from  the  analysing  magnet.  These  slits  feed  the 
corona  current  central  system  which  stabilises  the 
energy . 

A  magnetic  quadrupole  lens  then  focuses  the  beam 
into  the  first  stage  of  the  Mobley  compression  system. 
This  stage  of  the  system  is  a  pair  of  deflection  plates 
fed  with  a  phase  locked  R.F.  oscillator  of  10  MHz. 

Phase  information  is  supplied  to  the  oscillator  from  a 
cylindrical  pick-off  in  the  beam  tube  just  before  the 
R.F.  plates. 

The  second  stage  is  a  90°,  75  inch  magnet  whose 
focal  points  are  at  the  target  and  the  R.F.  plates. 

The  beam  is  swept  by  the  R.F.  field  such  that  the  par¬ 
ticles  leaving  the  deflector  plates  first,  take  the 
longest  paths  through  the  magnet  while  those  leaving 
last  take  the  shortest  paths.  If  the  phase  and  ampli¬ 
tude  of  the  deflector  are  properly  adjusted,  all  parti¬ 
cles  will  arrive  at  the  target  at  almost  the  same  time. 
A  final  beam  pulse  with  width  of  0.4  ns  can  be  obtained 
at  the  target. 

4 . 2  The  Rotating  Target  System 


As  will  be  seen  later, 


the  centroids  of  the 


time  distribution  peaks  vary  with  the  energy  of  the 
detected  y-ray.  This  may  not  pose  a  problem  if  the 
so-called  prompt  time  distribution  centroid  curve  is 
defined  accurately  around  the  energy  region  of  interest 
by  a  sufficient  number  of  centroid  points  obtained  from 
y-rays  whose  lifetimes  are  much  shorter  than  the  time 
resolution  of  the  system  (less  than,  say,  1  ps)  .  How¬ 
ever,  in  most  cases  the  isotope  under  study  does  not 
have  a  sufficient  number  of  'prompt'  deexcitation  y-ray 
in  the  required  energy  region.  It  is  essential  to  the 
accuracy  of  the  results  that  this  problem  be  overcome. 

A  simple  solution  is  the  use  of  two  targets,  the  one 
producing  the  energy  level  whose  lifetime  is  to  be  mea¬ 
sured  and  the  second,  supplementary  'prompt'  reference 
y-rays . 

In  most  previous  delayed  coincidence  experiments  the 
prompt  and  delayed  time  distribution  curves  were  ob¬ 
tained  in  alternate  experimental  runs.  The  errors  due 
to  electronic  drifts  and  shifts  were  minimised  by 
frequent  alternation  of  the  measurements  with  the  two 
sources.  This  method  is  obviously  cumbersome  and  thus 
it  is  very  desirable  to  obtain  the  prompt  and  delayed 
time  distributions  simultaneously  in  order  to  avoid 
errors  due  to  electronic  drifts  and  shifts  and  changes 
in  pulsed  beam  conditions  from  one  run  to  the  next. 
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Another  important  point  to  be  discussed  in 
Chapter  5  is  the  observed  centroid  shift  for  a  given 
energy  between  the  time  peaks  obtained  from  events  in 
which  the  full  energy  of  the  'prompt'  y-rays  is  depo¬ 
sited  in  the  Ge(Li)  detector,  and  those  obtained  from 
Compton  events  in  which  part  of  the  energy  of  the 
incident  'prompt'  y-rays  is  deposited.  Thus  the  peak 
to  background  ratio  is  very  important  in  the  determina¬ 
tion  of  the  true  centroid  of  the  background  subtracted 
time  distributions  and  the  associated  error  in  this 
result.  Background  contributions  can  be  reduced  by  the 
use  of  self-supporting  thin  targets.  The  possibility  of 
target  fracture  by  the  incident  beam  is  alleviated  by  the 
rotation  of  the  target,  effectively  increasing  the  tar¬ 
get  area  to  the  incident  beam. 

With  these  points  in  mind  it  became  necessary  to 
design  a  twin  target  rotating  system.  Two  photographs 
of  the  rotating  target  system  taken  from  different 
angles  are  shown  in  figs.  27  and  28.  Rotation  in  the 
vacuum  target  chamber  is  provided  by  a  magnetic  coupl¬ 
ing  capable  of  550  r.p.m.  which  is  driven  by  a  1/8  h.p. 
variable  speed  motor.  A  system  of  geared  pulleys  then 
increases  the  maximum  r.p.m.  of  the  target  assembly  to 
7000,  the  final  drive  being  accomplished  by  the  fric¬ 
tion  of  an  O-ring  pulley.  The  target  rotator  fits  into 
the  inner  race  of  a  5.08  cm  bearing  which  can  be 
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Fig. 27.  A  photograph  of  the  Rotating  Target  System. 
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Fig. 28.  A  photograph  of  the  Rotating  Target  System. 


screwed  in  or  out  of  the  main  assembly  mount  to  provide 
the  necessary  and  optimum  friction  for  maximum  rota¬ 
tion.  It  was  found  that  best  performance  was  obtained 
when  the  bearings  were  lubricated  with  a  small  amount 
of  Molybdenum  Disulphide  powder.  The  target  rotator, 
bearings  and  bearing  mount  were  made  of  the  same 
material,  steel, to  overcome  difficulties  due  to  heat 
expansion  which  might  cause  the  inner  and  outer  races 
of  the  bearings  to  seize.  The  system  is  enclosed  in  a 
target  chamber  which  has  a  6.35  cm  diameter  window  of 
0.12  mm  steel  so  designed  that  targets  can  be  viewed 
from  50°  to  130°  with  respect  to  the  beam  axis. 

4 . 3  Target  Discriminator 

When  using  two  targets,  each  rotating  in  and  out 
of  the  beam,  it  is  useful  to  store  the  information  from 
each  target  separately  in  different  regions  of  the  on¬ 
line  computer.  A  separation  of  the  energy  spectra  from 
each  target  means  that  the  peak  to  background  ratio  is 
not  adversely  affected  by  the  presence  of  the  second 
target.  If  the  background  should  also  contain  the 
Compton  tail  of  a  y-ray  transition  from  a  state  whose 
lifetime  is  several  nanoseconds  then  the  separation  of 
the  energy  spectra  eliminates  the  problematic  effects 
of  this  lifetime  on  the  time  distribution  peaks 
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corresponding  to  energy  windows  set  on  the  other  spectrum. 

To  produce  this  separation  of  information  from 
each  target  two  slot  apertures  are  machined  in  the  tar¬ 
get  rotator  diametrically  opposite  at  different  radii. 

The  targets  are  made  semicircular  and  screw  onto  the 
target  rotator.  In  conjunction  with  these  slots,  a 
light  source  and  two  photodiodes  are  arranged  such  that 
as  the  targets  rotate  light  can  pass  through  either  of 
the  slots  alternately  and  be  detected  by  the  photodiodes. 
The  beam  is  collimated  by  a  system  of  movable  slits  up¬ 
stream  from  the  target  and  is  only  incident  on  one 
target  at  a  time  since  the  axis  of  rotation  of  the 
targets  is  0.63  cm  off  the  beam  axis. 

The  outputs  from  the  two  photodiodes  pass  through 
a  kovar  connection  from  the  target  vacuum  chamber  and 
are  applied  to  the  so-called  Target  Discriminator  through 
high  input  impedance  F.E.T.  circuits  since  the  photo¬ 
diodes  have  a  high  output  impedance.  Fig.  29  shows  the 
electronic  circuit  of  the  Target  Discriminator  and 
fig.  30  a  photograph  of  the  actual  NIM  module.  The 
fundamentals  of  the  Target  Discriminator  are  illustrated 
in  fig.  31.  Basically  the  signals  from  the  photodiodes 
are  amplified,  inverted  and  sent  into  a  flip-flop  cir¬ 
cuit.  The  outputs  from  this  and  the  original  inverted 
pulses  are  then  fed  into  two  separate  AND  gates  as 
illustrated  in  fig.  31.  The  final  outputs  are  high, 


Fig.  29.  Circuit  diagram  of  the  Target  Discriminator 
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Fig . 30 . 


A  photograph  of  the  Target  Discriminator  Module. 
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Fig.  31.  A  schematic  diagram  of  the  pulse  shape  logic 
in  the  Target  Discriminator  circuit. 
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amplitude  5  V,  or  low  when  the  targets  are  in  or  out 
of  the  beam,  respectively.  This  high  output  has  an 
amplitude  sufficient  to  trigger  a  Universal  Coincidence 
box  and  is  thus  compatible  with  other  logic  signals 
used  for  coincidence  purposes.  It  is  noticed  that  both 
signals  are  low  when  the  beam  is  switching  targets  and 
thus  a  complete  separation  of  the  energy  spectra  from 
each  target  can  be  achieved. 

4 . 4  Experimental  Electronics  and  Method 

Fig.  32  shows  the  block  diagram  of  the  experimen¬ 
tal  electronic  arrangement.  The  deexcitation  y-rays, 
detected  using  a  23  c.c.  coaxial  Ge(Li)  detector,  produce 
two  signals,  one  of  which  is  sent  through  a  linear 
amplifier  and  into  ADC  P  consisting  of  two  regions, 
each  2048  channels.  The  other  output  is  fed  into  the 
Triple  Constant  Fraction  Discriminator,  TCFD,  the 
output  of  0.3  fraction  being  used  to  'START'  the  time 
TAC.  The  'STOP'  signal  is  provided  by  a  capacitive 
beam  pick-off.  The  output  of  the  time  TAC  is  then  fed 
into  ADC  R,  size  1024  channels. 

The  other  two  outputs  of  the  TCFD  at  fractions 
0.1  and  0.5,  respectively,  'START'  and  'STOP'  the 
'shape'  TAC  to  produce  a  'shape'  spectrum  in  ADC  Q, 
size  1024  channels.  The  on-line  computer,  a  Honeywell 
516,  requires  a  coincidence  signal  for  all  three  ADC's 
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Fig.  32.  Block  diagram  of  the  experimental  electronic 
arrangement . 
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which  is  provided  by  the  three  timing  SCA's  and  a 
Universal  Coincidence  box.  The  coincidence  box  is 
set  to  require  4  out  of  5  coincidence  signals  since 
the  Target  Discriminator  provides  alternate  signals 
when  the  beam  is  on  one  target  or  the  other.  The 
signal  from  the  second  coincidence  box  is  used  to 
rout  ADC  P  so  that  the  energy  spectrum  from  each  target 
is  stored  in  the  two  separate  regions  of  ADC  P.  Fig. 31 
does  not  show  the  several  Delay  Amplifiers  and  Delay 
boxes  necessary  to  establish  a  coincidence  between  the 
final  coincidence  signal  and  those  signals  fed  into  ADC 
P,  Q  and  R. 

The  acquisition  of  data  into  the  various  regions 
of  memory  of  the  on-line  Honeywell  516  computer  is 
accomplished  by  means  of  a  two  dimensional  sort  routine 
written  by  Dr.  D . A .  Hutcheon.  Using  this  routine  the 
total  time  spectrum  of  ADC  R  is  sorted  internally  by  the 
computer  according  to  digital  windows  set  on  the  energy 
spectra  of  ADC  P  and  on  that  region  of  ADC  Q,  the  'shape' 
spectrum,  which  produces  improved  time  resolution.  A 
maximum  of  32  digital  window  combinations  can  be  set  on 
the  spectra  of  ADC  P  and  Q,  and  the  final  sorted  time 
spectra  corresponding  to  each  window  are  stored  in  region 
X  of  the  computer,  a  maximum  size  of  16384  channels. 

During  an  experimental  run  windows  are  set  on 
the  full  energy  y-ray  peaks  to  be  investigated,  including 


those  peaks  to  act  as  'prompt'  references,  and  on  the 
backgrounds  above  and  below  these  peaks.  The  centroids 
of  the  resulting  sorted  time  distribution  curves  were 
obtained  by  setting  'begin'  and  'end'  channels  at 
approximately  5%  and  7^%  of  the  peak  height  of  each 
curve.  An  average  of  these  two  results  was  then  cal¬ 
culated  with  its  associated  error.  This  average  should 
be  a  more  accurate  value  of  the  centroid  of  the  curve 
because  setting  the  'begin'  and  'end'  channels  is  only 
approximate  since  the  statistics  are  poor  in  the  outer 
less  well-defined  parts  of  the  curves. 

The  centroids  obtained  corresponding  to  energy 
windows  set  on  the  background  were  then  plotted  as 
illustrated  by  crosses  in  fig.  33.  A  background  time 
centroid  curve  is  drawn  through  these  points  and  the 
values  of  the  background  centroids  at  the  energies  of 
the  y-ray  peaks  can  be  interpolated  so  that  background 
subtractions  can  be  executed.  The  resulting  centroids 
are  then  plotted  and  the  'prompt'  reference  centroids 
are  fitted  with  the  function  described  below. 

In  most  CFPHT  units,  zero-crossing  detection  is 
based  upon  the  use  of  a  tunnel  diode  which  does  not 
switch  instantaneously  but  requires  an  accumulation 
of  a  certain  amount  of  charge  to  activate  the  diode. 
This  charge  is  represented  by  the  area  under  the  curve 
as  illustrated  in  fig.  34,  and  leads  to  the  timing 
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Fig.  33. 


Centroids  of  the  time  distributions  vs.  energy 
for  the  ^^Ni-l^lTa  experimental  run.  Crosses 
represent  centroids  from  energy  windows  set  on 
the  Compton  continuum  and  circles  the  centroids 
for  the  time  distributions  of  'prompt*  full 
energy  peaks. 
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Illustration  of  the  timing  error.  At, 
introduced  by  the  requirement  that  a 
small  amount  of  charge  is  needed  to 
activate  the  cross-over  diode. 


Fig.  34. 
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error,  At.  This  is  the  main  factor  leading  to  the  walk 
in  the  prompt  centroid  time  curve.  By  considering 
pulses  of  different  amplitudes  i.e.  energy  E,  and  the 
same  amount  of  charge  to  activate  the  diode,  different 
triggering  times  occur  as  shown  in  fig.  34.  A  simple 
relation  can  be  calculated  between  At  and  E  using 
linearly  rising  pulses  with  equal  rise-times: 

At  =  —  .  4-1 

/E 

Hence,  Constant  Fraction  triggering  occurs  at  a  time 

t  =  a  +  —  .  4-2 

/E 

This  relation  was  used  to  fit  the  prompt  centroid  points 
with  surprising  success.  An  extra  term/C  loge  E,was 
added  to  improve  the  fit  at  lower  energies  since 
the  walk  of  the  prompt  time  centroid  curve  is 
more  pronounced  in  this  energy  region. 

However,  in  most  cases  the  parameter  C  had  a  small 
value.  Any  centroid  which  is  shifted  from  this  curve 
corresponds  to  a  y-ray  with  a  measurable  lifetime.  The 
time  spectrum  ADC  R  was  calibrated  using  a  precision 
Time  Calibrator,  ORTEC  model  462. 

One  of  the  greatest  advantages  of  the  above  tech¬ 
nique  is  the  simultaneous  acquisition  of  data  for  both 
the  prompt  and  delayed  time  distribution  curves.  Errors 
due  to  electronic  drifts  and  shifts  which  may  occur 
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between  alternate  experimental  runs  as  was  the  case  in 
previous  delayed-coincidence  measurements,  are  thus 
avoided.  Also  in  previous  measurements  background 
subtractions  were  inaccurate  due  to  the  use  of  scintilla¬ 
tor  detectors  with  their  poor  energy  resolution.  Our 
technique  goes  a  long  way  to  solving  this  problem  by  the 
use  of  the  target  rotator  system.  Separation  of  the 
spectra  from  the  two  targets  by  the  Target  Discriminator 
module  prevents  increasing  the  background  to  peak  ratio 
which  is  a  major  factor  in  producing  errors  in  background 
subtractions.  Due  to  the  extremely  good  energy  resolu¬ 
tion  of  Ge(Li)  detectors,  digital  windows  can  be  placed 
on  the  background  continuum  above  and  below  the  y-ray 
peak  to  within  a  few  kev  of  the  peak,  to  ascertain  the 
true  background  contribution  to  the  time  distribution 
curves  corresponding  to  energy  windows  on  the  y-ray  peaks. 

Another  advantage  is  that  several  levels  can  be 
investigated  for  lifetimes  simultaneously  and  hence  the 
influence  on  the  lifetimes  of  cascading  y-rays  to  the 
levels  of  interest  can  be  studied.  This  is  especially 
important  if  one  of  the  states  cascading  to  the  level 
of  interest  has  a  measurable  lifetime.  Coincidences 
from  other  cascades  has  greatly  distorted  measurements 
in  previous  delayed-coincidence  techniques. 

Another  plus  factor  of  the  technique  is  that 
centroid  shifts  are  measured  relative  to  the  prompt 
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time  centroid  curve  at  the  appropriate  energy.  This 
is  important  due  to  the  walk  in  the  prompt  time  cen¬ 
troid  curve.  An  ultimate  refinement  of  the  technique 
would  be  the  development  of  the  timing  electronics 
which  produce  a  walk  free  prompt  time . centroid  curve. 

Simms  et  al .  (Si  61)  and  Kim  and  Milner  (Ki  71) 
also  recognized  the  necessity  of  obtaining  the  prompt 
and  delayed  time  distribution  curves  simultaneously. 

Simms  et  al. obtained  the  desired  result  by  a  triple 
coincidence  technique  which  included  g-y  and  y-y 
coincidences  using  plastic  scintillators  and  a  Nal 
detector.  However,  the  method  is  limited  by  the  number 
of  states  that  can  be  excited  by  g  and  y-ray  transitions. 
Another  important  point  which  was  not  considered  is  the 
walk  of  the  centroid  of  the  prompt  time  distributions 
with  the  energy  of  the  detected  radiation  although  this 
effect  is  less  significant  with  plastic  scintillators 
than  with  Ge(Li)  detectors.  Kim  and  Milner,  on  the 
other  hand,  using  a  pulsed  beam  obtained  the  prompt 
and  delayed  time  distributions  simultaneously  by  having 
the  energy  and  time  distributions  accumulate  in  a  two- 
parameter  analyzer  operated  in  the  512x32  channel  mode; 
512  channels  for  energy  analysis  and  32  channels  for 
the  time  analysis.  Supplementary  prompt  reference  y- 
rays  were  obtained  from  a  second  target  in  a  separate 
experimental  run.  The  main  problem  of  the  two  parameter 
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analysis  is  that  the  time  dispersion  [0.567  ns/channel] 
is  large  compared  to  lifetimes  of  the  order  of  10-11 
secs  and  hence  the  resolution  of  the  time  distributions 
must  be  good  if  centroid  shifts  of  the  order  of  a  tenth 
of  a  channel  are  to  be  measured  with  accuracy.  The 
energy  resolution  of  the  detector  must  also  be  excellent 
to  resolve  y-ray  peaks  differing  in  energy  by  a  few 
keV  since  there  are  only  512  channels  for  energy  analy¬ 
sis.  The  present  technique  overcomes  these  problems  by 
allowing  the  experimenter  more  freedom  in  determining 
the  size  of  the  individual  regions  of  spectra.  Also 
the  ability  to  choose  only  those  y-ray  peaks  of  interest 
and  to  be  able  to  sort  the  data  with  respect  to  another 
parameter,  are  additional  advantages  of  the  present 
technique . 

There  are  two  major  limitations  to  the  technique. 
The  first  is  the  availability  of  suitable  'prompt' 
reference  y-rays.  The  second  is  the  peak  to  background 
ratio.  It  was  found  that  when  this  ratio  decreased  below 
1:2  large  uncertainties  were  introduced  in  background 
subtractions.  Hence  only  y-rays  with  reasonable  inten¬ 
sities  can  be  studied. 

‘i 

Another  background  problem  is  overcome  by  setting 
the  detector-target  distance  to  about  1  metre  so  that 
by  time  of  flight  neutrons  produced  in  nuclear  reaction 
at  the  target  can  be  differentiated  from  y-rays.  Other- 


wise  the  detection  of  the  neutrons  leads  to  a  time  peak 
which  falls  close  to  or  under  the  time  distribution 
curve  obtained  from  y-ray  detection.  This  situation 
is  obviously  unacceptable.  Although  the  large  detector- 
target  distance  reduces  the  count  rate  considerably, 
most  experimental  runs  were  less  than  3  hours. 


CHAPTER  5 


THE  FULL  ENERGY  -  COMPTON  DILEMMA 


In  all  the  experimental  runs  performed  an  interest¬ 
ing  feature  has  been  consistently  observed.  It  is  found 
that  for  a  given  pulse  height  there  is  a  centroid  shift 
between  the  time  peaks  obtained  from  events  in  which  the 
full  energy  of  the  'prompt'  y-rays  is  deposited  in  the 
Ge(Li)  detector,  and  those  obtained  from  Compton  scatter¬ 
ing  of  higher  energy  'prompt'  y-rays.  An  ideal  'prompt' 
y-ray  is  that  which  originates  from  a  level  whose  life¬ 
time  is  much  less  than,  say,  5  ps .  An  example  of  this 
effect  is  shown  in  fig.  33  from  experimental  data  ob¬ 
tained  in  the  study  of  lifetimes  in  ^Cu.  For  future 
reference  we  define  the  background  and  prompt  centroid 
time  curves  as  the  'best  fit'  curves  to  the  centroid 
points  from  energy  windows  set  on  the  Compton  background 
and  'prompt'  y-ray  peaks,  respectively. 

It  is  noted  that  for  very  low  energies,  those 
less  than  80  kev,  the  background  centroid  time  curve 
crosses  the  prompt  centroid  time  curve.  This  can  be 
explained  by  the  realisation  that  the  front  face  of 
the  detector  is  effectively  'black'  to  these  low  energy 
y-rays  and  X-rays.  Thus  total  absorption  would  occur 
in  a  region  of  poor  field  characteristics  of  the 
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detector.  Due  to  these  poor  fringing  fields  the  collec¬ 
tion  time  of  the  total  charge  formed  by  the  y-ray- 
detector  interaction  would  be  increased  and  hence  the 
centroids  of  the  time  distribution  curves  for  these 
y-ray  events  would  be  delayed.  However,  Compton  events 
at  the  same  energy  could  occur  throughout  the  volume  of 
the  detector  and  hence  charge  collection  times  would  not 
be  delayed  in  the  same  manner. 

To  investigate  the  phenomenon  at  higher  energies 
the  timing  properties  of  the  coaxial  Ge(Li)  detector 
were  studied  by  using  collimated  y-ray  sources.  The 
first  interesting  results  are  shown  in  fig.  35.  In 

22 

this  case,  a  timing  coincidence  was  set  up  using  a  Na 

source  and  collimating  the  radiation  to  the  Ge(Li) 

detector  by  a  6.3  cm  thick  lead  brick  with  a  hole  0.4  cm 

in  diameter  through  it.  The  collimated  beam  was  directed 

parallel  to  the  axis  of  the  detector  crystal  and  moved 

in  steps  of  0.5  cm  from  the  centre  to  the  outer  radius 

of  the  crystal.  Experiments  were  done  with  a  23  c.c. 

coaxial  Ge(Li)  detector  with  the  intrinsic  region  of 

the  crystal  having  dimensions  0.34  cm  inner  radius  and 

1.54  cm  outer  radius.  The  results  shown  are  for  the 

22 

system  gated  on  the  0.511  MeV  peak  from  the  Na  source. 

A  significant  improvement  in  the  time  resolution  is  seen 
at  the  centre  of  the  intrinsic  region  as  compared  to  the 
regions  near  the  inner  and  outer  electrodes  of  the 
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Fig.  35.  Time  distributions  obtained  as  a  collimated 
y-ray  beam  from  a  “^Na  source  was  moved  in 
steps  of  0.5  cm  from  the  centre  of  the 
detector  to  the  outer  radius,  the  timing 
system  being  gated  by  an  energy  window  on 
the  0.511  MeV  peak. 
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Fig.  36. 


Time  distributions  obtained  as  a  collimated 
y-ray  beam  from  a  ^Co  source  was  moved  in 
steps  of  0.5  cm  from  the  centre  of  the 
detector  to  the  outer  radius,  the  timing 
system  being  gated  by  an  energy  window  on 
the  1.332  MeV  y-rays  of  a  ^®Co  source. 
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detector.  As  shown  in  fig.  5,  the  centre  of  the  intrin¬ 
sic  region  corresponds  to  pulses  that  have  the  shortest 
rise-times . 

Similar  results  are  shown  in  fig.  36  for  a 

6  0 

collimated  y-ray  beam  from  a  Co  source  with  the 
electronics  gated  on  the  1.33  MeV  y-rays. 

In  this  case,  however,  due  to  the  higher 
energy  y-rays  the  lead  collimator  was  not  as  effective 
but  the  same  trends  are  seen  that  are  also  observed  in 
fig.  35. 

An  interesting  feature  of  these  curves  is  that 
even  at  the  centre  of  the  detector,  in  a  region  which 
should  be  insensitive,  radiation  is  detected.  The 
deterioration  of  the  early  side  of  the  time  resolution 
curves  using  the  ARC  timing  technique  can  be  explained 
by  events  which  occur  in  the  central  core  of  the  crystal 
and  also  in  the  intrinsic  region  near  the  inner  elec¬ 
trode.  Our  technique  of  pulse  shape  selection  which 
eliminates  these  pulses,  therefore  selects  those  events 
which  occur  in  the  centre  of  the  intrinsic  region  i.e. 
about  1.0  cm  from  the  centre  of  the  crystal. 

Fig.  37  shows  a  plot  of  the  centroids  of  the 
time  resolution  curves  of  figs.  35  and  36  versus  the 
collimated  beam  position.  The  points  illustrate  that 
y-ray-detector  interactions  occurring  towards  the  outer 
diameter  of  the  crystal  result  in  events  which  are  timed 
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A  plot  of  the  centroids  of  the  time 
distributions  of  figs.  35  and  36  vs. 
the  collimated  beam  position. 
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later.  This  is  consistent  with  the  theoretical  results 
shown  in  figs.  6  and  7  which  show  that  the  zero¬ 
crossing  time,  t  ,  occurs  later  as  the  interaction  posi- 
tion  moves  towards  the  outer  radius  of  the  detector. 

/■A  1  ^  *7 

Using  collimated  y-ray  beams  from  uCo,  /Cs, 

22 

Na  and  RdTh  sources  separately,  and  with  no  timing 
coincidence  in  this  case,  a  comparison  between  the 
number  of  Compton  and  full  energy  events  at  a  given 
pulse  height  was  considered  for  the  four  collimated 
beam  positions.  The  results  are  shown  in  fig.  38. 

They  were  obtained  by  setting  energy  windows  on  the 
full  energy  peaks  of  a  given  source  and  using  the  same 
windows  to  select  events  from  the  Compton  continua  of 
the  other  sources  where  applicable.  The  graphs  are 
plotted  so  that  at  the  collimated  beam  position  which 
shows  a  maximum,  is  normalised  to  unity.  From  these 
graphs  it  is  observed  that  there  is  a  greater  proportion 
of  Compton  to  full  energy  events  in  the  regions  near 
the  inner  and  outer  electrodes  of  the  detector.  In 
each  case  an  approximate  background  subtraction  was 
done  to  account  for  scattering  in  or  transmission  through 
the  lead  collimator,  since  perfect  collimation  is  effec¬ 
tively  impossible. 

These  results  are  not  entirely  unexpected.  In 
this  region  most  events  selected  from  a  Compton  contin¬ 
uum  almost  entirely  represent  single  interactions  in 
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Fig.  38.  The  normalised  number  of  Compton  and  full  energy 
events,  for  a  given  pulse  height,  plotted  as  a 
function  of  the  collimated  y-ray  beam  position. 
6°co ,  13 7Cs ,  22Na  and  RdTh  sources  were  used 
to  obtain  the  data. 
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the  detector.  However,  this  is  not  the  case  for  the 
full  energy  events  as  indicated  in  fig.  39  which  shows 
the  calculated  percentage  contribution  to  full  energy 
peaks  by  multiple  scattering  events.  Bengtson  and 
Moszyriski  (Be  72)  calculated  these  curves  using  a  Monte 
Carlo  technique  and  assuming  in  the  Compton  scattering 
problem  that  90°  was  the  only  possible  scattering  angle 
to  simplify  the  calculations. 

Thus  due  to  this  multiple  scattering  process 
there  is  a  greater  probability  of  obtaining  a  full 
energy  event  if  the  first  interaction  occurs  in  the 
centre  of  the  intrinsic  region  of  the  detector.  Other¬ 
wise  subsequent  interactions  would  usually  not  occur 
within  the  detector  volume  before  total  absorption  of 
the  energy  of  the  initial  y-ray  is  obtained.  This  would 
therefore  account  for  the  greater  percentage  of  Compton 
events  occurring  near  the  electrodes  of  the  detector. 

Using  these  results  combined  with  the  observa¬ 
tions  of  fig.  37  one  would  expect  to  find  the  centroids 
of  the  time  distribution  curves  from  Compton  events 
delayed  with  respect  to  those  obtained  from  full  energy 
events  at  the  same  pulse  height. 

Since  the  time  resolution  is  poor  near  the  elec¬ 
trodes  one  would  expect  to  observe  a  difference  in  the 
shape  spectra  across  the  face  of  the  detector  using 
collimated  y-ray  beams.  Figs.  40,  41  and  42  show  the 
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The  calculated  percentage  contribution 
to  full  energy  peaks  by  multiple 
interacting  events  (Be  72). 
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Fig.  40.  'Shape'  spectra  obtained  as  a  collimated  y-ray 

beam  from  a  RdTh  source  was  moved  from  the  centre 
of  the  detector  in  steps  of  0.5  cm  to  the  outer 
radius,  the  system  being  gated  by  an  energy  window 
on  the  0.352  MeV  y-rays  from  the  source. 
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Fig.  41. 


'Shape'  spectra  obtained  as  a  collimated 
y-ray  beam  from  a  l^Cs  SOUrce  was  moved 


of  0.5  cm  to  o 
being  gated  by 
0.661  MeV  y-rays 


tine  outer  raaius,  tne  sys 
an  energy  window  on  the 
from  a  ^7Cs  source. 
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Fig.  42.  'Shape'  spectra  obtained  as  a  collimated 
y-ray  beam  from  a  60Co  source  was  moved 
from  the  centre  of  the  detector  in  steps 
of  0.5  cm  to  the  outer  radius,  the  system 
being  gated  by  an  energy  window  on  the 
1.332  MeV  y^-rays  of  a  °^Co  source. 
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shape  spectra  obtained  from  the  four  collimated  beam 

positions  at  various  energies.  The  figures  show 

results  for  the  0.352  MeV  y-ray  from  the  RdTh  source, 

13  7 

the  0.661  MeV  y-ray  from  the  Cs  source  and  the 

6  0 

1.332  MeV  y-ray  from  the  Co  source.  Similar  results 
were  obtained  for  other  full  energy  peaks.  The  arrows 
on  the  diagrams  illustrate  the  region  used  in  our 
pulse  shape  selection  technique  and  correspond  appro¬ 
ximately  to  bins  4  to  7  in  figs.  17  and  19,  that  is 
those  bins  which  give  the  most  improved  time  resolution. 

These  figures  show  the  true  significance  of  the 
pulse  shape  selection  technique.  The  technique  effec¬ 
tively  selects  those  pulses  which  occur  due  to  inter¬ 
actions  in  the  centre  of  the  intrinsic  region  while 
eliminating  a  great  proportion  of  those  events  which 
occur  near  the  electrodes. 

The  left  hand  hump,  which  could  not  be  theoreti¬ 
cally  obtained  as  illustrated  in  figs.  10  and  11,  is 
seen  to  come  mostly  from  events  which  occur  in  the 
central  core,  or  in  what  should  be  the  insensitive 
region  of  the  detector.  Obviously  events  which  occur 
in  this  region  produce  electron-hole  pairs  which  are 
still  collected  by  the  electrodes  but  with  very  long 
collection  times  which  give  rise  to  pulses  with  poor 
time  resolution  characteristics.  Thus  by  drilling  out 


this  central  core,  coaxial  Ge(Li)  detector  should  have 
better  time  resolution. 

Similar  shape  spectra  were  obtained  for  energy 
windows  set  on  the  Compton  continua  of  the  sources, 
however,  no  significant  visual  difference  could  be  seen 
between  these  spectra  and  those  obtained  from  full 
energy  events.  The  slight  difference,  however,  shows 
up  when  one  normalises  the  shape  spectra  for  a  given 
energy  and  collimated  beam  position  and  then  calculates 
the  total  number  of  events  in  that  region  of  the  shape 
spectra  which  gives  good  time  resolution.  Fig.  43  shows 
the  results  of  these  calculations  plotted  as  the  full 
energy  to  Compton  ratio,  for  a  given  pulse  height,  versus 
the  collimated  beam  position.  The  same  trends  are 
shown  again  as  seen  in  fig.  38  that  there  is  a  greater 
percentage  of  Compton  to  full  energy  events  in  those 
regions  near  the  electrodes  of  the  detector. 

These  slight  differences  are  further  accentuated 
when  one  considers  the  volume  contributions  at  a  given 
radius  of  the  detector.  The  full  face  of  the  detector 
was  irradiated  with  the  sources  separately  and  fig.  44 
shows  the  total  shape  spectra  obtained  from  energy 
windows  set  on  the  full  energy  peaks  and  Compton  con¬ 
tinuum  for  the  same  pulse  height.  Curves  are  shown  for 
the  energies  0.661  MeV,  1.172  MeV  and  1.332  MeV,  and  in 
each  section  of  fig.  44  the  two  curves  are  normalised 
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Relative  peak  to  Compton  ratios,  for  given 
pulse  heights,  plotted  as  a  function  of  the 
collimated  beam  position  with  only  events 
in  the  'good  shape'  region  being  selected. 
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'Shape'  spectra  obtained  by  irradiating  the  full 
face  of  the  detector  and  gating  the  system  on  full 
energy  events  from  one  source  and  the  Compton  con¬ 
tinuum  of  another,  using  the  same  energy  window. 


Fig.  44. 
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to  the  same  number  of  counts.  It  is  seen  that  there  is 
a  marked  difference  in  the  shape  spectra  obtained  from 
Compton  events  compared  to  those  obtained  from  full 
energy  events  especially  at  higher  energies.  It  is 
expected  that  the  shape  spectra  become  similar  at  low 
energies  since  in  this  case  y-rays  have  a  greater  pro¬ 
bability  of  depositing  all  their  energy  in  the  detector 
in  a  single  interaction. 

Similar  differences  between  Compton  and  full 
energy  events  have  been  seen  before  (St  67,  Mo  70, 

Mo  72)  using  different  techniques.  The  greater  pro¬ 
portion  of  Compton  events  in  the  left  side  of  the 
shape  spectra  i.e.  the  region  associated  with  poor 
time  resolution,  can  be  explained  by  noting  that  the 
regions  near  the  electrodes  of  a  Ge(Li)  detector  are 
often  poorly  compensated,  the  field  therefore  being 
weaker  and  enhancing  carrier  recombination  and  charge 
trapping  effects.  Thus  y-rays  interacting  in  these 
regions  often  give  rise  to  defective  slow-rising  pulses 
with  poor  time  resolution  characteristics.  Strauss  and 
Larsen  (St  67)  have  stated  that  much  of  the  continuum 
at  1.33  MeV  was  due  to  these  defective  pulses. 

To  finalise  the  study  of  the  timing  properties 

of  coaxial  Ge(Li)  detectors,  a  timing  coincidence  was 

22 

set  up  using  a  collimated  beam  from  a  Na  source, 
directed  at  0.5  cm  intervals  along  the  side  of  the 
detector,  perpendicular  to  the  axis  of  the  crystal. 
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Fig.  45  shows  the  time  resolution  curves  obtained  with 
the  electronics  gated  on  the  0.511  MeV  peak  at  1.0  cm 
intervals  from  the  front  face  to  the  back  face  of  the 
crystal.  No  shape  selection  was  involved  in  these 
results.  It  is  observed  that  the  time  resolution  is 
poor  at  both  ends  of  the  detector  where  the  electric 
field,  a  fringing  field,  is  poor. 

A  plot  of  the  centroids  of  these  time  distribu¬ 
tion  curves,  and  those  at  other  energies,  against  the 
position  of  the  collimated  beam  are  shown  in  fig.  46. 
This  result  is  quite  unexpected,  the  timing  coming 
earlier  at  the  back  face  of  the  detector.  It  could 
be  that  the  geometry  of  the  particular  detector  used 
resulted  in  a  variation  in  the  electric  field  along 
the  length  of  the  detector  crystal.  However,  as  no 
measurements  have  been  done  to  see  whether  a  similar 
result  occurs  with  other  crystals  no  conclusion  is 
possible . 
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Fig.  45.  Time  distributions  obtained  as  a  collimated 
y-ray  beam  from  a  22Wa  source  was  moved  in 
steps  of  1.0  cm  from  the  front  face  to  the 
back  face  of  the  detector  crystal,  the  timing 
system  being  gated  on  the  0.511  MeV  peak. 
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Fig.  46.  A  plot  of  the  centroids  of  the  time 
distributions  vs.  collimated  beam 
positions  along  side  of  detector 
crystal,  the  system  being  gated  by 
energy  windows  at  200,  300  and  511 
keV. 


CHAPTER  6 


ELECTROMAGNETIC  TRANSITION  RATES  AND  PROBABILITIES 


By  Heisenberg's  uncertain  principle,  the  meanlife 
of  an  excited  nuclear  state,  is  defined  in  terms  of  the 
energy  width  T  for  a  ydecay  as 

Tt  =  fi  =  0.6582  x  10  eV  sec  .  6-1 

The  radiative  transition  probability  is  the  reciprocal 
of  the  meanlife 


which  is  the  sum  of  the  partial  transition  probabilities 


T  =  E  T (oL) 
oL 


6-2 


where  L  is  the  multipolarity  of  the  decay  and  a  repre¬ 
sents  the  mode,  either  electric  or  magnetic. 

The  partial  transition  probability  of  a  gamma 
transition,  angular  momentum  LM  and  energy  E,  from  the 
initial  level  spin  JL  to  the  final  level  with  spin 
is  given  by  (Mo  65,  Sk  66) 


T(aLM)  = 


8  7T  (L  +  1) 

L [ (2L+1) ! ! ] 


-  i  (^-) 

2  ft  'tic' 


2L+1  p 

<f  MTm  i> 
i  i  lm 1  1 
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where  HI.,  is  the  corresponding  multipole  operator. 
LM 
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Usually  one  is  not  interested  in  the  orientation  of 
either  the  initial  or  final  state,  so  one  sums  over 
magnetic  substates  M^  and  averages  over  JVb  .  Defining 
the  so-called  reduced  matrix  element 


B  (  gL  ,  Ji^Jf)  =  (2Ji+  1) 


-1 


M.M, 
l  f 


a  ,  ,2 

<f  MJt-  i> 

LM 1  1 


The  partial  transition  probability  becomes 
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T(aL,  -  .8".<L  +  11  „  l  (",2L+1  B (gL) 


L[  (2L+  1)  !  !  ] 


2  h  tie' 
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The  electric  and  magnetic  operators  for  a  single 
particle  are  given  by 


M 


LM 


t  *  u  q  ^  T 

er  Ylm  -  i  — - - (L+l)  axr*[grad(r  Y^)] 
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and 


MLM  =  (lTT  g£  !:  +  gs  !)-[grad(rL  ylm)]*  , 
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respectively.  The  first  term  in  each  denotes  the  con¬ 
tribution  due  to  the  orbital  motion,  while  the  second 
term  arises  from  the  intrinsic  particle  spin.  is 

the  nuclear  magneton  and  the  g's  are  the  gyromagnetic 
ratios . 

From  these  equations  it  follows  that  any  2L  pole 
multipole  radiation  field  carries  angular  momentum  L 
and  z-component  M  and  its  parity  being  (-1)L  or  (-1)L  ^ 
for  EL  or  ML  radiation,  respectively.  Hence  from  the 


Ill 


conservation  of  angular  momentum  2^  pole  radiation  can 
occur  only  between  states  whose  angular  momenta  obey 
the  selection  rules 


J  Ji -  Jf |  <  L  <  Ji +  j 

and 

M  =  Mi  -  Mf 


From  the  conservation  of  parity,  one  has 


ATT  =  7Ti/7Tf  =  (“)L  for 

EL 

radiation 

=  (-1)L  for 

ML 

radiation 

where  tt^  and  are  the  parities  of  the  nuclear  states 

involved  in  the  transition. 

In  order  to  set  a  standard  for  comparison  of  all 

experimental  transitions,  regardless  of  any  model  of 

the  states  involved,  the  so-called  Weisskopf  and 

Moszkowski  units  were  defined.  Since  one  does  not 

distinguish  neutron  from  proton  transitions,  for  EL 

transitions  involving  a  single  proton,  the  main  con- 

E 

tribution  to  is  provided  by  the  orbital  term. 

For  a  central  potential  the  wavef unctions  of  the  initial 
and  final  states  can  be  written  as  a  product  of  a  radial 
function  which  depends  on  the  details  of  the  potential 
and  a  two-component  function  which  contains  angle  and 
spin  variable  but  is  independent  of  the  details  of  the 
potential.  Thus  the  electric  multipole  matrix  element 
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can  be  written  as 


<f  M 


E 

LM 


i> 


e 


rL  r2  dr 

0 


Mf  *  M. 

0JlfJf  YLM  e£^Jf 


where  dft  refers  to  integration  over  all  angles  and 
sum  over  spins.  To  evaluate  the  reduced  matrix 
element  a  'statistical  factor'  S  is  defined  as  follows: 


-1 

S(J±L  Jf)  =  47T(2Ji+l) 


E 

M.M* 
1  f 


Mr*  M.  2 

°£fJf  YLM  0£iJf  ^1 
4  7T 
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which  simplifies 


2 

S(J.LJf)  =  (2Jf  +  1)  (J-J^  -  ^|L0) 
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Hence  the  reduced  matrix  element  becomes 


B  (EL) 


2  2L 
e  R 

4  77 


(  (Rf  (|)L  ^  r2  dr)2  S(J.LJf) 
0 
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In  order  to  obtain  simple  standard  units  the  radial 
wavef unctions  (ft  for  both  initial  and  final  states  are 
assumed  constant  throughout  the  interior  of  the  nucleus 
(for  r  <  a)  and  vanish  outside  (for  r  >  a) .  On  the 
basis  of  this  constant  density  model,  the  values  of 
the  radial  integrals  are  given  as  follows: 

oo 

{  dif  (^-)L  <ftL  r2  dr  =  3/  (L  +  3)  . 

0 
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. 
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A  further  assumption  is  made  that  the  radiation  of 
multipole  L  occurs  between  the  states  =  L  +  \  and 
Jf  =  k.  Hence 

S(L+h,  L,k)  =  1  •  6-12 


Thus  substituting  eqs.  6-11  and  6-12  in  eq.  6-10  one 
obtains  the  Weisskopf  or  Moszkowski  unit  for  an  EL- 
transition 


2  2L  _  2 

VEL)  -  t!-  <lT3> 


4 it  (L+3 


2 

)  1.2 


1  ,  3  _2L  2L/3  2L  2 


A 


[fm  e  ] 
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2L  2 


in  units  of  fm  e  where  the  following  relation  has  been 
used  for  the  nuclear  radius 


R  =  r  A 
o 


1/3 


and  where  r  =1.2  fm  and  A  is  the  atomic  number  of  the 

o 


nucleus . 

Similarly  a  Moszkowski  unit  for  magnetic  tran¬ 
sitions  can  be  calculated: 


bm(ml)  =  !  V-ETi  (lT2>  1-2“"  '  A 

2L-2  2 


2L-2  _  (2L-2 ) / 3 


[  f  m 


y  ] 
o 
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The  Weisskopt  unit  for  magnetic  transition  is  a  numerical 
approximation  of  this  Moszkowski  unit 
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B  (ML)  = 
w 


10_ 

7T 


(LT3}  X‘2 


2L-2  _  (2L-2)/3  r -  2L-2  2, 

A  [fm  yQ] . 
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It  must  be  emphasized  that  the  above  theoretical  es¬ 
timates  are  crude  and  are  not  expected  to  represent 
accurately  transitions  in  actual  nuclei.  However, they 
provide  a  starting  point  for  interpretation  of  known 
results . 

Several  convenient  units  have  been  tabulated 
(Sk  66)  for  use  in  the  above  formulae: 

"he  =  1.973  x  10  ^  MeV  cm 

e2  =  1.440  MeV  fm 

y  =  0.01589  MeV  fm2  . 
o 

Another  parameter  of  significant  importance  in 
the  calculation  of  partial  transition  probabilities  is 
the  E2/M1  multiple  mixing  ratio  6  which  may  be  defined 
by  the  following  relations: 

r2  _  T  (E2) 

6  T (Ml ) 

T  (Ml)  =  — y  T 

1  + 

*2 

T  (E2  )  =  — - - y  T 

1  +  6 

and 


T  =  T (Ml)  +  T  (E2 )  . 


If  the  nuclear  state  should  decay  by  other  competing 
processes  such  as  internal  conversion,  the  radiative 


decay 


lifetime  t 

T  =  T  , 

y  obs 


must  be  corrected 

v 

(1  +  aT) 


by 


where 

total 


t  ,  is  the  observed  lifetime  and  a 
obs 

internal  conversion  coefficient. 


is 


the 


CHAPTER  7 


EXPERIMENTAL  RESULTS 

7 . 1  Establishment  of  Suitable  ’Prompt1  Reference  y-rays 

One  of  the  most  important  problems  faced  in  this 
research  was  the  definition  of  the  prompt  time  centroid 
curve.  A  sufficient  number  of  centroid  points  corres¬ 
ponding  to  y-ray  transitions  from  states  whose  lifetimes 
are  short  compared  to  the  time  resolution  of  the  system, 
less  than  say  5  ps ,  are  required.  The  usefulness  of 
the  technique  and  the  validity  of  the  results  hinge  on 
finding  a  solution  to  this  problem. 

The  shortest  lifetimes  that  can  be  measured  by 
this  technique  are  of  the  order  of  10  ps ,  and  lifetimes 
of  this  magnitude  or  greater  are  most  frequently 
encountered  in  low  lying  energy  levels  of  nuclei, 
typically  at  energies  of  less  than  300  keV.  In  this 
energy  region  the  walk  of  the  prompt  time  centroid 
curve  is  pronounced  and  becomes  worse  as  the  energy 
decreases.  Hence,  it  is  absolutely  necessary  to  have 
sufficient  reference  'prompt*  y-rays  in  this  region 
where,  unfortunately,  there  are  very  few  'prompt'  y- 
ray  transitions.  Unless  an  isotope  is  found  which  has 
several  'prompt'  y-ray  transitions  one  has  to  resort  to 
a  'smorgasbord'  target  consisting  of  several  different 
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isotopes,  each  having  one  or  two  'prompt*  Y“raY  transi¬ 
tions,  in  order  to  adequately  define  the  prompt  time 
centroid  curve.  This  would  lead  to  problems  in  target 
fabrication  and  also  the  possibility  of  many  unwanted 
y-ray  transitions  and  nuclear  reactions. 

In  some  cases  states  which  have  measured  life¬ 
times  less  than  60  ps  can  be  used  to  help  define  the 
prompt  time  centroid  curve  by  adjusting  their  centroids 

amounts  equivalent  to  their  measured  lifetimes.  A 

181 

literature  search  yielded  two  suitable  isotopes,  Ta 
and  ^~*Ho,  each  with  many  y-ray  transition  from  states 
with  lifetimes  of  60  ps  or  less.  Another  advantage  of 
these  two  isotopes  is  that  they  occur  as  100%  of  the 
natural  elemental  abundance. 

181 

Fig.  47  shows  the  decay  schemes  of  Ta  (Ro  70, 
16  6 

NDS  73)  and  Ho  (NDS  74)  for  the  low  energy  states. 

For  proton  beams  of  less  than  5.5  MeV  the  states  which 

are  strongly  populated  are  the  136  and  301  keV  levels 

in  "^^Ta  and  the  95  and  210  keV  levels  in  ^~*Ho.  The 

y-ray  transitions  of  interest,  the  136,  165  and  301  keV 

181 

transitions  in  Ta  and  the  95,  115,  and  210  keV  tran- 

sitions  in  °Ho,  are  shown  in  bold  ink  in  fig.  47. 

181 

The  lifetime  of  the  136  keV  level  in  Ta  has 
been  measured  using  several  techniques.  Values  of 
60±5  ps  (Bl  60)  and  61.4±5.4  ps  (Ar  70)  were  obtained 
using  a  microwave  beam  pulsing  system.  Steiner  et  al. 
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Fig.  47. 
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Experimental  level  schemes  ;  or 

1  O  rr>i  v  —  ./  franci  '-in^c;  q  f 


’Ho.  The  7- ray  transition£ 
are  shown  Ly  bolder  lines 


Ta  and 
interest 


(St  69)  obtained  a  value  of  55±3  ps  using  the  Mossbauer 
Effect  and  the  weighted  mean  of  several  Coulomb  Excita¬ 
tion  measurements  yields  a  reduced  transition  probabi¬ 
lity  B(E2)t  =  (2.0±0.4)  x 104  e2fm4  from  which  a  lifetime 
of  55±11  ps  is  calculated  (NDS  73).  The  weighted  mean 
of  these  results,  57.2±2.3  ps ,  has  been  used  in  this 
work . 

The  reported  lifetime  of  the  301  keV  level  in 

181 

Ta  is  21.6±1.7  ps  (NDS  73)  which  has  been  calculated 

from  an  average  value  of  B(E2)f  =  (0 . 50±0 . 04 ) xio4  e2fm4. 

The  lifetime  of  the  95  keV  level  in  ^~*Ho  has 

also  been  measured  using  several  techniques.  Armbruster 

et  al. (Ar  70)  obtained  a  value  of  40.6±3.4  ps  with  a 

microwave  beam  pulsing  system.  Blunck  et  al.(NDS  74) 

using  the  Mossbauer  Effect  obtained  a  value  of  31.7±0.6 

which  is  in  agreement  with  the  value  of  32±1  ps  (NDS  74) 

calculated  from  the  B(E2)1  result  (01  60).  A  weighted 

mean  of  32.0±0.5  ps  was  adopted  for  this  work. 

165 

The  lifetime  for  the  210  keV  state  in  Ho  is 
18.8±1.2  ps  (NDS  74)  calculated  from  the  measured  B(E2)f 
value  of  (01  60)  of  (0.6310.04)  x 104  e2fm4. 

To  define  the  prompt  time  centroid  curve  at  even 
lower  energies  the  X-rays  of  both  4^4Ta  and  4^~*Ho  were 

used.  These  include  the  K  X-rays  at  56.28  and  57.54 

a 

181 

keV  and  the  K0  X-ray  of  65.20  keV  in  Ta  and  the  K 

3  a 

X-rays  of  46.70  and  47.55  keV  in 


Ho.  X-rays  of  these 
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energies  have  very  short  lifetimes  which  are  illustrated 
by  the  measurement  of  the  width  of  the  X-ray  transi¬ 
tion  in  19"^Au  at  68.81  keV  to  be  54  eV  (He  54).  This 
is  equivalent  to  a  lifetime  of  0.012  fs. 

These  X-rays  together  with  the  y-ray  transitions 
discussed  above  define  the  prompt  time  centroid  curve 
very  well  below  165  keV  as  illustrated  in  fig.  48. 

These  results  were  obtained  using  a  5.0  MeV  pulsed 

2 

proton  beam  on  two  targets,  one  being  a  1  mg/cm  en- 
64 

riched  Ni  foil  and  the  other  composed  of  two  sections 

2  163  2 

consisting  of  a  110  mg/cm  DHo  foil  and  a  156  mg/ cm 

l^Ta  foil.  The  centroids  at  95,  115,  136  ,  165  ,  210  and 

310  keV  were  adjusted  to  take  into  account  the  measured 

lifetimes  of  the  states  from  which  these  transitions  occur. 

To  define  the  prompt  time  centroid  curve  at 

energies  above  150  keV  more  accurately, a  detailed 

64 

investigation  of  the  lifetimes  of  the  states  in  Cu 

64 

was  performed.  Cu  was  chosen  because  its  y-ray 

decay  scheme  (Gr  73)  has  many  transitions  between  159 

and  600  keV,  an  important  energy  region  for  lifetime 

investigations.  Recently  a  detailed  investigation  of 

the  properties  of  the  excited  states  and  y-ray  tran- 
64 

sitions  in  Cu  was  done  in  this  laboratory  for  a 

Ph.D.  dissertation  by  Green  (Gr  73).  Fig.  49  shows 

64 

the  level  scheme  for  Cu  below  927  keV.  The  states 
in  ^Cu  are  excited  via  the  ^Ni  (p ,n)  ^Cu  reaction 
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Fig.  48. 


Centroids  of  the  time  distributions  for 
full  energy  events  using  the  ^Ni  foil 
and  the  foil  combination. 


122 


|  (3)  -59  32  9 


3 

0 


3 

(2,3) 

I 


1+) 


4(2) 


3 
l  + 


2(*L± 


l  + 


49  40 


5Z355. 


4  13  6  7? 


6  3  5  2732 


6  4  8  82 


98  2 


4  96 


927 

895 

878 


746 

739 


663 


609 


574 


362 

344 

278 


A _ 1  w  w  Q 


159 


64 


Cu 


(Gr  73) 


Fig.  49. 


.experimental  level  scheme  of 


o4Cu  below  927  keV. 
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which  has  a  Q  value  of  -2.458  MeV.  An  illustration  of 
the  experimental  results  for  the  159,  203,  212  and 
278  keV  transitions  are  shown  in  fig.  48,  the  centroids 
being  indicated  by  crosses. 

Several  other  second  targets  were  used  in  order 

to  provide  'prompt'  y-rays  to  determine  the  lifetimes 

64  3  9  7 

of  the  levels  in  Cu.  The  278  keV  level  in  Au  which 

decays  to  the  ground  state  has  a  lifetime  of  21.6±1.2  ps 

(NDS  72)  which  is  calculated  from  the  B(E2)t  and  the 

mixing  ratio  measurement  by  McGowan  et  al . (Me  71). 

197 

Thus  by  using  a  Au  foil  target  the  lifetime  of  the 

64 

level  at  278  keV  m  Cu  could  be  measured. 

Two  other  targets  used  were  NaCl  evaporated  on 

181 

Ta  and  natural  copper  foil.  The  first  excited  state 
23 

at  439  keV  in  Na  has  a  value  cited  by  Endt  and  van  der 

Leun  (En  73)  of  1600180  fs  and  the  first  three  excited 

states  of  ^Cu  at  669,  962  and  1327  keV,  and  ^Cu  at 

771,  1115  and  1482  keV  all  have  lifetimes  less  than 

1  ps  (Wo  70).  A  final  comparison  target  used  was  a 
2  4  7 

1  mg/cm  enriched  Ti  foil.  The  (p,n)  reaction  excites 

4  7 

the  states  in  V.  The  level  at  660  keV  has  been  measured 

+  2  8 

to  have  a  lifetime  of  4_2*  ps  by  the  recoil  distance 
technique  (B1  73).  This  state  has  three  decay  modes; 
a  48%  branch  to  the  ground  state,  a  15%  branch  to  the 
first  excited  state  via  a  572  keV  y-ray  and  a  37%  branch 
to  the  third  excited  state  via  a  400  keV  y-ray  transition. 


An  example  of  an  experimental  run  using  the 
64  .  47 

Ni  and  Ti  foils  is  shown  in  fig.  50.  Eighteen 

experimental  runs  were  performed  to  obtain  the 

64 

lifetimes  in  Cu  and  the  results  are  shown  in 
table  1.  The  formulation  used  to  calculate  the 
weighted  average  is  that  employed  in  Nuclear  Data 
Sheets:  If  x^±Ax^,  •  *  *Xn±Axn  are  n  independent 

measurements  of  a  given  quantity,  the  weighted  average 
of  these  measurements  is  x±Ax  where 


and 


x  =  W 


I  xi/(Axi)2 


and  Ax  is  the  larger  of  (W) ^  and  [w£{(Ax^)  (x-x^)  }/n-l] 

In  those  cases  where  no  centroid  shift  is  observable 
the  result  is  quoted  as  an  upper  limit  using  the  error 
Ax  as  this  limit. 

64 

The  results  shown  m  table  1  indicate  that  Cu 


is  an  excellent  source  of  'prompt'  reference  Y-rays. 
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Table  1 


Lifetime  results  for 


64Cu 


Level  Transition  to  Branching  Ratio  Lifetime 

Ei(keV),  1 1  Ef(keV),  l£  (%)  (ps) 


159 

,  2  + 

0 

/ 

i+ 

100 

29 . 8 ±5 . 4 

278 

,  2{+) 

0 

/ 

i+ 

100 

6  +  7 

6-6 

344 

,  1+ 

0 

r 

i+ 

96 

<5 

362 

,  3 

159 

9 

2  + 

98 

<6 

574 

,  4(2) 

362 

9 

3 

100 

9  +  15 
-9 

609 

2  + 
f  ^ 

0 

9 

1+ 

82 

<13 

633 

,  i 

'  0 

9 

1  + 

32 

159 

9 

2 

27 

<12 

278 

9 

2  (  +  ) 

35 

739 

,  (2,3) 

159 

9 

2+ 

77 

<16 

746 

,  3 

278 

9 

2  (  +  ) 

100 

<18 

878 

,  o 

344 

9 

1+ 

42 

<22 

895 

,  3 

278 

9 

2  (  +  ) 

49 

<29 

927 

,  1(3) 

278 

9 

2  (  +  ) 

59 

<15 
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7 . 2  A  Lifetime  Study  of  the  Isotopes  of  Titanium  and 

and  Vanadium 

The  isotopes  of  Titanium  and  Vanadium  have  many 
low  energy  excited  states  whose  lifetimes  have  not  been 
measured  or  where  previous  measurements  show  disagree¬ 
ments  in  results.  It  was  felt  therefore  that  an 
investigation  of  these  levels  would  be  useful  in  obtain¬ 
ing  several  lifetimes,  illustrating  the  potential  of  the 
technique,  and  would  furnish  more  information  on  those 
discrepancies  which  occur  in  previous  reported  measure¬ 
ments  . 

1  3 

A  complete  set  of  1^-"  *  rolled  isotopically 

enriched  Titanium  foils  were  obtained  from  Micromatter, 

64 

Seattle.  In  addition  a  Ni  foil  was  also  obtained. 

The  relevant  information  on  the  isotopic  analysis  and 
thicknesses  is  given  in  table  2. 

45  . 

7.2a  °Ti 

45  . 

The  decay  scheme  for  Ti  as  given  by  Jett  et  al 
(Je  68)  and  Iyengar  and  Robertson  (Iy  71)  is  shown  in 
fig.  51.  Lynch  et  al  (Ly  70)  have  measured  the  life¬ 
times  of  the  first  and  second  excited  states  to  be 
4 . 5±0 . 5  ys  and  17.211.0  ns,  respectively,  using  a  post 
acceleration  pulsed  proton  beam.  The  second  excited 
state  has  also  been  measured,  using  the  pulsed  proton 
beam  technique,  by  Blasi  et  al  .  (B1  69)  to  be  16.612.2  ns. 
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Table  2 


Target 

Percent  Isotopic  Composition 

Thickness 

A=4  6 

47 

48 

49 

50 

mg/ cm2 

46Ti 

83.8 

5.0 

9.76 

0.73 

0 . 73 

1.01 

4  7Ti 

1.87 

80.1 

15.8 

1.11 

1.1 

1.001 

4  8Ti 

0.17 

0.2 

99.36 

0.17 

0.11 

0.981 

49  . 

Ti 

1.33 

1.25 

14.09 

81.61 

1.73 

1.01 

50m . 

Ti 

1.7 

1.3 

12.4 

1.5 

83.2 

1.01 

A=5  8 

60 

61 

62 

64 

64AT. 

Ni 

0.92 

0 . 73 

0.05 

0.38 

97.92 

0.99 
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Fig.  51.  Experimental  level  scheme  and  previous 

lifetime  results  for  the  states  in  ^5Ti. 
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The  lifetimes  of  the  higher  energy  levels  have 
been  measured  by  Iyengar  and  Robertson  (Iy  71)  and 
Zuk  et  al . (Zu  72),  both  using  the  Doppler  Shift  Atten¬ 
uation  Method.  Their  results  are  shown  in  fig.  51. 

The  levels  of  interest  in  this  work  are  those  which 
have  lower  limits  i.e.  the  levels  at  330,  744  and 
1227  keV.  However,  these  levels  have  been  investigated 
by  Blasi  et  al  .  (B1  71)  using  the  recoil  distance  tech¬ 
nique  through  the  ~^p  (^o  ,pny)  ^~*Ti  reaction  at  an  "^0 
beam  energy  of  36  MeV.  The  following  lifetimes  were 
found;  330  keV  level,  1.72±0.10  ns;  744  keV  level, 
14.5±2.5  ps ;  1227  keV  level,  4.0±0.8  ps . 

The  two  targets  used  in  the  investigations  were 

45  2  2  64 

a  Sc  foil,  48  mg/ cm  ,  and  the  1  mg/cm  enriched  Ni 

foil.  Pulsed  proton  beams  of  energies  4.2  to  4.5  MeV 

45  64 

were  used  to  excite  the  states  in  Ti  and  Cu  via 

45  45 

the  (p,n)  reaction,  the  Q  value  for  the  Sc(p,n)  Ti 
reaction  being  -2.83  MeV.  An  illustration  of  the 
results  is  shown  in  fig.  52.  The  weighted  mean  of 
six  experimental  runs  gave  a  lifetime  for  the  330  keV 
level  of  1.585±0.020  ns.  This  result  is  slightly 
smaller  than  that  quoted  by  Blasi  et  al  (B1  71)  of 
1.7210.10  ns.  The  293  keV  y-ray  transition  from  the 
330  keV  level  has  been  found  to  be  isotopic  (Iy  71, 

Zu  72)  and  therefore  no  mixing  ratio  has  been  found. 
Assuming  the  transition  to  be  pure  El,  our  result 
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which  yields  an  experimental  El  transition  width 
T (El)  =  4.15  x 10  eV  which  using  the  Weisskopf  esti- 

_5 

mate  corresponds  to  an  El  strength  of  1.9  x  10  W.u. 

For  a  pure  M2  mode  of  decay,  the  M2  transition  strength 
is  1033  W.u.  M2  enhancements  of  this  magnitude  are 
not  found,  however,  the  El  strength  is  not  an  unreason¬ 
able  result  and  thus  the  transition  is  expected  to  be 
mainly  El. 

No  observable  shift  was  seen  for  the  414  keV 

Y-ray  transition  from  the  744  keV  level  and  an  upper 

limit  of  <16  ps  was  calculated  for  this  level.  The 

parity  of  the  5/2  spin,  744  keV  level  is  not  known, 

however,  Blasi  et  al . (B1  71)  have  suggested  a  positive 

parity  since  this  state  could  be  the  5/2+  member  of  a 

possible  rotational  band.  The  414  keV  transition  has 

been  found  to  be  isotropic  (Iy  71,  Zu  72)  and  thus  no 

mixing  ratio  has  been  established.  Assuming  the  lowest 

mode  of  decay,  a  positive  parity  for  the  744  keV  level 

would  result  in  a  Ml  strength  of  0.03  W.u.  and  a  nega- 

-4 

tive  parity  would  lead  to  an  El  strength  of  7.5  x  10 
W.u.  Both  these  strengths  are  not  uncommon  and  hence 
the  parity  of  the  state  cannot  be  established  from  life¬ 
time  considerations. 

47 

7.2b  Ti 

4  7 

Fig.  53  shows  the  decay  scheme  of  Ti  (We  72). 
Also  indicated  are  the  results  of  lifetime  measurements 
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using  the  Doppler  Shift  Attenuation  Method  obtained 

by  Weaver  et  al . (We  72).  Of  the  higher  energy  levels, 

the  states  at  1.794,  1.825,  2.365  and  2.416  MeV  are 

possible  candidates  for  lifetime  measurements  by  our 

technique.  However,  the  two  levels  at  2.365  and  2.416 

MeV  did  not  provide  sufficient  yield  using  the  (p,p') 

reaction  reaction  to  make  measurements  possible. 

For  the  160  keV  -*  g .  s .  transition  Weaver  et  ah 

(We  72)  obtained  a  mixing  ratio  of  6  =  0.02  7±0.00  7. 

The  B(E2)f  value  for  the  transition  has  been  determined 

using  Coulomb  Excitation  by  Temmer  and  Heydenburg 

(Te  56)  and  Ritter  et  al  «  (Ri  62)  who  obtained  respec- 

2  4  2  4 

tive  values  of  396  e  fm  and  280  e  fm  with  ±15%  errors. 

These  values  correspond  to  a  weighted  mean  value  of 

B(E2)l  of  239±41  e^fm^.  Using  this  result  and  their 

mixing  ratio.  Weaver  et  al 0 calculated  a  lifetime  of 

25^12  Ps  f°r  t*ie  state.  However,  Holland  and 

4 

Lynch  (Ho  61)  using  a  post  acceleration  pulsed  He 
beam  and  a  Nal(Tl)  scintillator  detector  obtained  a 
lifetime  of  320±100  ps  for  the  160  keV  state.  A  re¬ 
measurement  of  the  lifetime  of  this  level  is  therefore 
indicated. 

Measurements  were  done  on  the  160  keV  level 
4  7  64 

using  the  Ti  foil  and  the  Ni  foil  attached  to  the 
2  181 

156  mg/cm  Ta  foil.  The  results  of  one  of  the 

experimental  runs  are  shown  in  fig.  54.  The  two  points 
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at  172  and  260  keV  correspond  to  the  decay  of  the  third 

4  7 

excited  state  in  V  which  will  be  discussed  in  a  later 
section.  The  weighted  mean  of  eight  runs  gave  a  life¬ 
time  of  301±8  ps  for  the  160  keV  state  which  is  in 
agreement  with  the  value  obtained  by  Holland  and  Lynch. 

Using  the  weighted  mean  value  for  the  B(E2)1 
and  then  our  result  for  the  lifetime, a  mixing  ratio  of 
6  =  0.095±0.003  is  obtained.  From  these  results  the 
160  keV-^-g.s.  transition  has  an  E2  strength  of  24  W.u. 
and  a  Ml  strength  of  0.03  W.u. 

The  levels  at  1.795  and  1.825  MeV  were  studied 
4  7 

using  the  Ti  foil.  Two  second  targets  were  used,  a 

^Ti  foil  and  a  ~^V  target  prepared  by  mixing  powdered 

V2O5  with  benzene  and  polyurethane  and  smearing  the 

resulting  glue  thinly  on  a  piece  of  gold  leaf.  This 

target  is  commonly  referred  to  as  a  gloop  target. 

Kavaloski  and  Kossler  (Ka  69)  have  studied  the  levels 
4  8 

in  Ti  up  to  3.74  MeV  excitation  energy  using  the 

Doppler  Shift  Attenuation  Method.  Their  results  and 

level  scheme  are  shown  in  fig.  55.  Szoghy  et  al  (Sz  73) 

51 

used  the  same  technique  to  study  the  levels  in  Cr 

below  4.1  MeV  excitation.  Their  results  and  level 

51  51 

scheme  are  also  shown  in  fig.  55.  The  V(p,n)  Cr 

51 

reaction  Q  value  is  -1.532  MeV,  so  that  the  V  target 
4  8 

plus  the  Ti  target  provide  many  'prompt'  reference 

4  7  . 

y-rays  for  the  study  of  the  two  levels  m  Ti .  As 
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shown  in  fig.  53  supplementary  'prompt'  y-rays  are 

4  7 

also  available  from  the  Ti  target. 

The  experimental  results  gave  no  observable 
centroid  shift  for  the  1.795  and  1.825  MeV  y-rays 
and  hence  upper  limits  of  <23  ps  and  <10  ps  were 
calculated,  respectively,  for  these  levels.  Weaver 
et  al  have  tentatively  set  a  value  of  spin  1/2  to 
the  1.795  MeV  level.  The  Weisskopf  E2  transition 
width  for  the  1.795  MeV  transition  to  the  ground  state 

_4 

spin  5/2  is  1.51  x 10  eV  and  hence  the  measured 
limits  on  the  lifetime  correspond  to  an  E2  strength 
of  0.19  W.u.  <  (E2)  <  2.17  W.u.  Weaver  et  al  assign  a 
spin  of  3/2+  to  the  1.825  MeV  level  with  mixing  ratios 
for  the  ground  state  transition  of  either  6 =  0 . 219± 0 . 037 
or  6  =  2.18^q*^.  Based  on  the  experimental  limits  for 

the  lifetimes  the  first  mixing  ratio  corresponds  to 

-5 

limits  on  the  El  and  M2  strengths  of  1.16  x  10  W.u.  < 

(El)  <5.96  x 10  W.u.  and  0.54  W.u.  <  (M2 )  <  5.22  W.u., 
respectively.  Similarly  for  6  =  2.18^*^  the  follow¬ 
ing  limits  are  obtained;  1.88  x 10  ^  W.u.  <  (El)  < 

1.21  xlO  W.u.  and  13.6  W.u .  <  (M2 )  <71.9  W.u.,  res¬ 
pectively.  The  M2  enhancements  obtained  for  6  =  2.18+q*^ 
are  rather  large  so  that  lifetime  considerations  indi¬ 
cate  that  6  =  0.219±0.037  is  the  more  acceptable  mixing 
ratio  for  this  transition. 


■ 
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49 

7.2c  ^Ti 

No  lifetime  information  at  all  has  been  reported 
49 

for  Ti.  The  level  scheme  of  fig.  56  is  that  given 

49  49 

in  Nuclear  Data  Sheets  (NDS  70).  Since  the  Ti(p,n)  V 

49 

Q  value  is  -1.384  MeV,  the  excited  states  of  Ti  are 

not  observed  using  proton  inelastic  scattering.  Hence 

a  4.0  MeV  pulsed  deuteron  beam  was  used  to  investigate 

49  47  48 

the  levels  in  Ti.  The  Ti  and  Ti  foils  were  used 

as  the  two  targets.  'Prompt'  reference  y-rays  are 

4  7  48  48  49 

provided  by  the  Ti(d,p)  Ti  and  the  Ti(d,n)  V 

49 

reactions  lifetime  measurements  in  V  are  discussed 

in  section  7.2e.  The  results  of  one  of  the  experimental 

runs  are  shown  in  fig.  56.  No  shifts  were  observed  for 

the  1.381,  1.585  and  1.762  MeV  y-ray  transitions  in 
49  .  . 

Ti.  Thus  upper  limits  for  the  lifetimes  of  these 
levels  are: 

1.381  MeV  level  <  5  ps 
1.585  MeV  level  <  11  ps 
1.762  MeV  level  <  14  ps. 

The  0.342  MeV  y-ray  transition  from  the  1.723  MeV 
level  to  the  first  excited  state  at  1.381  MeV  could  not 
be  studied  since  no  suitable  'prompt'  reference  y-rays 
could  be  found  in  this  energy  region  using  the  pulsed 
deuteron  beam. 

All  the  experimental  results  obtained  for  the 
Titanium  isotopes  are  shown  in  table  3  together  with 
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previous  results.  The  experimental  transition  widths 
are  also  shown  with  the  transition  strengths  in 
Weisskopf  units  for  both  electric  and  magnetic  y-ray 
transition  decay  modes. 

47 

7. 2d  *  'v 

Blasi  et  al  (B1  73)  investigated  the  low  lying 

47  47  47 

levels  in  V  through  the  Ti(p,n)  V  reaction  with 

proton  energies  between  3.6  and  5.3  MeV  and  through 
40  10  47 

the  Ca (  B,2pn)  V  reaction  at  a  beam  energy  of 

25  MeV.  Their  level  scheme  is  shown  in  fig.  57.  The 

4  7 

only  lifetime  measurements  on  the  levels  in  V  have 

been  performed  by  Blasi  et  al  .  (B1  73)  using  the  recoil 

distance  technique  to  study  the  levels  at  260  keV  and 

above  and  a  pulsed  proton  beam  to  obtain  upper  limits 

of  <lns  for  both  the  levels  at  88  and  146  keV.  Their 

results  are  also  shown  in  fig.  57. 

The  level  at  260  keV  was  studied  using  the 

^Ti  foil  and  the  ^Ni-^^Ta  foil  combination  as  the 

4  7  4  7 

second  target.  The  Ti(p,n)  V  reaction  has  a  Q 
value  of  -3.690  MeV.  An  example  of  the  experimental 
results  is  shown  in  fig.  54  of  section  7.2b.  The 
weighted  mean  of  six  experimental  runs  yielded  a 
lifetime  of  78±11  ps  for  the  260  keV  level  which  is 
in  agreement  with  the  result  of  90±12  ps  obtained  by 
Blasi  et  al.  The  spins  of  the  260  keV  level  and  ground 
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Fig.  57.  Experimental  level  scheme  and  previous 
lifetime  results  for  ^V. 
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state  are  3/2+  and  3/2  ,  respectively,  and  hence 

assuming  a  pure  El  decay  mode  the  present  result 

-4 

leads  to  an  El  transition  strength  of  5.46  x 10  W.u. 

The  levels  at  88  and  146  keV  were  studied 
using  as  the  second  target  the  P^PTa-P6^Ho  foil  com¬ 
bination  and  a  5.0  MeV  pulsed  proton  beam.  Typical 
experimental  results  are  shown  in  fig.  58.  The 
resulting  lifetime  of  the  146  keV  level  observed  by 
the  58  keV  y-ray  transition  to  the  88  keV  level  is 
786195  ps .  The  centroid  shift  of  the  88  keV  y-ray  is 
equivalent  to  1173±37  ps.  However,  the  contribution 
to  the  centroid  shift  by  the  cascading  58  keV  y-ray 
must  be  established  at  the  angle  of  observation  in 
order  to  determine  the  lifetime  of  the  88  keV  level. 

To  establish  the  contribution,  an  angular 
distribution  of  the  58  and  88  keV  y-rays  was  performed 
using  a  d.c.  proton  beam  at  an  energy  of  5.0  MeV. 

The  results  which  have  been  corrected  for  absorption 
in  the  walls  of  the  target  chamber  and  the  relative 
detector  efficiency  are  shown  in  fig.  59.  The  angular 
distributions  are  seen  to  be  nearly  isotropic. 


The  angular  distribution  of  nuclear  y-rays  can 
be  written  in  the  form  (Ro  67) 

W (0)  =  l 

k  even 


akPk^cos  0)  * 


7.1 
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Rose  and  Brink  (Ro  67)  have  shown  that,  if  only  two 
multipolarities  L  and  L'  are  present  in  the  decay, 


the  coefficients  a 


k 


can  be  written 


7.2 


where  is  the  spin  of  the  decaying  state,  the  spin 
of  the  final  state,  and  6^  is  the  mixing  ratio  of  the 
radiation.  The  R^  coefficients  are  "geometric'1  factors 
containing  vector  coupling  coefficients  and  are  deter¬ 
mined  completely  by  the  spins  and  multipolarities 
involved  in  the  decay.  The  B^(J^)  coefficient  contains 
the  physics  of  the  formation  of  the  state  J^,  and  is 
related  to  the  population  parameters  P(m^)  of  the  sub¬ 
states  of  J^. 


.T,  ,T  —  m 


7.3 


The  B^.(J-^)  coefficients  have  been  calculated  (Gr  71) 


using  the  Compound  Nuclear  Statistical  Model  (Sh  66) 
for  the  analysis  of  the  (p,ny)  reaction.  In  this 
model  the  substates  in  the  residual  nucleus  are  almost 
equally  populated  to  m^  =  J  (target)  +  1. 


If  the  state  of  spin  J ^  should  cascade  to  another 


state  spin  and  if  this  cascade  is  observed  then  the 
nuclear  alignment  parameters  B^(J^)  in  eq.  7.2  must  be 
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replaced  by  B^  ^JiJ2 )  where 


Uk(JlJ2)  =  IUk(I,JxJ2)  +  6X2Uk(L'JlJ2)  ]/(1  +612) 


7.4 


L,  L’  and  6-^  refer  to  the  cascade  from  state  spin 
to  state  spin  and  U^(LJ-^J2)  are  coefficients  con¬ 
taining  vector  coupling  coefficients  determined  com¬ 
pletely  by  the  spins  and  multipolarities  involved  in 
the  decay.  6-^  in  e<3*  7. 2  must  also  be  replaced  by 

<$23* 


Assuming  the  maximum  limit  for  E2  enhancements 
as  used  in  Nuclear  Data  Sheets  of  1000  VI. u.,  an  upper 
limit  for  the  mixing  ratio  of  |6|  <  0.08  is  obtained 
for  the  146  -+  88  keV  y-ray  transition  using  the  present 
result  for  the  lifetime  of  the  146  keV  level.  Using 
this  value  and  the  slope  of  the  58  keV  y-ray  angular 
distribution  to  obtain  a  value  of  [calculations  yield 
a^  ss  0],  eq.  7.2  gives  (J-^)  =  -0.076.  This  small 

value  of  (J^)  is  expected  since  the  target  nucleus 

4  7  - 

Ti  has  a  ground  state  spin  5/2  and  thus  equal  sub¬ 
state  populations  are  expected  up  to  m^  =  7/2  which 
is  the  spin  of  the  146  keV  state.  Equal  substate  popu¬ 
lations  result  in  an  isotropic  angular  distribution 
since  (J^ )  =  0 . 


' 
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Because  p  0,  the  angular  distribution  of 

the  88  keV  y-ray  should  also  be  isotropic  and  thus 
using  the  relative  aQ  coefficients  as  obtained  from 
fig.  59  the  lifetime  for  the  88  keV  level  was  calcu¬ 
lated  to  be  951±52  ps .  Menti  (Me  67)  has  done  Internal 

4  7 

Conversion  Measurements  on  V  and  according  to  him 
the  88  keV  state  decays  98%  of  the  time  via  the  Ml  mode 
to  the  ground  state.  Hence  assuming  a  pure  Ml  decay 
the  present  result  corresponds  to  a  Ml  strength  of 
0.05  W.u.  Again,  assuming  a  pure  Ml  decay  for  the 
146  +88  keV  y-ray  transition,  the  present  result  of 
786±95  ps  for  the  lifetime  of  the  146  keV  level  yields 
a  Ml  strength  of  0.21  W.u. 


49 

The  level  scheme  of  V  is  shown  in  fig.  60  from 

the  work  of  Blasi  et  al  (B1  71A) .  The  lifetimes  of  the 

levels  above  the  second  excited  state  have  been  studied 

using  the  Doppler  Shift  Attenuation  method  by  Kiuru 

4  8  4  9 

(Ki  72)  using  the  Ti(p,y)  V  reaction  at  a  maximum 

beam  energy  of  1.57  MeV  and  recently  by  Tabor  and 

46  49 

Zurmiihle  (Ta  74)  using  the  Ti(a,p)  V  reaction  at 
a  beam  energy  of  10  MeV.  Their  results  are  shown  in 
fig.  60.  The  lifetimes  measured  by  Tabor  and  Zurmiihle 
for  the  748  and  1140  keV  levels  are  substantially  larger 
than  those  reported  by  Kiuru.  However,  most  of  the  other 
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Fig.  60.  Experimental  level  scheme  and  previous  lifetime 
results  for  states  in  49V. 
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lifetimes  are  in  reasonable  agreement.  Brown  et  al 

(Br  74)  using  the  recoil  distance  technique  and  the 
46  49 

Ti(a,p)  V  reaction  at  a  beam  energy  of  11  MeV 
obtained  a  value  of  5.1±1.0  ps  for  the  1021  keV  level. 

The  lifetime  of  the  level  at  153  keV  has  been 
measured  using  several  techniques.  Currie  (Cu  63) 
using  a  3+-y  delayed  coincidence  technique  obtained 
a  value  of  29.7±0.8  ns.  Blasi  et  al  (B1  71A)  using 
a  pulsed  proton  beam  obtained  a  lifetime  of  28.9±2.9  ns 
and  finally  Vingiani  et  al  (Vi  72)  using  a  y-y  delayed 
coincidence  technique  obtained  a  value  of  28.7±0.5  ns. 

Two  measurements  of  the  lifetime  of  the  first 
excited  state  at  90  keV  have  been  reported  both  using 
the  3+-y  delayed  coincidence  technique.  Cheung  and  Mark 
(Ch  71)  obtained  the  value  of  620±29  ps  while  Okon  et  al 
(Ok  73)  measured  the  lifetime  to  be  476±29  ps.  These 
results  are  seen  to  be  in  disagreement  and  a  remeasure¬ 
ment  of  the  lifetime  of  the  90  keV  level  would  provide 

more  information  about  this  state. 

49 

The  Ti  foil  was  used  to  study  the  levels  in 
49  49  49 

V  via  the  Ti(p,n)  V  reaction  which  has  a  Q  value 
of  -1.384  MeV.  Supplementary  'prompt’  reference  y-rays 
were  obtained  using  as  the  second  targets  the  ^^Ta-^^Ho 
foil  combination  and  the  Ni-  xTa  foil  combination. 

Due  to  the  long  lifetime  of  the  153  keV  level  the  Slope 
Method  was  used  to  calculate  the  lifetime.  Fig.  61  shows 
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Fig 


SINTOO 


.  61.  Logarithmic  plot  of  the  delayed  time  distribu¬ 
tion  of  the  153  keV  -*g.s.  transition  in  ^V. 
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a  logarithmic  plot  (base  10)  of  the  delayed  time 
distribution  obtained  from  the  153  keV  +  g.s.  decay. 

A  background  subtraction  has  been  performed  on  the 
result  shown  in  fig.  61.  A  least  squares  fit  to  the 
slope  from  seven  experimental  runs  yielded  a  lifetime 
of  28.71±0.35  ns.  This  result  is  in  excellent  agree¬ 
ment  with  the  previous  results. 

The  153  keV  transition  is  expected  to  be  solely 
E2  in  character  since  the  transition  proceeds  from  the 
3/2  153  keV  level  to  the  7/2  spin  ground  state  via 

a  58%  branch.  The  experimental  transition  width  for 

_  8 

this  decay  is  1.33  x  10  eV  and  so  our  experimental 
determination  of  the  lifetime  corresponds  to  an  E2 
strength  of  18.5  W.u. 

The  calculation  of  the  lifetime  of  the  90  keV 
level  is  complicated  by  the  42%  branch  from  the  153  keV 
level  and  thus  the  delayed  time  distribution  for  the 
90  keV  y-ray  transition  has  a  tail  as  illustrated  in 
fig.  62.  The  analysis  entailed  first  the  subtraction 
of  the  background  contributions  from  the  time  distri¬ 
butions  of  the  90  and  153  keV  y-rays.  Using  the  delayed 
time  distribution  of  the  153  keV  y-ray,  a  fit  was  made 
to  the  tail  of  the  time  distribution  for  the  90  keV  y- 
ray  and  a  subtraction  performed.  Account  was  taken  of 
the  walk  with  energy  of  the  prompt  time  centroid  curve 
between  the  two  energies,  90  and  153  keV.  The  resulting 
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time  distribution  for  the  90  keV  y-ray  is  shown  in 
fig.  62.  Since  the  tail  contribution  is  relatively 
small  compared  to  the  time  contribution  to  the  90  keV 
y-ray  decay  from  direct  population  of  the  state  by 
the  (p,n)  reaction,  a  100%  error  in  the  determination 
of  the  walk  of  the  prompt  time  centroid  curve  would 
only  lead  to  an  error  of  15  ps  in  the  final  result. 

Another  point  considered  was  the  effect  of  the 
small  change  in  FWHM  between  prompt  time  distributions 
at  90  and  153  keV.  An  investigation  of  this  effect 
indicated  that  a  maximum  possible  error  of  20  ps  could 
occur. 

An  example  of  the  experimental  results  using  a 
5.0  MeV  pulsed  proton  beam  and  the  two  targets,  the 
^Ti  foil  and  the  ^^Ta-^^Ho  foil  combination,  are 
shown  in  fig.  63.  The  weighted  mean  of  six  experimental 
runs  yielded  a  lifetime  for  the  90  keV  level  of  329±19  ps. 
This  result  does  not  agree  with  either  of  the  previous 
results . 

The  total  internal  conversion  coefficient  for 
the  90  keV  transition  has  been  measured  by  Menti  (Me  67) 
to  be  0.034.  According  to  him  the  90  keV  state  decays 
more  than  99%  of  the  time  via  the  Ml  mode  to  the  ground 
state.  Thus  assuming  a  pure  Ml  transition  the  Ml  strength 
for  this  transition  is  0.13  W.u. 
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Very  little  experimental  work  has  been  reported 
50 

for  V.  However,  Blasi  et  al  (B1  69A)  studied  the 
decay  and  level  scheme  of  the  low  energy  states  and 
their  results  are  shown  in  the  insert  to  fig.  64. 

They  found  that  there  was  no  direct  feeding  of  the 
levels  at  226  and  320  keV  using  the  (p,n)  reaction  at 
proton  energies  below  4.5  MeV  and  that  these  levels  are 
fed  by  cascading  y-rays  from  the  levels  at  355  and 
388  keV. 

Both  Temmer  and  Heydenburg  (Te  56)  and  Fagg  et 

al  (Fa  56)  obtained  by  Coulomb  Excitation  measurements 

a  value  for  the  reduced  transition  probability  B(E2)f 
2  4 

of  110  e  fm  for  the  226  keV  level  which  is  equivalent 
to  an  E2  partial  lifetime  of  t(E2)  =  10.6  ns.  Blasi 
et  al  (B1  69A)  using  a  pulsed  beam  gave  an  upper  limit 
of  3  ns  to  this  level,  which  suggests  a  strong  magnetic 
dipole  contribution. 

Experimental  investigations  on  the  low  lying 

50  2 

energy  levels  of  V  were  performed  using  the  1  mg/cm 

50  64 

enriched  Ti  foil  and  two  second  targets:  The  Ni 

foil  to  study  the  226  keV  Y-ray  transition  and  the 
l^Ta-^^Ho  foil  combination  to  study  the  94  keV  tran¬ 
sition.  Pulsed  proton  beams  of  4  to  4.4  MeV  were  used 

50  50  50 

to  excite  the  states  in  V  via  the  Ti(p,n)  V 
reaction  which  has  a  Q  value  of  -2.997  MeV.  Typical 
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results  are  shown  in  fig.  64.  Within  the  error  bars 
an  identical  centroid  shift  was  observed  for  the 
94  keV  and  226  keV  transitions  in  all  experimental 
runs.  Since  the  226  keV  level  is  fed  entirely  by  the 
94  keV  y-ray  cascade  from  the  320  keV  level  the  life¬ 
time  of  the  226  keV  level  must  be  short.  The  present 
result  is  t  < 16  ps  which  is  the  calculated  error  in  the 
difference  of  the  weighted  means  of  the  centroid  shifts 
for  the  94  and  226  keV  transitions. 

The  320  keV  level  is  also  fed  entirely  from  y-ray 
cascades  from  the  355  keV  level  and  hence  only  an  upper 
limit  of  <450  ps  for  lifetime  of  the  320  keV  level  could 
be  calculated. 

The  Weisskopf  estimate  of  the  Ml  partial  life¬ 
time  for  the  226  keV  transition  is  x  (Ml)  =  2.7  ps  and 

03 

hence  our  result  of  <16  ps  indicates  a  strong  magnetic 
dipole  contribution  to  this  transition.  The  Weisskopf 
estimates  of  the  Ml  partial  lifetimes  for  the  94  keV 
transition  from  the  320  keV  level  and  the  35  keV  tran¬ 
sition  from  the  355  keV  level  are  x  (Ml)  =  38  ps  and 

03 

x  (Ml)  =  741  ps ,  respectively.  Our  upper  limit  of 

03 

<450  ps  for  the  lifetime  of  the  320  keV  level  suggests 
that  the  apparent  lifetime  is  mostly  due  to  the  feeding 
of  the  320  keV  level  by  the  35  keV  y-ray  transition 


from  the  355  keV  state. 
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7 . 2g  51V 
51 

V  has  received  extensive  experimental  and 
theoretical  study  due  to  its  closed  shell  of  28 
neutrons.  The  level  scheme  shown  in  fig.  65  is  from 
the  work  of  Goodman  and  Donahue  (Go  72)  and  Horoshko 
et  al  (Ho  70).  Also  shown  are  the  lifetimes  obtained 
by  Goodman  and  Donahue  using  the  Doppler  Shift  Attenua¬ 
tion  technique.  Only  one  measurement  has  been  reported 
for  the  second  excited  state  at  928  keV,  that  by  Singh 
et  al  (Si  70),  who  obtained  a  lifetime  of  101±36  ps 

using  the  delayed  coincidence  technique  and  recording 

51 

the  delays  between  the  8  rays  from  the  5-8  min  Ti 

51 

feeding  levels  in  V  and  the  following  y-rays.  They 
also  obtained  a  lifetime  of  274±43  ps  for  the  first 
excited  state  at  320  keV.  Many  other  measurements  on 
the  lifetime  of  the  first  excited  state  have  been 
reported  and  the  results  are  shown  in  the  insert  to 
fig.  65.  The  weighted  average  of  these  results  is 
290 ±20  ps. 

51 

For  our  experiments  a  gloop  V  target  was 
64 

prepared  and  the  Ni  foil  was  used  as  the  second 
target.  One  set  of  results  is  shown  in  fig.  66.  A 
weighted  mean  of  three  experimental  runs  yielded  a 
lifetime  for  the  320  keV  level  of  266±9  ps  which  is  in 
agreement  with  the  previous  results. 
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Fig.  65.  Experimental  level  scheme  and  previous 
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Fig.  66. 


Experimental  result 
of  the  lifetime  of 


s  for  tne  determination 
the  32  0  keV  level  in  ?'^V. 


The  solid  line  is  a  theoretical  least  squares 
fit  to  the  experimental  centroid  points. 
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Horoshko  et  al  (Ho  70)  have  calculated  the 

weighted  mean  of  the  reduced  transition  probability 

B(E2)1  from  several  Coulomb  Excitation  measurements 

2  4 

and  obtain  a  value  of  154.0±7.6  e  fm  for  the  320  keV 

transition.  Thus  using  this  result  and  our  lifetime 

result  a  value  for  the  mixing  ratio  of  |6|=  0.45±0.02 

is  obtained.  The  value  is  in  agreement  with  the  two 

experimental  determinations  of  the  mixing  ratio; 

Ritter  et  al  (Ri  62)  obtained  6  =  -0.52±0.07  from 

20 

Coulomb  Excitation  measurements  using  18  MeV  Ne  ions, 
and  Krause  (Kr  63)  obtained  6  =  -0.43±0.03  from 
Resonance  Fluorescence  measurements.  The  experimental 

E2  and  Ml  transitions  widths  are  T(E2)  =  4.16x10”^  eV 

—  6 

and  T (Ml)  =  2.04x10  eV  and  hence  our  result  corresponds 

_3 

to  strengths  of  13.7  W.u.  and  3.0x10  W.u. ,  respectively. 

The  928  keV  level  was  not  observed  from  proton 

inelastic  scattering  so  this  state  was  investigated 

50 

using  a  pulsed  deuteron  beam  incident  on  the  Ti  foil. 

Additional  'prompt*  reference  y-rays  were  obtained 

4  7 

using  a  second  target,  the  Ti  foil,  and  the  reaction 
4  7  4  8 

Ti (d,p)  Ti.  An  illustration  of  the  results  is  shown 
in  fig.  67.  No  observable  centroid  shift  was  found  for 
the  928  keV  y-ray  decay  and  an  upper  limit  of  x  <  14  ps 
was  calculated  for  the  lifetime.  This  result  is  in 
marked  disagreement  with  that  obtained  by  Singh  et  al 
(Si  70)  who  obtained  a  value  of  101±36  ps . 
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Fig.  67. 


Experimental  results  for  tlie  determination 
of  the  lifetime  of  the  928  keV  level  in  ^ly. 
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Horoshko  et  al  (Ho  70)  have  measured  the  reduced 
transition  probability  using  the  Coulomb  Excitation 
Method  to  be  B(E2)I  =  76±5  e^fm^  for  the  928  keV  -*  g .  s . 
transition.  They  also  obtained  a  mixing  ratio  of 
6  =  0.03±0.09  for  this  transition  which  indicates  a 
pure  E2  transition.  Their  measured  branching  ratio 
was  86%.  Thus  the  calculated  lifetime  from  these  mea¬ 
surements  is  x  =  13.4±1.0  ps  which  agrees  with  our  upper 
limit . 

For  the  3/2  928  keV  to  the  5/2  320  keV  level 

y-ray  transition  Horoshko  et  al  obtained  two  possible 

mixing  ratios  of  6  =  -0.33+^’^  and  6  =  -8.8+^*''"  which 

are  in  qualitative  agreement  with  the  results  of 

Robinson  et  al  (Ro  63)  who  obtained  the  values  of  6  = 

-0.18±0.02  and  6  =  +9±2.  From  theoretical  considerations 

+4  1 

Horoshko  et  al  strongly  favour  the  value  of  6  =  -8.8 

—  oo 

which  indicates  nearly  a  pure  E2  transition.  They  also 

found  agreement  between  the  measured  and  calculated 

+4  1 

B(E2)  values  only  for  6  =  -8.8  *  .  Using  this  value 

—  00 

2  4 

and  the  B(E2)+  value  of  107±9  e  fm  calculated  by 

Horoshko  et  al,  the  calculated  value  for  the  lifetime 
+3  3 

is  t  =  12.7_2*q  ps ,  again  in  agreement  with  our  results. 

Table  4  shows  all  the  experimental  results 
obtained  for  the  Vanadium  isotopes  and  a  comparison 
with  previous  results.  Also  included  in  table  4  are 
the  experiment  transition  widths  and  the  transition 
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strengths  in  Weisskopf  units  for  both  magnetic  and 
electric  decay  modes. 

n  59,..  -.65,. 

7 . 3  Ni  and  Zn 

These  two  isotopes  were  initially  investigated 
to  see  whether  they  could  be  used  as  sources  of  'prompt' 
reference  y-rays.  However,  both  isotopes  yielded  low 
energy  y-ray  transitions  from  states  whose  lifetimes 
were  measurable  by  our  technique.  The  level  schemes  for 

C  Q  C  C 

Ni  and  Zn  from  the  work  by  Hutton  et  al  (Hu  73)  and 
Ezell  and  Scott  (Ez  74),  respectively,  are  shown  in 
fig.  68. 

_  _  59... 

7 . 3a  Ni 

59 

A  gloop  Co  target  was  prepared  and  the  second 

64 

target  used  was  the  Ni  foil.  A  4.4  MeV  pulsed  proton 

59  .  64 

beam  was  used  to  excite  the  states  in  Ni  and  Cu 

by  the  (p,n)  reaction,  the  Q  value  being  -1.858  MeV 
59  59 

for  the  Co(p,n)  Ni  reaction.  Results  obtained  in 
one  of  the  experimental  runs  are  shown  in  fig.  69  which 
clearly  show  the  centroid  shifts  of  the  340  and  465  keV 
y-ray  transitions.  The  weighted  mean  of  four  experimen¬ 
tal  runs  gave  lifetime  results  for  the  340  and  465  keV 
levels  of  93±10  ps  and  35±18  ps,  respectively. 
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Fig.  68.  Experimental  level  schemes  for  9Ni  and  5Zn 
and  previous  lifetime  results  for  excited 
states  of  ^^Ni. 
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Fig.  69.  Experimental  results  for  lifetime  determina¬ 
tions  in  jyNi. 
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Hutton  et  al  (Hu  73)  studied  the  lifetimes  of 

59  56  59 

Ni  using  the  Fe(a,n)  Ni  reaction  (E  =  8.5-10.5 

a 

MeV)  and  the  Doppler  Shift  Attenuation  and  Recoil 
Distance  techniques.  They  obtained  a  lifetime  of 
120±20  ps  for  the  340  keV  level  and  29±5  ps  for  the 
465  keV  level.  Within  the  experimental  errors  these 
results  are  in  very  good  agreement  with  the  present 
results.  Both  the  levels  at  340  keV  and  465  keV  are 
populated  by  cascading  y-rays  from  higher  excited  states. 
However,  these  higher  energy  levels  have  lifetimes  less 
than  2  ps  as  illustrated  in  fig.  68  from  the  results 
obtained  by  Hutton  et  al  (Hu  73) ,  and  therefore  they 
have  no  effect  on  the  measured  lifetime  values. 

The  mixing  ratio  of  6  =  -0.11±0.03  was  obtained 
by  Hutton  et  al  for  the  340  keV  y-ray  transition. 

Using  our  result  for  the  lifetime,  the  experimental 
transition  widths  for  the  Ml  and  E2  modes  of  decay 
are  T(Ml)  =  6.99xl0~6  eV  and  T(E2)  =  8.46xl0-8  eV 
which  correspond  to  a  Ml  and  E2  strength  of  0.009  W.u. 
and  1.7  W.u.,  respectively. 

The  mixing  ratio  has  not  been  measured  for  the 
spin  1/2  465  keV  level.  If  a  pure  Ml  transition  is 

assumed  a  Ml  strength  of  0.009  W.u.  is  obtained. 

Taking  the  other  extreme  of  a  pure  E2  transition  an 
E2  strength  of  78  W.u.  is  obtained.  This  E2  enhancement 
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is  of  much  greater  magnitude  than  generally  found  and 
thus  this  transition  is  expected  to  be  mostly  Ml  in 
character. 

7.3b  65Zn 

The  lifetime  of  the  first  excited  state  at 
54  keV  has  been  measured  by  August  and  Friichtenicht 
(Au  60)  to  be  1.65±0.05  ys  using  a  3+-Y  delayed 
coincidence  technique.  Using  a  delayed  coincidence 

technique  Li-Scholz  and  Bakhru  (Li  68)  obtained  a  life¬ 
time  for  the  115  keV  level  of  663±101  ps .  The  level  at 
207  keV  has  been  assigned  an  upper  limit  of  500  ps  by 
Tubbs  (Tu  66)  and  hence  this  level  was  a  candidate  for 
our  technique. 

The  target  used  was  a  natural  copper  foil, 

2 

48  mg/cm  ,  which  has  a  natural  isotopic  abundance  of 

30.9%  ^Cu,  the  rest  being  ^Cu.  A  4  MeV  pulsed 

65 

proton  beam  was  used  to  excite  the  states  of  Zn  via 
6  3  6  3 

the  DCu  (p  ,n)  D:jZn  reaction  which  has  a  Q  value  of 

/TO  /T  O 

-2.13  MeV.  The  Q  value  of  the  Cu(p,n)  Zn  reaction 

6  3 

is  -4.14  MeV  so  that  no  Y~raY  transition  in  Zn  are 
observed.  The  second  target  used  was  either  the 
^Ni-^^Ta  foil  combination  or  the  ‘*'^'*'Ta-^^^Ho  foil 
combination.  Two  examples  of  the  experimental  results 
are  shown  in  figs.  70  and  71. 
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An  unexpected  centroid  shift  was  observed  for  the 
201  keV  y-ray  decay  from  the  9/2+  1066  keV  level  to 
the  7/2  state  at  865  keV.  The  two  decays  at  202  and 
212  keV  in  64Cu  provide  very  good  'prompt'  references 
for  the  201  keV  transition.  The  lifetime  was  found  to 
be  829±37  ps  for  the  1066  keV  level.  Ezell  and  Scott 
(Ez  74)  have  assumed  a  pure  El  decay  for  this  transition 

which  has  an  experimental  transition  width  of  T  = 

-7 

7.94x10  eV  using  our  lifetime  result.  This  corresponds 

-5 

to  an  El  strength  of  9  x 10  W.u.  Assuming  a  pure  M2 

4 

transition  a  strength  of  10  W.u.  is  calculated.  Enhance¬ 
ments  of  this  magnitude  are  not  found  so  that  an  El 
decay  mode  is  predominant  for  the  201  keV  y-ray  tran¬ 
sition  . 

The  state  at  207  keV  decays  by  two  branches;  one 
via  a  153  keV  transition  to  the  first  excited  state  at 
54  keV  and  a  second  to  the  ground  state.  The  weighted 
mean  of  four  experimental  runs  yielded  lifetimes  of 
213±11  ps  for  the  153  keV  transition  and  227±22  ps  for 
the  207  keV  transition.  The  adopted  lifetime  for  the 
207  keV  level  was  then  calculated  as  the  weighted  mean 
of  these  two  results  which  gives  216±10  ps.  Robert  et 
al  (Ro  70A)  have  studied  the  low  energy  states  of  ^~>Zn 
and  deduce  the  mixing  ratio  of  the  153  keV  transition 
to  be  6  =  0.24±0.06.  Using  this  result  and  their 
branching  ratios  for  the  207  keV  state,  our  lifetime 
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result  leads  to  experimental  transition  widths  of 
r(Ml)  =  2.19xl0"4 * 6  eV  and  T(E2)  =  (1.2610.65)  x 10"7  eV, 
respectively,  for  the  153  keV  y-ray  transition.  Hence 
the  Ml  and  E2  strengths  are  0.03  W.u.  and  120±62  W.u., 
respectively . 

The  result  for  the  lifetime  of  the  115  keV  level 
was  calculated  to  be  641113  ps.  The  previous  result 
by  Li-Scholz  and  Bakhru  (Li  67)  of  663+101  ps  is  in 
agreement  with  the  present  result.  Assuming  a  pure 
Ml  decay  mode  for  this  transition  a  Ml  strength  of 
0.03  W.u.  is  obtained. 

The  above  levels  are  fed  in  part  from  cascading 

65 

y-rays  from  higher  energy  levels  in  Zn.  Thus  a 
study  of  the  higher  excitation  energy  levels  was  per¬ 
formed.  Two  second  targets  were  employed  to  provide 

the  necessary  'prompt'  reference  y-rays.  These  were 
4  7  4  9 

the  Ti  foil  and  the  Ti  foil.  The  lifetimes  of  the 

4  7  49 

levels  in  Ti  and  V  are  discussed  in  sections  7.2b 

and  7.2e,  respectively.  Using  the  natural  copper  foil 

to  obtain  the  reaction  ^Cu (p,n) ^~*Zn,  the  reaction 

6  3  6  3 

JCu  (p,p  '  )  D  JCu  also  provides  'prompt'  reference  y-rays 
at  energies  669,  962  and  1327  keV  (Wo  70).  An  illus¬ 
tration  of  the  experimental  results  is  shown  in  fig.  72 

4  7  49 

for  the  case  of  the  copper  foil  and  the  '  Ti  target. 

No  observable  shifts  were  found  except  for  the 
865  keV  y-ray  transition  from  the  865  keV  level.  This 
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Fig.  72. 


Experimental  results  for  lifetime  determinations 
of  the  higher  energy  levels  in  ^5zn. 


ENERGY  (MeV) 


level  is  fed  by  the  201  keV  y-rays  from  the  1066  keV 
level  and  it  was  found  that  the  shift  was  entirely 
due  to  the  feeding  of  the  865  keV  level  by  cascading 
y-rays  from  the  1066  keV  level. 

Table  5  shows  the  upper  limits  which  have  been 

assigned  to  the  higher  energy  levels  in  ^Zn  plus  the 

other  results  obtained  for  the  low  energy  levels  of 
59  65 

Ni  and  Zn.  Previous  results  are  also  given  where 
applicable  and  the  table  also  includes  the  experimen¬ 
tal  transition  widths  and  the  strengths  in  Weisskopf 
units  for  both  electric  and  magnetic  decay  modes. 
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CHAPTER  8 


THEORETICAL  INTERPRETATIONS 

7 

The  1  f  shell  region  has  been  considered  with 
particular  interest  since,  following  the  shell  model 
picture,  nucleons  are  filling  a  fairly  well  isolated 
single  particle  level  of  quite  large  j  and  thus  several 
simple  shell  model  states  are  allowed.  This  argument 

was  based  on  the  relatively  large  energy  gap  between 

7  .  3 

the  1 f  ^  and  the  neighbouring  2s-ld  and  2  p  ^  shells, 

40  4  8 

and  by  the  presence  of  the  doubly  magic  Ca  and  Ca 
cores  which  have  been  found  to  be  particularly  stable. 

7 

The  properties  of  the  nuclei  in  the  If  j  shell 

have  been  studied  by  McCullen,  Bayman  and  Zamick  (Me  64) 

[MBZ  model]  in  terms  of  the  spherical  shell  model  with 

a  suitable  residual  interaction  between  nucleons  in  a 

pure  (If  ~)n  configuration  outside  the  closed  shells. 

They  assumed  that  the  residual  interaction  was  a  two- 

body  interaction  where  the  matrix  elements  are  completely 

determined  by  the  interaction  energy  differences  between 

7  ^ 

the  eight  states  (If  ^  ,  I  =  0,1,...  7.  They  presumed 

that  these  energy  differences  describe  the  observed 

42  42  50  . 

states  in  Ca,  Sc  and  Ti. 

Though  many  of  the  observed  properties  are  accoun¬ 
ted  for  by  this  model,  there  are  serious  discrepancies 
which  cannot  be  avoided  without  more  realistic 
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181 


assumptions.  One  test  of  the  pure  configuration 
assumption  is  a  comparison  of  the  spectra  of  cross¬ 
conjugate  pairs  of  nuclei,  that  is  nuclei  with  Z=20+a, 

N  =  20+b  and  Z  =  28-b  and  N  =  28-a.  From  particle- 

4  7 

hole  theorems  these  cross-conjugate  pairs,  e.g.  Ti 
49 

and  V,  should  exhibit  identical  spectra.  However, 
experimental  level  schemes  are  not  identical.  Equiva¬ 
lent  nuclei  which  are  nuclei  having  the  same  number  of 

43  45  51 

particles  and  holes  i.e.  Ca,  Ca  and  V,  are  also 

predicted  to  have  approximately  the  same  level  spectrum 

but  experimentally  this  is  not  so. 

Probably  the  most  serious  weaknesses  of  the  model 

are  its  inability  to  explain  the  ground  state  triplets 

45  .  47  49 

exhibited  by  Ti ,  V  and  V  and  that  Ml  transitions 

between  states  are  forbidden  although  they  occur 

between  experimental  levels.  The  inclusion  of  other 

5  3  1 

configurations,  including  the  If  j,  2p  ^  and  2  p  j  con¬ 
figurations,  is  expected  to  be  important  even  in  low 
lying  states  of  nuclei  in  this  region.  Electromagnetic 
multipole  matrix  elements,  particularly  the  B(M1)  values, 
could  be  strongly  affected  by  such  admixtures.  However, 
such  configuration  mixing  has  only  been  attempted  with 
the  Ca  isotopes  (Mi  61,  Ra  65,  En  66)  and  the  N  =  28 
isotones  (Au  67,  Li  70,  Li  71).  Lips  and  McEllistrem 

(Li  70,  Li  71)  used  a  model  which  permits  protons  out- 

4  8 

side  the  closed  Ca  core  to  occupy  states  of  the 
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7  n  7  n-1  _  n-1 

configurations  (If  ,  (i  f  ^  ip  j, )  and  (lfj 

5 

If  |)  using  a  modified  effective  interaction  calcula¬ 
tion  to  obtain  the  eigenstates. 

These  calculations  yielded  appreciable  improve¬ 
ments  over  the  pure  configuration  model.  However, 
although  further  admixtures  of  different  configurations 
may  explain  other  discrepancies,  extensive  computation 
is  required  as  well  as  a  large  number  of  parameters. 

An  increasing  amount  of  evidence  indicates  that  strong 

7 

collective  behaviour  occurs  in  nuclei  m  the  1 f  ^  shell, 

so  to  overcome  the  problems  of  the  shell  model,  Malik 

and  Scholz  (Ma  66,  Sc  66,  Sc  67)  tried  an  alternative 

approach  to  the  theoretical  description  of  these  nuclei. 

They  used  the  Strong  Coupling  Symmetric  Rotator  model  (RPC) 

including  the  Coriolis  coupling  between  rotational  bands. 

Many  of  the  difficulties  encountered  in  the  shell 

model  are  avoided  in  this  treatment,  namely  the  repro- 

45  47 

duction  of  the  ground  state  triplets  in  Ti ,  V  and 
49 

V,  the  correct  number  of  levels  below  2.5  MeV  for 

7 

nuclei  in  the  upper  half  of  the  If  ^  shell,  and  a  better 
prediction  of  the  spectra  for  cross-conjugate  nuclei. 

They  found  that  the  Coriolis  coupling  strongly  mixed 
the  different  bands  so  that  the  original  band  structure 
and  the  spacings  of  the  single  particle  Nilsson  levels 
is  destroyed.  Also  the  level  schemes  of  the  even-even 
adjacent  nuclei  of  Titanium  and  Chromium  (NDS  73A)  as 
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shown  in  fig.  73  do  not  exhibit  the  characteristic 

rotational  J  (J+l)  type  spectra  with  spin  values 

0+  2+  4+  ....  However,  band  mixing  could  destroy 

the  rotational  type  spectra  in  these  nuclei. 

Recently  evidence  has  been  established  that 

7 

many  of  the  odd-even  nuclei  in  the  1 f  j  shell  have 
several  low  lying  positive  parity  states  which  were 
considered  as  hole  states  by  the  excitation  of  a  par¬ 
ticle  out  of  the  Id -2s  shell.  Blasi  et  al  (B1  70, 

Bl  71)  pointed  out  that  the  positive  parity  states 
4  5  4  5 

of  Ti  and  Sc  show  the  characteristic  J  (J+l)  sequence 

3 

of  a  rotational  band  built  on  a  d-^-  hole  state.  Similar 

47  49 

rotational  effects  have  been  observed  in  V  and  V 
(Bl  71A,  Bl  73) .  Fig.  74  shows  the  characteristic 
J (J+l)  sequence  for  these  nuclei  where  the  energy  of 
the  positive  parity  levels  are  plotted  against  J(J+1). 

Due  to  this  strong  evidence,  the  Symmetric 
Rotator  model,  as  described  in  Appendix  B,  was  used  to 
compare  the  level  schemes  and  electromagnetic  transition 
rates  of  the  odd-even  isotopes  of  Titanium  and  Vanadium. 
A  computer  program  written  by  Dr.  D.A.  Hutcheon  at 
Oxford  University  was  modified  and  corrected  by  myself 
and  Dr.  Hutcheon  to  run  on  the  University  of  Alberta's 
IBM  360  computer.  The  model  has  several  parameters, 
however,  an  attempt  was  made  to  develop  a  'prescription ' 
to  obtain  the  parameters  of  the  model  from  experimental 
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Fig.  73.  Experimental  level  schemes  of  the  even-even 
Titanium  and  Chromium  isotopes . 
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The  excitation  energies  of  the  positive 
parity  levels  in  ^-*Ti,  ^Sc,  4  7 y  and  49 v 


plotted  as  a  function  of  J  (J+l)  . 


Fig.  74. 
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data.  Thus  none  of  the  parameters  were  allowed  to 

vary,  except  the  deformation  parameter,  6.  Hence  it 

was  felt  that  a  meaningful  test  of  the  usefulness  and 

validity  of  the  model  for  reproducing  the  properties 

7 

of  the  nuclei  m  the  1 f  ^  region,  could  be  established. 

Values  of  k  =  0.066  and  y  =  0.32  were  obtained 
4  8  3  49 

from  the  Ca (  He,d)  Sc  data  (NDS  70)  by  calculating 
the  unperturbed  single-particle  energies  E  .  given  by 


Zj 


the  centre-of-gravity  relation 

l  E.  S. 

4*  i  l 


E£j  “  l  S± 


i 

where  denote  the  excitation  energy  of  a  level  belong¬ 
ing  to  the  Z j  shell  model  states  and  are  the  corres¬ 
ponding  spectroscopic  factors.  These  values  are 
similar  to  those  originally  used  by  Nilsson  (Ni  55)  of 
x  =  0.05  and  y  =  0.35.  The  moment  of  inertia  parameter 
fi  /2^  for  each  nucleus  studied  was  calculated  from  the 
average  energies  of  the  first  2+  state  of  the  adjacent 
even-even  nuclei,  assuming  a  rotational  character  to 
these  levels.  For  each  nucleus  only  one  rotational 
constant  was  used  for  all  bands  to  eliminate  further 
parametrization . 

The  computer  program  allows  particle  and  hole 
excited  states  of  the  Nilsson  model  by  permitting  an 
unpaired  nucleon  to  occupy  any  orbit  in  the  N  =  3 
oscillator  shell,  with  the  remaining  pairs  occupying 
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the  lowest  available  orbitals.  For  the  band  mixing 
calculations,  all  bands  based  on  the  ten  available 
single  particle  and  hole  states  in  the  N  =  3  shell 
have  been  included.  However,  overlap  between  bands 
based  on  the  particle  excited  states  and  the  hole 
excited  states  are  not  allowed. 

The  predicted  spectra  for  the  odd-even  isotopes 
of  Titanium  and  Vanadium  are  shown  in  figs.  75  to  80. 

Also  shown  in  each  figure  are  the  corresponding  experi¬ 
mental  level  schemes,  only  negative  parity  states 

being  shown,  and  the  theoretical  predictions  of  the 
7  n 

pure  (If  j)  configuration  model,  the  MBZ  model  (Me  64) , 

where  calculations  have  been  done.  The  predicted 
45  47  49  51 

spectra  for  Ti,  V,  V  and  V,  each  with  three 

nucleons  of  the  odd  species  outside  the  closed  shell, 

are  in  good  agreement  with  the  experimental  results. 

51  - 

However,  in  V  the  first  excited  5/2  state  is  predic¬ 
ted  too  high  and  the  9/2  and  11/2  spin  states  are 

45  47  49 

predicted  slightly  low.  In  Ti ,  V  and  V,  the 

ground  state  triplets  are  very  well  reproduced  by  the 

theory  with  the  correct  spin  assignments.  The  pure 
7  n 

(1  f  2-)  configuration  model  is  seen  only  to  predict 

a  doublet  for  these  nuclei  with  incorrect  spin  assignments 

45  47 

for  the  ground  states  of  Ti  and  V.  The  Strong 
Coupling  model  is  also  seen  to  reproduce  the  large 
energy  gap  between  the  ground  state  triplets  and  the 
other  excited  states. 
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Fig.  75.  Experimental  and  theoretical  level  schemes 
of  45Ti. 
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Fig.  76. 


Experimental  and  theoretical  level  schemes 
of  ^Ti. 
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Fig.  77. 


Experimental  and  theoretical  level  schemes 
of  49Ti. 
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Fig.  78. 


Experimental  and  theoretical  level  schemes 
of  47V. 
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Fig.  79. 


Experimental  and  theoretical  level  schemes 
of  49V. 
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Fig.  80.  Experimental  and  theoretical  level  schemes 
of  51V. 
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The  Strong  Coupling  model  produces  only  a  fair 
fit  to  the  experimental  level  scheme  of  47Ti.  How¬ 
ever,  the  anomalous  5/2  ground  state  spin  is  predic¬ 
ted  plus  the  low  lying  7/2  first  excited  state.  The 
9/2  second  excited  state  as  well  as  the  3/2*"  state 
are  predicted  too  low.  An  extra  state,  spin  1/2**,  is 
also  predicted  at  558  keV  which  is  not  observed  in  the 

experimental  level  scheme. 

49 

For  the  case  of  Ti ,  where  there  is  little 
experimental  data,  the  ground  state  spin  of  l/2~  is 
reproduced  with  the  large  energy  gap  to  the  first 
excited  state. 

The  electromagnetic  transition  rates  predicted 
by  the  model  are  shown  in  Tables  6  and  7,  for  the 
isotopes  of  Titanium  and  Vanadium,  respectively.  It 
is  seen  that  good  agreement  between  theory  and  experi¬ 
ment  is  obtained  in  most  cases. 

In  conclusion,  although  only  the  deformation 
parameter,  6,  was  allowed  to  vary,  the  values  used  in 
the  theoretical  calculations  are  in  good  agreement  with 
the  values  calculated  from  B(E2)  measurements  as  listed 
by  Wapstra  (Wa  65)  for  even-even  nuclei  in  this  region. 
As  shown  by  Davidson  (Da  68),  the  B(E2)  value  for  the 
ground  state  to  2+  first  excited  state  in  even-even 
nuclei,  as  predicted  by  the  Rotational  model,  is  given 
by 
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Table  6 


Comparison  of  experimental  reduced  transition  probabilities 

45  4  7 

with  the  Strong  Coupling  model  calculations  for  Ti ,  Ti 
49  . 

and  Ti. 


Isotope 

Multipole 

B  (A) 

a 

E.  (keV)  ,1?- 

l  '  l 

+Ef  (keV)  ,l£ 

Experiment^ 

Theory 

37,  3/2*" 

0, 

7/2" 

E2 

139±14C 

129 

40,  5/2" 

o. 

7/2" 

Ml 

0.052±0.003C 

0.124 

1354,  9/2" 

0, 

7/2" 

Ml 

+0.039d 

°-105-0.034 

0.0002 

45 

Ti 

E2 

^  +280d 

J  -145 

91 

1468,11/2" 

0, 

7/2" 

E2 

178±43e 

127 

1521,  3/2" 

37, 

3/2" 

Ml 

0 . 183  ±0 . 030e 

0.103 

160,  7/2" 

0, 

5/2 

Ml 

0.0458±0.0001f 

0.0281 

E2 

239±41g 

253 

1253,  9/2" 

160  , 

7/2" 

Ml 

+  0  022*1 
°-i83-o;oii 

0.116 

E2 

49+56h 

-33 

258 

1253,  9/2" 

0, 

5/2" 

E2 

126I121 

88 

^Ti 

1445,11/2" 

160, 

7/2" 

E2 

146-3f 

162 

1445,11/2" 

1253, 

9/2" 

Ml 

0.32±0.081 

0.33 

1550,  3/2" 

0, 

5/2" 

Ml 

(1 . 4±0 . 4  )  xlO"4  3 

0.225 

E2 

.+10.03 
16.5  ~ 

-2.6 

0.25 

1550,  3/2" 

160, 

7/2" 

E2 

39 . 8±11 . 41 

0.92 

1381,  3/2" 

o. 

7/2" 

E2 

>32 

24 

1585,  5/2" 

o, 

7/2" 

Ml 

>0.0013 

0.0506 

49Ti 

E2 

>7.4 

15 

1762,  7/2" 

o, 

7/2" 

Ml 

>0,00074 

0.0157 
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Table  6  (cont'd) 

a)  Units  are  y2  for  B(M1)  and  e2fm4  for  B(E2). 

b)  Assumed  mixing  ratio  6=0  except  as  noted 

c)  Lifetimes  from  (Ly  70)  and  for  37  keV  transition  ol_  =17 . 8 

(Je  68) .  ^ 

d)  6  =  +0.625^*2^  from  (Zu  72). 

e)  Lifetimes  from  (Zu  72). 

f)  6  =  0.095-0.003  (see  text). 

g)  Weighted  mean  of  B(E2)  values  (Te  56,  Ri  62). 

h)  6  =  0.15-0.06  and  lifetime  from  (We  72). 

i)  Lifetimes  from  (We  72) . 

j)  6  =  -4.4*2*g  an<^  lifetime  from  (We  72). 
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Table  7 


Comparison  of  experimental  reduced  transition  probabilities 

4  7  4  9 

with  the  Strong  Coupling  model  calculations  for  V,  V 

,  51„ 
and  V. 


Isotope 

Multipole 

B  (A) 

a 

Ei(keV)  ,I?-*E 

f  (keV)  ,l£ 

Experiment*"* 

Theory 

88,  5/2" 

0, 

3/2“ 

Ml 

0.088±0.005 

0.0366 

47v 

146,  7/2" 

88, 

5/2" 

Ml 

0 . 371±0 .040 

0.204 

1295,11/2“ 

146, 

7/2" 

E2 

204-685° 

242 

90,  5/2" 

0, 

7/2" 

Ml 

0 .230±0 .013 

0.0033 

153,  3/2" 

0, 

7/2" 

E2 

19  7±3 

243 

153,  3/2" 

90, 

5/2" 

Ml 

0.0035±0.0001 

0.181 

49v 

1021,11/2" 

0, 

7/2" 

E2 

172±59d 

189 

1155,  9/2“ 

0, 

7/2" 

Ml 

0 . 016±0 . 004e 

0.0454 

E2 

106±28e 

135 

1155,  9/2" 

90, 

5/2" 

E2 

126±17f 

98 

320,  5/2" 

0, 

7/2" 

E2 

154±7.6g 

82 

Ml 

.  nn(.+0.0008h 
°-006-0.0014 

0.0147 

928,  3/2" 

o, 

7/2“ 

E2 

76±5g 

41 

928,  3/2" 

320  , 

5/2" 

Ml 

>  0.00011 

0.5231 

> 

1 — 1 

LO 

E2 

10  7±9g 

44 

1609,11/2" 

0, 

7/2" 

E2 

83±8g 

59 

1813,  9/2“ 

o, 

7/2" 

Ml 

2  96+25*^0^ 

-  2.26 

0.0312 

E2 

28±5g 

48 

1813,  9/2" 

320  , 

5/2“ 

E2 

24±5g 

48 

■ 
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Table  7  (cont'd) 

a)  Units  are  y2  for  B(M1)  and  e2fm4  for  B(E2). 

b)  Assumed  mixing  ratio  6=0  except  as  noted. 

c)  Lifetime  from  (B1  73) . 

d)  Lifetime  from  (Ta  74) . 

e)  6  =  “°*78^q*22  and  lifetime  from  (Ta  74). 

f)  Lifetimes  from  (Ta  74). 

g)  B(E2)  from  (Ho  70). 

h)  6  =  -0.32+°,2°  from  (Ho  70). 

—  U  .  D  / 

i)  6  =  -8 . 8+4  * 1  from  (Ho  70). 

'  —OO 

j)  6  =  +3*75^o*51  fr°m  (Ho  70)  and  lifetime  from  (Go  72). 
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The  intrinsic  quadrupole  moment  Qq  is  given  by  the 
relation 

Q„  =  —  6(1  +  0.16  6) 

°  /57  ° 

and  the  deformation  parameter  3  is  related  to  6 
6  a  0.953  . 

Using  these  formula,  the  average  values  of  \&\  obtained 

from  the  two  adjacent  even-even  isotopes  are  0.35, 

4  7  4  q 

0.28  and  0.19  for  V,  V  and  V,  respectively. 

The  corresponding  values  used  in  the  theoretical 
calculations  were  0.342,  0.291  and  0.191,  respectively. 
These  are  in  excellent  agreement  with  the  experimen¬ 
tally  calculated  values. 

The  agreement  is  not  as  good  for  the  Titanium 

isotopes  where  the  mean  values  obtained  from  the  B(E2) 

4  5  4  7 

measurements  are  0.30,  0.28  and  0.21  for  Ti,  Ti  and 

49  . 

Ti,  respectively.  The  corresponding  values  used  m 

the  theoretical  calculations  were  0.188,  0.236  and 

0.08,  respectively. 

As  shown  in  fig.  73  the  even-even  isotopes  of 
Titanium  and  Chromium  (NDS  73)  exhibit  spectra  typically 
associated  with  quadrupole  vibrations  of  the  core. 

The  Vibrational  model  predicts  a  one  phonon  excited 
2+  state  with  energy  -fim ,  and  a  degenerate  two  phonon 


excited  state  with  energy  2hoo  and  spins  0  +  /  2+  and  4  +  . 
The  level  spectra  in  fig.  73  approximately  show  this 
type  of  spectra  and  thus  it  would  seem  that  the  Inter¬ 
mediate  Coupling  Model,  as  described  in  Appendix  A, 
would  be  useful  in  reproducing  the  properties  of  the  odd 
even  isotopes  of  Titanium  and  Vanadium.  However,  using 
the  expected  single  particle  energy  separations  for  the 
1  f  |  -  1  f  -j,  2p  j  -  1  f  ^  anc^  2  p  i  -  1  f  j  orbitals  a 
poor  fit  to  the  experimental  data  was  obtained.  Only 
for  unreasonably  large  coupling  parameters  could  the 
experimental  level  schemes  and  spectroscopic  factors  be 
reproduced  with  fair  agreement.  This  is  consistent  with 
the  use  of  the  Strong  Coupling  model  where  large  coupl¬ 
ing  parameters  lead  to  permanent  deformation  of  the  core 
The  Nickel  isotopes  have  received  much  theore¬ 
tical  attention  in  recent  years.  Both  Auerbach  (Au  67A) 
and  Glaudemans  et  al  (G1  72)  described  those  isotopes 

in  terms  of  strongly  admixed  spherical  shell  model  con- 

3  5 

figurations  where  neutrons  occupy  the  2p  j,  1  f  anc* 

1  56 

2 p  2  orbits  outside  the  Ni  core.  Glaudemans  et  al 

also  calculated  the  electromagnetic  transition  rates 

using  an  effective  neutron  charge  of  e^  =  (1.70±0.08)e 

and  neglecting  any  “*^Ni  core  excitations.  In  both 

59 

calculations,  the  level  scheme  of  Ni  is  well  repro¬ 
duced  for  the  first  three  excited  states,  however,  for 
higher  energy  levels  a  poorer  fit  to  experimental  data 


was  obtained. 
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Since  the  even-even  adjacent  Nickel  isotopes 
(NDS  73A)  exhibit  spectra  typical  of  the  vibrational 
model,  it  was  of  interest  to  see  whether  the  Intermed¬ 
iate  Coupling  Model  could  be  applied  to  describe  the 

59 

properties  of  Ni.  Although  there  are  several  variable 
parameters  within  the  model  it  was  hoped  that  these 
could  be  fixed  from  experimental  data  and  then  only  the 
coupling  parameter,  £,  allowed  to  vary. 

The  one  phonon  core  excitation  energy,  fioo  = 

1.393  MeV,  was  obtained  from  the  average  of  the  ener¬ 
gies  of  the  first  2+  states  of  the  adjacent  even-even 

Nickel  isotopes  (NDS  73A) .  The  odd  neutron  has  avail- 

3  5 

able  to  it  the  single  particle  orbitals  2p  j,  1  f  and 
2  p  i-.  The  single-quasiparticle  energies  of  these  shell 

model  levels,  e, ^  5  -  e0  3  =  0.48  MeV  and  1  -  3  = 

lf2  2pJ  2pJ  2p2 

2.07  MeV,  were  obtained  using  the  spectroscopic  infor¬ 
mation  from  a  (d,p)  experiment  (Ch  73).  The  energy 
levels  and  transition  rates  were  then  obtained  using 
the  computer  program  CORPAR  (Ca  70)  and  the  resulting 
level  scheme  is  shown  in  fig.  81.  Using  up  to  a  possi¬ 
ble  3  phonon  core  excitations,  the  best  fit  to  the 
experimental  data  was  obtained  for  a  coupling  parameter 
£  =  2.135.  Also  shown  in  fig.  81  are  the  experimental 
level  scheme  (Hu  73)  and  the  two  theoretical  level 
schemes  obtained  by  Auerbach  (Au  67A)  and  Glaudemans 
et  al  (G1  72)  for  comparison. 


ENERGY  (MeV) 
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Fig.  81.  Experimental  and  theoretical  level  schemes 
of  ^Ni. 
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The  first  three  states  are  well  reproduced  by 
the  Intermediate  Coupling  Model  and  although  the  fit 
is  poorer  for  higher  energy  levels,  corresponding  ex¬ 
perimental  levels  can  be  found  for  the  predicted 
levels.  The  number  of  levels  below  2.0  MeV  is  also 
in  agreement  with  experimental  results.  Although  no 
effective  charge  was  used  for  the  odd  neutron,  the 
electromagnetic  transition  rates  are  in  fairly  good 
agreement  with  experimental  results  as  shown  in  Table  8. 

No  shell  model  calculations  exist  for  ^~*Zn. 

6  5  6  3 

However,  3Zn  should  be  quite  similar  to  Ni,  for 

which  extensive  theoretical  studies  are  available 

(Au  67A,  G1  72).  Both  nuclei  have  seven  neutrons 

7 

outside  the  closed  1  f  ^  shell  and  although  the  proton 

6  3 

shell  is  closed  for  Ni,  to  a  first  approximation  the 
two  extra  protons  in  ^Zn  may  be  assumed  not  to  con¬ 
tribute  to  the  structure  of  the  low  lying  states  in 
6 5Zn  (We  70)  . 

Weidinger  et  al  (We  70)  have  done  Weak  Coupling 
Model  calculations  on  ^Zn  considering  ^Zn  as  the 
core  with  a  single  neutron  coupled  to  the  core. 

Excited  states  of  the  core  were  restricted  to  a  single 
phonon  excitation  and  the  core-particle  interaction 
was  the  Thankappan  and  True  (Th  65)  spin-spin,  quadrupole- 
quadrupole  interaction.  Since  this  model  is  similar  to 
the  Intermediate  Coupling  Model,  it  seemed  worthwhile 


. 


to  apply  the  Intermediate  Coupling  Model  to  Zn. 

The  single-quasiparticle  energies  of  the  shell 
model  states  e.  3  -  e.  c  5  =  0.130  MeV  and  1  -  ^5 

1 P  2"  lf2  ^2  lf2 

0.450  MeV,  were  obtained  from  the  (d,p)  data  of 

Ehrenstein  and  Schiffer  (Eh  67).  A  value  of  -fioo  =  0.8 

MeV  was  used  for  the  one  phonon  core  excitation  energy, 

which  is  slightly  smaller  than  the  average  value  of  the 

energies  of  the  first  2+  states  in  ^Zn  and  ^Zn  (NDS 

73A) .  A  best  fit  to  the  experimental  data  was  obtained 

for  a  coupling  parameter  £  =  1.8  and  the  predicted 

level  scheme  is  shown  in  fig.  82.  For  comparison,  the 

experimental  level  scheme  (Ez  74)  and  the  predicted 

6  3 

level  schemes  of  Ni  from  the  work  of  Auerbach  (Au  67A) 
and  Glaudemans  et  al  (G1  72),  are  also  shown  in  fig.  82. 

The  first  three  levels  are  correctly  reproduced, 
however,  the  3/2  level  at  207  keV  is  not  predicted. 

The  experimental  spectroscopic  factor  for  this  state 
is  small  (El  67)  and  hence  the  3/2  state  produced  by 

5 

the  coupling  of  the  ground  state  1  f  level  to  the 
first  excited  state  of  the  core  may  be  the  corresponding 
theoretical  level.  The  large  energy  separation  between 
the  first  set  of  states  and  the  second  which  is  charac¬ 
teristic  of  the  Intermediate  Coupling  Model  is  also  well 
reproduced. 

The  predicted  electromagnetic  transition  rates 
are  shown  in  Table  8.  It  is  seen  that  for  the  low 
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Table  8 


Comparison  of  experimental  reduced  transition  probabilities 

59 

with  the  Intermediate  Coupling  model  calculations  for  Ni 

,  65„ 
and  Zn. 


Isotope 

Multipole 

■  B(a) 

a 

E.  (keV)  ,  I?+E 
i  '  l 

f  (keV) ,1* 

Experiment13 

Theory 

340, 

5/2” 

o, 

3/2" 

Ml 

0 . 0154  ±0 . 0016° 

0.0019 

E2 

23±12c 

16 

465, 

1/2” 

o, 

3/2" 

Ml 

0 . 016±0 . 005 

1.44 

E2  ( 6=°°) 

78±22 

72 

878, 

3/2" 

o. 

3/2" 

Ml 

0.13510.023d 

0.091 

59.7. 

Ni 

E2 

20±17d 

40 

1189, 

5/2" 

o. 

3/2" 

Ml 

0 . 065±0 . 041e 

0.038 

E2 

+470e 

122 

x  -122 

62 

1302, 

1/2" 

0, 

3/2" 

Ml 

0.101±0.010f 

0.124 

1339, 

7/2" 

o, 

3/2" 

E2 

43±12 

89 

1339, 

7/2" 

340  , 

5/2" 

Ml 

0 .0005  +  0. 0005g 

0.28 

E2 

+  240^ 
447-116 

0.76 

54, 

1/2” 

0, 

5/2" 

E2 

-,00+41  h 
l80-30 

27.9 

115, 

3/2" 

0, 

5/2" 

Ml 

0 . 0584±0 . 0012 

0.00011 

207 , 

3/2" 

0, 

5/2" 

Ml 

0.0071±0 .0003 

0.0023 

65^ 

Zn 

207, 

3/2" 

54, 

1/2" 

Ml 

0.053±0.0021 

0.0004 

E2 

1864±9601 

21 

769  , 

5/2" 

o, 

5/2" 

Ml 

>0 .00359 

0.02 

E2 

>1 . 699 

2.33 

1264, 

9/2" 

o, 

5/2" 

E2 

>10 

28 

207 


Table  8  (cont'd) 

a)  Units  are  y2  for  B(M1)  and  e2fm4  for  B(E2). 

b)  Assumed  mixing  ratio  6=  0  except  as  noted. 

c)  6  =  0 . 11±0 . 03  from  (Hu  73). 

d)  6  =  -0.09±0.04  and  lifetime  from  (Hu  73). 

e)  6  =  0.43±0.90  and  lifetime  from  (Hu  73). 

f)  Lifetime  from  (Hu  73) . 

g)  6  =  -8.1^2*^  and  lifetime  from  (Hu  73). 

h)  Lifetime  from  (Au  60)  and  aT~aK=6±1  from  (Ro  70). 

i )  6  =  0 . 24±0 . 06  from  (Ro  70). 


j)  6  =  0 . 20±0 . 06  from  (Ez  74). 


lying  energy  levels,  the  transition  rates  are  in  poor 
agreement  with  the  model  predictions.  However,  the 
mixing  ratios  for  these  transitions  are  not  known  and 
assuming  the  lowest  multipoles  for  the  decay  of  these 
states  may  be  incorrect.  The  use  of  an  effective 
charge  parameter  will  also  change  the  predicted  tran¬ 
sition  rates.  However,  for  the  higher  energy  y-ray 
transitions,  better  agreement  is  obtained  with  experi¬ 
mental  results.  Hence  the  Intermediate  Coupling  Model 

g 

cannot  be  completely  ruled  out  in  its  application  to 
since  some  of  the  properties  of  this  isotope  are  repro 
duced  by  the  model. 


CHAPTER  9 


CONCLUSIONS 


This  research  has  yielded  a  new  method  of  improv¬ 
ing  the  time  resolution  of  Ge(Li)  detectors  by  a  pulse 
shape  discrimination  technique.  The  technique  gives  a 
significant  improvement  in  the  time  resolution  without 
seriously  affecting  the  counting  efficiency  of  the 
system.  The  results  are  considerably  better  than  those 
of  previous  selection  techniques.  The  improvement  in 
the  time  resolution  leads  to  greater  accuracy  in  the 
determination  of  the  centroids  of  the  time  distributions, 
which  is  the  main  limitation  of  the  Centroid  Shift  tech¬ 
nique  for  measuring  lifetimes  in  the  tens  of  picosecond 
region . 

A  further  refinement  of  the  technique  was  the 
development  of  the  Triple  Constant  Fraction  Discrimina¬ 
tor  which  provides  a  convenient  method  of  obtaining  the 
improved  time  resolution.  It  also  allows  the  experimen¬ 
ter  to  reject  certain  events  which  are  mistimed  by  the 
electronics . 

One  of  the  greatest  advantages  of  the  present 
delayed  coincidence  technique  is  the  simultaneous 
collection  of  data  for  both  prompt  and  delayed  time 
distributions.  Thus  errors  are  avoided  that  may  occur 
due  to  electronic  drifts  and  shifts  between  alternate 
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experimental  runs  to  obtain  the  time  distributions 
separately.  The  use  of  two  targets  to  obtain  supple¬ 
mentary  'prompt'  reference  y-rays  in  order  to  define 
the  prompt  time  centroid  curve,  recognises  the  effect 
of  the  walk  of  the  centroids  of  the  prompt  time  dis¬ 
tribution  with  the  energy  of  the  detected  y-ray. 

The  use  of  the  two  target  rotating  system  to 
obtain  these  supplementary  'prompt'  reference  y-rays 
is  another  advantage  of  the  experimental  method.  The 
Target  Discriminator,  which  allows  the  experimenter 
to  separate  the  information  obtained  from  each  target, 
reduces  errors  in  background  subtractions  since  the 
peak  to  background  ratio  is  not  adversely  affected  by 
the  presence  of  the  second  target. 

The  method  of  analysis  where  the  centroid  shifts 
of  the  delayed  time  distributions  are  measured  relative 
to  the  prompt  time  centroid  curve  at  the  appropriate 
energy  eliminates  those  errors  that  would  occur  due  to 
the  walk  of  the  prompt  time  centroid  curve.  Analysis 
by  this  method  ensures  the  accuracy  of  the  present 
results  and  the  validity  of  the  technique. 

The  present  delayed  coincidence  technique  compares 
very  favourably  with  other  techniques  to  measure  life¬ 
times  in  the  picosecond  region  and  in  some  cases  distinct 
advantages  can  be  seen.  Using  gas  targets,  the  experi¬ 
mental  upper  limit  for  lifetime  measurements  by  the 
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Doppler  Shift  Attenuation  technique  can  be  extended 
into  the  region  of  10  ^  secs.  However,  the  use  of 
gas  targets  is  a  distinct  disadvantage  and  usually 
high  energy  heavy  ion  beams  are  required  to  obtain 
sufficient  recoil  velocities  to  observe  a  Doppler  shift 
in  energy  of  the  emitted  y-ray.  The  method  is  also 
restricted  to  higher  y-ray  energies  so  that  measurable 
Doppler  shifts  are  observed.  In  the  analysis  of  Doppler 
shift  experiments  the  function  describing  the  energy 
loss  of  the  residual  recoiling  nuclei  is  not  accurately 
known  and  thus  an  important  source  of  error  occurs  in 
the  results. 

The  Recoil  Distance  technique  is  again  generally 
restricted  to  reactions  induced  by  high  energy  heavy 
ion  beams  since  high  recoil  velocities  are  required  to 
measure  lifetimes  of  the  order  10  ^  secs.  Recently, 
improvements  in  the  technique  have  made  lifetime 
measurements  using  ^He  beams,  possible  (Br  73)  . 

Another  problem  with  this  technique  is  the  fabrication 
of  thin  targets  to  permit  the  residual  nuclei  to  recoil 
freely  into  vacuum.  The  mounting  of  these  targets  to 
ensure  an  optically  flat  surface  requires  very  precise 
engineering  and  measuring  techniques. 

The  above  problems  are  not  encountered  in  the 
present  delayed  coincidence  technique.  However,  there 
are  a  few  major  limitations  to  the  technique.  The 
first  is  the  time  resolution  of  Ge(Li)  detectors  which 
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becomes  increasingly  worse  with  decreasing  y-ray  energy. 
Thus  lifetime  measurements  in  the  region  of  50  ps  or 
less  become  more  difficult  for  low  energy  y-ray  tran¬ 
sitions.  Another  problem  is  the  availability  of  suitable 
'prompt'  reference  y-rays.  However,  once  these  have  been 
established  for  a  given  energy  region  the  problem  is 
overcome  to  large  extent.  The  third  major  limitation 
is  the  peak  to  background  ratio.  It  was  found  that  when 
this  ratio  decreased  below  1:2  large  uncertainties  were 
introduced  in  background  subtractions.  Hence  only  y- 
rays  with  reasonable  intensities  can  be  studied.  This 
problem  also  occurs  in  the  other  two  techniques  mentioned 
above  but  to  a  lesser  extent. 

From  the  theoretical  standpoint,  the  Strong  Coupl¬ 
ing  Symmetric  Rotator  Model  is  seen  to  predict  the  level 
schemes  of  the  odd-even  Titanium  and  Vanadium  isotopes 
very  well,  even  though  only  the  deformation  parameter 
was  allowed  to  vary.  This  provided  a  meaningful  test 

of  the  model.  The  low  lying  ground  state  triplets  in 
45  4  7  49 

Ti,  V  and  V  are  all  predicted  and  with  the  correct 

spin  assignments.  This  was  not  possible  with  the  pure 
7  n 

(1  f -^ )  configuration  shell  model.  The  observed  large 

energy  gaps  between  these  triplets  and  the  next  set  of 

negative  parity  states  are  also  well  reproduced  by  the 

model.  The  only  isotope  where  a  relatively  poor  fit 

4  7 

was  obtained  is  Ti  although  the  predicted  B(E2) 


, 
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transition  rate  for  the  160  keV-*g.s.  transition  is 
in  very  good  agreement  with  the  weighted  mean  of 
experimental  results.  Other  electromagnetic  transi¬ 
tion  rates  predicted  by  the  model  are  in  fairly  good 
agreement  with  experimental  results. 

Although  little  success  was  obtained  using  the 
Intermediate  Coupling  Model  to  predict  the  level  schemes 
of  the  Titanium  and  Vanadium  isotopes,  the  model  is 

reasonably  successful  in  its  description  of  the  pro- 
59  65 

perties  of  Ni  and  Zn.  However,  the  predicted 

electromagnetic  transition  rates  for  the  low  energy 

65 

transitions  in  Zn  are  in  poor  agreement  with  experi¬ 
mental  results. 

In  conclusion,  there  are  two  areas  where  further 
experimental  investigations  would  be  interesting  using 
the  present  technique  to  study  the  Titanium  and  Vanadium 
isotopes.  These  are  the  two  odd-odd  isotopes  of 
Vanadium,  ^V  and  ^V.  In  ~^V  the  determination  of 
the  lifetimes  of  the  320,  355  and  388  keV  levels  would 
be  possible  if  a  suitable  detector  could  be  obtained 
for  detecting  the  low  energy  y-rays  at  35  and  33  keV 
from  the  two  levels  at  355  and  388  keV,  respectively. 

In  the  present  work,  it  was  found  that  although  theore¬ 
tically  better  timing  of  low  energy  y-rays  should  be 
obtained  with  small  volume  planar  Ge(Li)  detectors  as 
compared  to  coaxial  Ge(Li)  detectors,  experimentally 
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this  was  not  the  case.  The  problem  here  may  be  that 
at  these  low  energies  the  noise  is  comparable  in 
magnitude  to  the  real  event  signals.  The  electrode 
area  to  detector  volume  is  greater  in  planar  detectors 
so  noise  contributions  are  larger  in  these  detectors  as 

compared  to  coaxial  detectors. 

48 

For  the  case  of  V  several  investigations  into 
the  lifetimes  of  the  levels  of  this  isotope  have  been 
reported  (Ba  63,  Au  67B,  Bo  71,  Hu  73A,  Ha  73).  The 
low  energy  level  scheme  shown  in  fig.  83  is  from  the 
work  of  Huber  et  al.(Hu  73A)  with  the  lifetime  results 
indicated  with  the  appropriate  references.  The  levels 
at  309,421  and  520  keV  are  obvious  candidates  for  mea¬ 
surements  by  the  present  technique.  Unfortunately,  the 

4  8  4  8 

Q  value  for  the  reaction  Ti(p,n)  V  is  -4.795  MeV  and 
it  was  found  that  using  d.c.  proton  beams,  beam  energies 
greater  than  6.0  MeV  were  required  to  excite  these 
states  to  acceptable  strengths.  Although  in  principle, 
pulsed  proton  beams  of  this  energy  are  obtainable  from 
the  University  of  Alberta's  van  da  Graaff,  many  attempts 
to  obtain  such  beams  were  unsuccessful  due  to  deterior¬ 
ating  ion  source  characteristics  and  accelerator  con¬ 
ditioning.  Thus  further  investigations  of  this  nucleus 
would  be  interesting  using  the  pulsed  beam  (p,n)  reac¬ 
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Fig.  83. 


Experimental  level  schemes  and  previous 
lifetime  results  for  the  states  in  4^V. 
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APPENDIX  A 


A  THEORETICAL  DESCRIPTION  OF  THE  ONE  AND  TWO  PARTICLE 

INTERMEDIATE  COUPLING  MODEL 


The  Intermediate  Coupling  in  the  Unified  Model 
was  first  outlined  by  Bohr  (Bo  52)  and  Bohr  and  Mottelson 
(Bo  53).  Later,  it  was  elaborated  by  Choudhury  and  T.F. 
O'Dwyer  (Ch  54,  Ch  67)  and  Ford  and  Levison  (Fo  55). 

The  model  is  of  considerable  interest  in  the  description 
of  the  low  energy  nuclear  properties  of  the  odd-mass 
nuclei  whose  doubly  even  neighbouring  nuclei  exhibit  a 
vibrational  spectrum.  As  described  extensively  in 
(Bo  53),  the  limiting  cases  of  this  model  correspond 
either  to  weak  coupling  (for  nuclei  near  double  closed 
shells)  or  strong  coupling  (case  of  deformed  nuclei) . 

1 .  One-Particle  Theoretical  Formulation 

In  this  model,  the  odd-mass  nucleus  is  treated 
as  a  coupled  system  comprising  of  a  doubly  even  core 
capable  of  performing  collective  quadrupole  oscilla¬ 
tions  of  frequency  and  an  extra-core  nucleon  which 
has  several  single  particle  levels  available  to  it. 

The  total  Hamiltonian  of  the  system  can  be  expressed 
in  the  form 
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H  +  H  +  H. 
c  p  int 


1 


where  H  is  the  Hamiltonian  for  the  core.  H  is  the 
c  '  p 

particle  Hamiltonian,  and  H^nt  represents  the  inter¬ 
action  energy  between  the  two  systems. 

The  quadrupole  oscillations  of  the  nuclear 
surface  of  the  core  can  be  described  by 


R(0<j))  =  R  [1  +  2  a  Y^(6  ,(p)  ]  2 

y  M 

where  a^  are  the  dynamical  variables  of  the  collective 
motion.  Provided  the  deviations  from  a  spherical  shape 
are  small,  the  collective  Hamiltonian  can  be  expressed 
in  the  form  (Ch  54) 


1 

2 


I 

y 


{b 


a 


y 


2  , 

+  C  a 
1  y 


3 


which  is  the  Hamiltonian  of  a  system  of  harmonic  oscil¬ 
lators  with  frequency 


0) 


4 


The  parameters  B  and  C  are,  respectively,  the  collective 
mass  and  the  nuclear  def ormability .  Quantising  the 
collective  Hamiltonian  by  defining  the  momentum  conju¬ 
gate  to  a  and  the  commutation  relations  for  a  and  n 
y  y  y  y 


3H 

n  =  — -  = 

fa  , n  .  ] 

=  ifiS 

y 

y 

yy 

M  da  M 

fa.  ,a.  .  1 

=  [n,n 

y 

y 

y 

y 

the  core  Hamiltonian  can  be  expressed  in  terms  of  creation 
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4* 

and  annihilation  operators  (b^  and  b^)  for  spin  2 
phonons  which  are  defined  by  the  equations 


b+  = 


y  (2BfiO) )  ^ 


1  -  [Boa  a,1"  -  in,,]  h  =  - - — r  [Boo  a  +  in+] 

I  i  11  I  i  -y  i  i  ii 


y 


y 


(2Bhoo) 


y 


y 


i  i 

and  solvinq  for  a  in  terms  of  b  and  b  yields 

y  y  y  y 


f 

a  = 

y 


2Bco 


[b+  +  (-) yb  ] 

y  -y 


The  core  Hamiltonian  then  becomes 


H  =  Z  (b+b  +  h  fiw 
c  y  y  2 

y  H 


and  has  eigenvalues 


EN  =  (N  +  J)  ho) 


where  N  can  be  interpreted  as  the  number  of  phonons 
associated  with  a  particular  state. 

For  a  given  single  particle  orbital,  the  parti¬ 
cle  Hamiltonian  H^  is  a  constant,  having  as  its  eigen¬ 
values  the  single  particle  energy  levels  which  are 
regarded  as  adjustable  parameters  in  the  theory. 

The  interaction  Hamiltonian  H.  ,  can  be  written 

mt 

to  first  order  in  a  (Mi  54) 


Hint  =  ’k<r)  E  “u  Y2(04>) 


8 


where  (r,0,<J>)  are  the  coordinates  of  the  particle.  In 
calculating  the  matrix  elements  of  8  one  needs  the 
values  of  the  radial  matrix  element 

k  =  <n 1  £  1  | k  (r)  I n£> 
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which  will  depend  in  general  on  the  quantum  numbers  n 
and  l  describing  the  initial  and  final  states  of  the 
particle.  However,  the  dependence  of  k  on  these  quantum 
numbers  is  not  critical  and  so  it  is  assumed  to  be 
constant.  It  is  convenient  to  write  the  Hamiltonian 
in  terms  of  phonon  creation  and  annihilation  operators, 

Hint  =-/f*u>£  E  (-)  u[b  +  (-J^  b+  ]  YtJ(e,,j,)  9 

y  y  M 

where  the  strength  of  the  interaction  is  characterised 
by  the  dimensionless  coupling  parameter 

?  -  (2^»"  k  •  10 

It  is  convenient  to  work  in  an  angular  momentum  repre¬ 
sentation.  The  basis  states  of  the  system  are  chosen 

as  the  eigenfunctions  of  H  +  H  ,  the  non-interacting 

c  p 

system 

|  jNR;  IM>  =  Z  (jRM-mm|lM)  |  j  M-m>  |Rm  >  11 

m 

where  j  is  the  angular  momentum  of  the  odd  nucleon, 
coupled  to  a  core  state  characterised  by  N  phonons 
coupled  to  a  core  spin  R,  to  give  a  state  of  total 
angular  momentum  I  and  projection  M. 

The  eigenfunctions  of  the  total  Hamiltonian 
(1)  will  be  linear  combinations  of  these  basis  states, 


of  the  form 
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| E  ; IM>  =  Z  A  (E  ; jNRI)  | jNR; IM>  .  12 

jNR 

The  expansion  coefficients  A (E^; jNRI)  and  the  corres¬ 
ponding  energy  eigenvalues  in  eq.  12  can  be  found  by 
calculating  the  matrix  elements  of  HT,  and  diagonalising 
the  resulting  matrix.  In  this  representation,  both  the 

core  Hamiltonian  H  and  the  particle  Hamiltonian  H  have 

c  p 

only  diagonal  matrix  elements 


<  jNR;  IM  |  H  +  H  |jNR;IM>  =  (N  +  +  £  . 

c  p  ^ 


13 


where  floo  is  the  phonon  energy  and  the  single  particle 
energy  for  the  particle  in  orbital  j . 

The  matrix  elements  of  the  interaction  Hamiltonian 


are : 


<  j  '  N  '  R '  ?  IM  |  Hint  |  jNR ;  IM>  =  -  J  j  fico£<  j  '  |  |  Y2  |  |  j  > 


/2  j  '+1  (-1) ^  '+R_I  W ( j j  'RR'  ;2I) 

{(-1)R'“R  /2R+1<NR| |b+ | |N'R' > 

+  /2 R  1  + 1  <N  *  R  '  |  |  b1"  |  |  NR>  }  14 

where  the  techniques  of  Racah  algebra,  using  the 
Clebsch-Gordon  and  Racah  coefficients,  and  the  Wigner- 
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Eckart  theorem  have  been  utilised  in  the  derivation 
(Ro  57)  . 


The  reduced  matrix  elements  for  the  spherical 
harmonic  and  the  boson  creation  operator  b  are 
the  only  terns  left  to  be  explicitly  defined.  The 
former  is  given  by 


/2ji+T<j'  |  |Y,|  |  j>  =  (-1)3'  hl  (2j+1^2j'+1>]  2  ( j  j  'h-h  |  2 


0) 


J5[l+(-)£+£']  . 


15 


* 

Since  the  interaction  Hamiltonian  must  be  a 
scalar,  using  the  definition  of  the  scalar  product  of 
two  irreducible  tensors  (Ro  57) 


H.  .  =  -k(Y^ 

mt  2 


V  - 


-k  z  <-)%*  Yt; 
w  2 


=  -k  Z  a  Y^  16 

y  2 

y 

where  since  R(0d)),  eq .  2,  is  real  =  (-)^  a  .  Eq.16 

y  ~y 

is  the  form  of  the  interaction  Hamiltonian  so  both  Y^ 
and  a*  are  chosen  as  irreducible  tensors  in  this 

y 

*  +  u 

calculation.  From  the  definition  of  a  ,  b  and  (-)  b 

y  y  ~y 

are  also  irreducible  tensors  but  not  b  alone.  This 

-  -y 

fact  is  very  important  when  using  the  Wigner-Eckart 
theorem  which  applies  only  to  irreducible  tensor  opera¬ 
tors.  The  definition  of  the  Wigner-Eckart  theorem  used 
is  that  of  Rose  (Ro  57) 

c-i'M'ItS  i  JM>  =  (  JKMk  |j'M,)<J'  |  |  T  |  I  J> 

-'IK1  1  '  K  1 
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Using  coefficients  of  fractional  parentage, 
detailed  calculations  of  the  method  used  in  obtaining 
the  reduced  matrix  elements  of  can  be  found  in 
(Ch  54)  and  (Fo  55).  The  result  is: 

<NR|  |bf |  |N'R' >  =  /N  [2N_1 (R' ) 2R|  }2N  R]  17 

where  [2N  ^ (R ' ) 2R j } 2N  R  ]  is  the  coefficient  of  frac¬ 
tional  parentage,  R'  and  R  are  the  total  angular 
momenta  of  the  (N-l)  and  N  phonon  systems,  respec¬ 
tively,  and  2  is  the  angular  momentum  of  each  phonon. 
A  list  of  the  reduced  matrix  elements  for  N  <  3  is 
given  in  Table  9  . 


Electromagnetic  Transitions 

The  transition  probability  for  radiation  of  a 
photon  of  multipole  order  L  and  frequency  cj  is 


T(qL) 


8tt  (L+l )  1 

- 2  fi 

L [  (2L+1)  !  ! ] z 


2L+1 

(|)  B(aL) 


18 


where  a  is  the  mode  of  the  transition,  either  electric 
or  magnetic,  and  B(oL)  is  the  reduced  transition  pro¬ 
bability  . 


B  (cjL  ) 


1 

2J.+1 

l 


E 

M .  Mj. 

l  f 


<f  I  M 


LM 


l> 


2 
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a 

The  values  of  the  multipole  operators,  used 

are  those  introduced  by  Bohr  and  Mottelson  (Bo  53)  for 
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Table  9 

Reduced  matrix  elements  for  the  boson  creation  operators. 


N 

R 

N' 

R’ 

<NR| |b+ | |N'R' > 

1 

2 

0 

0 

1 

2 

0 

1 

2 

v/2 

2 

2 

1 

2 

/2 

2 

4 

1 

2 

J2 

3 

2 

2 

0 

/7/5 

3 

0 

2 

2 

/3 

3 

2 

2 

2 

/  4/7 

3 

3 

2 

2 

/I577 

3 

4 

2 

2 

/II77 

3 

2 

2 

4 

Z36/35 

3 

3 

2 

4 

-/6/7 

3 

4 

2 

4 

/10/7 

3 

6 

2 

4 

/3 

the  system  resulting  from  the  coupling  of  a  single 
particle  to  the  quadrupole  oscillations  of  the  core 
surface.  In  this  formulation  only  the  electric 
quadrupole  and  magnetic  dipole  transitions  will  be 
considered.  These  multipole  operators  are: 


i)  Electric  quadrupole  operator: 


M 


Ze 


9,,  =  ( e«  -  -?)  r  Y  (GcJ))  +  -r—  Ze  R  a 
2  y  p  ..  2  Z  4tt  o 


A 


t 

o  -y 
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where  e^  is  equal  to  e  for  protons  and  0  for  neutrons 
r  is  the  radial  coordinate  and  R  the  nuclear  radius. 


o 

2 

The  term  -Ze/A  has  been  introduced  (Bo  53)  to  account 
for  the  orbital  motion  of  the  core  around  the  centre 
of  mass  of  the  whole  nucleus  which  will  give  a  matrix 
element  for  electric  transitions.  However,  this  will 
be  replaced  by  eefft  which  is  an  effective  charge  para 


meter. 


ii)  Magnetic  dipole  operator: 


M  /  3 

M,  =  y  h —  (g  s  +  g„  i  +  qn  R  ) 

ly  Ko  V 4 tt  ^s  y  ^  i  y  ^R  y 


21 


where  y  =  eh/2Mc  is  the  nuclear  magneton,  g0  and  g 
are  the  orbital  and  spin  gyromagnetic  ratios,  respec¬ 
tively,  and  gR  is  the  effective  gyromagnetic  ratio  of 

the  core  taken  as  equal  to  Z/A.  R  is  the  Z-component 

y 

of  angular  momentum  of  the  core. 
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The  initial  and  final  state  wavef unctions ,  | i> 

and  | f >  are  given  by 


|i>  =  I E ; IM>  =  I  A (E; jNRI) | jNR;IM> 

jNR 

and 

|f>  =  |  E  '  ;  I '  M  '  >  =  £  A(E;  j 'N'R'I' )  |  j 'N'R'I'M^. 

j  1 N  '  R  ' 

2 

In  the  formulation  the  matrix  elements  of  r  are  assumed 
independent  of  the  quantum  numbers  and  are  given  by  the 
constant  density  approximation 


> 


where  the  equilibrium  nuclear  radius  Rq  is  given  by 
R  =  1.2A1//3  fm  . 


Defining  X-^  and  X2  as  follows 


X 

X 


1 

2 


(e  +  eeff} 


_3_  /hw 
4  7T  J  2c 


Ze 


then  for  the  electric  quadrupole  operator  of  the  reduced 
transition  probability,  the  calculations  yield 

B(E2)  =  (21 '  +1)  2  A(E,  jNRI)A(E'  j  'N'R'I'  ) 

I+I'  jNR 

j  '  N  '  R  ' 

{6NN,<5RR'  (“}I  +j_R  Xl<j  '  I  I  Y2  i  I  j>y/2  j'+1  W(j  •  jI'I;2R) 
+  x26  j  j  ,6^,/Zffr+T  w(R,Ri'i;2j)  (-)I_j+R' 

[<N'R'  |  |bf  |  |NR>  +  (-)R'~R  <NRI  lfat  I  |N,R,>]  }  2. 
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The  reduced  transition  probability  for  magnetic  dipole 
operator  is: 


B (Ml )  =  y2  (21 ’+1)  Z  A (E ; jNRI ) A (E '  ; j ' NRI ' ) 
I+I'  °  jj'NR 


R+i-j-1, 


{gR(")XX,J*  J  J'6j  j  I6££'/R(R+1)  (2R+D  W  (RRI '  I ;  1  j  ) 


+  +  A2j+1)  (2  j  '+!) 


/£U+1)  (2A+1)  W(j 'jI'I;lR)W(££j ' j ; 1^ ) 


+  gs  (-)  "  '2'"  /  (2  j+1)  (2  j  ’+1)  W  ( j  1  j  I '  I ;  1R) 


Z-k+I '-R 


W  (¥sj  j  ;  1A)  } 
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2 .  The  Two-Particle  Theoretical  Formulation 

In  this  section  we  consider  a  system  of  A  nu¬ 
cleons  composed  of  a  collective  core  of  (A-2)  nucleons, 
coupled  to  a  single  proton  and  neutron  outside  the 
vibrating  core.  The  total  Hamiltonian  is 


H  +  H  +  H  +  H.  (p )  +  H .  (n)  +  V 

c  p  n  int  r  mt  pn 


24 


The  notation  is  same  as  in  section  1  but  the  extra 

term  V  is  the  residual  interaction  between  the  two 
pn 

extra  particles.  Again  it  is  convenient  to  work  in 
the  angular  momentum  representation  and  choose  as  the 
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basis  states,  the  eigenfunctions  of  the  core-plus- 

particles  Hamiltonian,  H  +  H  +  H  ,  and  represent 

c  p  n  c 

these  eigenfunctions  by  | j  j  J;  NR;IM>  where 

p  n 

|j  jnJ;NR;IM>  =  Z  ( JRM-p  p |  IM)  |  j  j  JM-p>|NRp> 

P  p  P  n 

and 


jpjnJM_p>=  E  (jpjnm  J  M-p)  |  jpm>  |  jn  M-p-m>  . 


The  eigenfunctions  of  the  total  Hamiltonian  will  be  a 
linear  combination  of  the  basis  states  and  are  repre¬ 
sented  as 


E;IM>  =  Z  A  E;i  j  JNRI  j  j  J;NR;IM>  .  25 

1  .  .  JpJn  1 Jp Jn 

JpJn 

JN  R 

In  this  representation  both  the  core  and  particle 
Hamiltonians  have  only  diagonal  matrix  elements. 


j  j  J ; NR ; IM>  = 
JpJn  ' 


(N  +  lotfco  +  eP  +  sn 
2  3  1 

26 


where  eP  and  £*?  are  the  single  particle  energies  for  a 
particle  in  the  orbital  j  . 

The  matrix  elements  of  the  interaction  Hamil¬ 


tonian  H .  ,  (p)  are 

int 
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< 


Hint(p)>  =  fioo£{<NR|  |b+  |  |N'R'>  +  (-)R*  R 


!R_!±1  <N'R'  |  |b+|  |NR>}<j  '  |  |Y2|  |  j  > 
f2R  +  1  P  Jp 


[  (2  j  '+1)  (2J'+1)  (2J+1 )  (2R+1 )  ]  ^  W  ( JJ  '  RR  1  ;  21 ) 

hr 


I-R'-j  +j 

W(JJ'j  j  '  ;  2  j  )  6TTI6..Ml6  .  (-)  n  p, 

JpJp  Jn  II'  MM'  i„"i 


-ln^n 
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A  similar  result  holds  for  (n)  except  that  the 

subscripts  p  are  replaced  by  n  and  the  whole  result 
is  multiplied  by  the  phase  factor 


(-) 


j  + j  ”J+ j ' + j ' -J 
Jp  Jn  Jp  Jn 


This  extra  phase  factor  comes  from  the  relation 


j  +  j  -J 

(j  j  m  m  I JM)  =  (-)  P  n  (j  j  m  m  |jM)  . 
JpJn  p  n1  JnJp  n  p1 


The  residual  interaction  Hamiltonian  V  is  assumed  to 

pn 

be  a  sum  of  the  scalar  products  of  two  tensors,  one  of 
which  acts  on  the  proton  coordinates  and  the  other  act¬ 
ing  on  the  neutron  coordinates  (Sh  63) 


00 

V  =  I  (Tk(p)  •  Tk(n)  ) 
pn  k=0 


=  l  M-t  T_r<p)  TI?(n)  •  28 


Using  this  relationship  the  matrix  elements  of  the 
residual  interaction  are 
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<Vpn>  6NN'  6RR'  6JJ'  ^MM'  5II  ' 


W(VpjnV*kJ)<^l  lTk(P}  I  I  V<jnl  'Tk(n)  M  jn>(-)Dp+Dn"J 


29 


However  this  term  can  be  assumed  to  be  small  and  can  be 
neglected  in  most  calculations. 

Electromagnetic  Transitions 

As  with  the  single  particle  case  the  following 
constants  are  defined: 


<r2> 


fR2  - 

5  o  ' 


and 


X- 


X- 


(e  .  +  R2 

l  eff  5  o 


3  /haj  _  2 

—  hr—  Ze  R 

4  n  J  2  c  o 


l  =  p  ,n 


The  electric  quadrupole  operator  for  the  two  particle 
system  then  becomes 

Mfy=  X1  Y2  <6pV  +  4  Y2  <6nV  +  X2  M  %  30 

and  the  magnetic  dipole  operator 


Mly=  yo  h  (gssy  +  gssy  +  g££y  +  gRV 


31 


The  reduced  transition  probabilities  can  then  be  cal¬ 
culated  to  give  the  two  horrendous  formulae 
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B(E2)  =  (21  1  +1 ) 
I+I ' 


2  A(E; j  j  JNRI)A(E' ; j ' j 'j'N'R'I') 

JNRjpjn  P  P 

J'N'R' j ' j ' 

JpJn 


XlWW^j^pl  ^2l|jp>W(j;jpJ’J,-2jn)W(I'IJ'J;2R) 


P"P 


i  j  “j  +J-J ' +R-I ' 
t(2j'+l)  (2J+1)  (2J,+1)]'S(-)  P  n 

IT 


n 


+  *lW6RR-6j  j-<inllY2ll3n>W<3^nJ’J;23p)W<I,IJ,J'-2R) 

PJP  ^ 


2  j  -j'-j  +R-I' 
Jn  Jn  Jp 


[  (2jp+l(2J+l)  (2J  '  +1 )  ]  2  (-) 


+  x26i  V6i  i'  6JJ'{~]1  J+R'v/  (  2Rt+1  )  W(R '  RI '  I ;  2  J ) 
•JpJp  -’n-'n 


[<N'R'|  Ib1^  |  NR>  +  (-)R'  Ry||—  <NR|  |  b+|  |  N'R'>]} 


and 

B(M1)  =  A  p2(2I'+l) 
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I+I 


.1.1 


JNRjpjn 

J'^n 


A (E ; j  j  JNRI) A (E ' ; j  j  J'NRI' ) 


P  n 


P  n 


{g  /R(R+1)  (2R+1)  W(RRI'I;1J)  (-)R+I  J  16JJ,6-  ,  »6  i 

-Ip-Ip  Jn-In 

+  gP  we-s55jpjp;Hfi)W(j^jpJ'J;ljn)W(I'IJ'J;lR) 

-  2  j  -J-J  '  +R-I  '  -ig-£ 

[|(2j  +1)  (2jp+l)  (2J+1)  (2J'+1)]'5  (-)  n  P 

S.  „  ,  6.  .  .  +gn  VHhhi  '  j  ;1£  )W(j’j  J'J;lj  )  W  (I '  IJ '  J;  1R) 

£  £  -\  I  JnJn  n  JnJn  Jp 

p  p  JnJn 

„  i  j'-j  +j  -2  J 1 +R-I ' 

[|(2jn+l)  (2j^+l)  (2J+1)  (2J'+1)]^(-)  n  n  P  P 


6  *  '  6  .  .  •  +  gPW(£  £  j'  j  ;13g)W(j  '  j  J  '  J ;  1  W  (I '  IJ  '  J ;  1R) 

n  n  JpJp  F  P  F  R  F  F 


[£  (£  +1)  (2£  +1)  (2  j  +1)  (2  j  ’+1)  (2J+1)  (2J'+1)  P 

P  P  P  P  P 


(_j>+j-R+i’-V3p+jp-jn-^} 
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APPENDIX  B 


THE  STRONG-COUPLING  SYMMETRIC  ROTATOR  MODEL 


Some  nuclei  show  an  energy  spectrum  characteris¬ 
tic  of  a  rotator  and  appear  to  have  a  non-spherical 
equilibrium  shape.  Such  nuclei  have  been  described 
successfully  by  the  Unified  Model,  which  exhibits  both 
collective  and  individual  particle  aspects  of  nuclear 
structure.  The  basis  of  the  model  is  the  strong¬ 
coupling  model  of  Bohr  and  Mottelson  (Bo  53)  where,  in 
odd-A  nuclei,  the  unpaired  nucleon  moves  in  a  field 
generated  by  the  deformed  core.  The  total  Hamiltonian 
consists  of  the  rotational  energy  of  the  core  and  the 
single  particle  Hamiltonian  H^  and  is  given  by 


H 


H  ,  +  H 

rot  p 


where 


H 


rot 


3 

l  2A 

k=l  Z9K 


R 


K 


1 


R^r  R2  and  R3  and  ^ and  ^3  are  resPectivelY r  the 
components  of  the  rotational  angular  momentum  R  ,  and 
the  components  of  the  moments  of  inertia  in  the  body 
fixed  coordinate  system.  In  general  R  is  not  equal  to 
the  total  angular  momentum  I  of  the  system  but  is  given 
by 
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I  =  R  +  J 


3 


where  J  is  the  intrinsic  angular  momentum  of  the  extra 
core  nucleon.  The  total  angular  momentum  will  be  a 
constant  of  motion,  however,  R  and  J  will  not  be 
because  of  the  deformation.  The  projection  of  I  upon 
the  laboratory  z  axis  is  M;  its  projection  on  the  body- 
fixed  3  axis  is  taken  as  K  whereas  the  projection  of  j 
on  this  axis  is  denoted  by  ft. 

If  the  system  possesses  an  axis  of  symmetry 
which  is  taken  as  the  3  axis,  then  K  and  ft  will  also 
be  good  quantum  numbers.  For  an  axial  symmetric  de¬ 
formation  we  take 


^1  ^2  ^ 


and  then  the  total  Hamiltonian  may  be  expressed  in  terms 
of  the  total  angular  momentum 


H  =  (J2  -  21 -j  +  J 2 )  +  H 


or 


„  tf2  /  r2  ,  .2 
H  =  2J  (J  +  J 


2I3 j3) 


h4 


2S 


(I+j_+  I_j  +  )  +  H 


P 

5 


where  the  quantities  I+  and  j+  are  the  operators  1^ ±  il? 
and  j1±ij2,  respectively.  In  the  last  two  equations 
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use  has  been  made  of  the  fact  that  since  a  quantal 
system  cannot  rotate  about  an  axis  of  symmetry  because 
it  is  impossible  to  define  an  azimuthal  angle  for  the 
system,  the  angular  momentum  of  the  rotator  about 
its  axis  of  symmetry  is  zero  and  hence  also  K  =  ft 
(Ro  70B) . 

The  second  term  in  eq.  5  is  the  Rotation- 
Particle-Coupling  term  or  the  Coriolis  interaction. 

In  its  absence  the  particle  would  be  unaware  of  the 
rotational  motion  of  the  axes  and  so  the  effect  of 
such  an  interaction  is  to  decouple  the  particle  from 
the  core.  In  general  this  term  possesses  only  off- 
diagonal  matrix  elements  except  for  the  case  K  =  ft  =  i- 
for  which  the  diagonal  matrix  elements  are  nonvanishing. 
This  follows  from  the  fact  that  since  the  nucleus  is 
assumed  to  possess  axial  symmetry  the  levels  are  dege¬ 
nerate  in  the  sign  of  K  and  ft. 

The  eigenfunctions  of  the  operator  (fi2/2$! )  J2 
are  the  normalised  D-functions  D(a$y)  which  are  the 
transformation  functions  for  the  spherical  harmonics 
under  finite  rotations.  These  satisfy  the  following 
relations 


1  °MK  i(i+1)  dmk 


Iz  DMK 


=  M  D 


MK 


I3  °MK 


=  K  D 


MK 


6 
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We  define  the  normalised  intrinsic  wave  function  as 


X 


fta 


(=|fta>)  where  a  is  an  additional  index  to  dis¬ 


tinguish  between  states  of  the  same  projection  quantum 

number  ft.  Because  the  potential  is  not  spherically 

symmetric,  the  average  force  on  the  particle  is  not 

2 

central  and  the  operator  j  in  eq.  5  is  not  a  cons¬ 
tant  of  motion  and  can  be  absorbed  in  the  particle 
Hamiltonian.  Thus  it  is  simpler  to  deal  with  states 
of  a  constant  angular  momentum  as  a  basis 


xfta  ^  cjfta  xjft 


There  is  an  alternative  representation  for  x^a  in  which 
the  constants  of  motion  are  the  orbital  angular  momen¬ 
tum  i,  its  component  A  on  the  3-axis,  and  £,  the  com¬ 
ponent  of  spin  on  the  3  axis 


xfta  j;A  a AAfta  XMft 


(xAAft“ lN£Afi>) 


8 


The  coefficients  c.^  are  connected  to  the  coefficients 

jfta 

a£Afta  ky  t^ie  angular  momentum  addition  relations 


c  .  n  =  £ 

j  fta 


( A^AE | j  ft) a 


AA 


AAfta 


In  the  absence  of  the  Coriolis  coupling  the  total 
wavefunction  $  is  a  properly  symmetrized  product  of  the 
D-function  and  single  particle  wavefunction  x^a 


il±l)  [d1  X 

. -  2’  1  MK  xfta 

16tt 


+  ^  °M-K  nx  X  — Orj 


•fta' 


<*>  =  ( 


10 


242 


£ 

where  IT  =  (-)  is  the  parity  of  x  (Pr  62)  .  0  is  not 

X 


an  adequate  wavefunction  when  the  Coriolis  interaction 
is  included  in  the  model  Hamiltonian.  This  couples 
states  which  differ  by  one  unit  in  K  and  hence  K  is  no 
longer  a  good  quantum  number.  One  may  use  $  as  a  suit¬ 
able  basis  and  expand  the  total  wave  function 


^(IrM;U)=  (2I±1)  z  A  (aK;  I)  [D 
16tt  Ka 


11 

where  u  is  an  additional  index  to  define  the  state. 

The  problem  remaining  is  that  of  the  extra 
nucleon  moving  in  the  potential  well  of  the  deformed, 
but  symmetric,  core.  It  should  be  expected  that  the 
form  of  the  potential  should  go  over  in  a  continuous 
fashion  to  the  shell-model  potential  when  the  deforma¬ 
tion  goes  to  zero  and  that  the  well-known  magic  numbers 
at  nucleon  number  2,  8,  20,  50,  82,  126  are  reproduced. 
This  last  condition  is  an  important  restriction  on  the 
form  of  the  potential.  In  our  case  the  Nilsson  (Ni  55) 
model  Hamiltonian  was  used  to  describe  the  single  par¬ 
ticle  motion 

9 

H  =  H  +  C  £  •  s  +  D  £ 


P 


o 


2 

Cl'S  is  the  usual  spin  orbit  tern  and  the  D  l  term 


gives  a  correction  to  the  oscillator  potential 


12 
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especially  at  large  distances.  C  and  D  are  parameters 
adjusted  to  reproduce  the  shell  model  states  when  the 
deformation  vanishes. 

Axial  symmetry  is  assured  by  taking  co-^  =  ^2  =  00 ' 
which  with  oo^  is  parameterized  by 

2  2  2 

co  =  coQ  (1  +  -j  6 ) 

13 

2  2  4 

C03  =  co^(1  -  j  6) 


where  6  gives  the  degree  of  deformation,  the  oscillator 
potential  becoming  spherically  symmetric  when  6=0. 
Requiring  the  volume  enclosed  by  any  given  spheroidal 
equipotential  surface  to  remain  constant  under  changes 
of  6  leads  to 


2p2  _  2p2 

W1R1  w2R2 


2p2 

U)3r3 


2 


V 

o 


and 


R1R2R3 


const 


where 


R.  are  the  semi  axes  of  the  surface. 
1 


Hence 


and 


2  2  2 

a)la32co3 


const 


u)o(6) 


16_ 

27 


o  -1/6 
6  J ) 


14 


where  w  is  the  value  of  co  (6)  for  6  =  0  and  is  related 
o  o 

O  *1/0 

approximately  to  the  nuclear  mass  by  hcoQ  =  41A 


MeV. 
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Introducing  new  coordinates 


i  =  1,2,3 


the  Hamiltonian  becomes 

H  =  H  +  +  C  l  -s  +  D  l2  15 

p  o  6 

where  Hq  is  a  spherically  symmetric  harmonic  oscillator 
Hamiltonian 

H  =  i(-A  +  p2)  16 

o  o  2  K 

and 

H5  -  -  I  71  *“o  P2  Y2(0*>  •  17 


The  representation  chosen  by  Nilsson  was  with 

2 

Hq  diagonal,  together  with  £  ,  £^  and  s^,  which  all 
commute  with  H  .  The  corresponding  quantum  numbers 
are  denoted  by  £ ,  A  and  £ .  The  base  vectors  used  are 
|^£A2>  where  N  represents  the  total  number  of  oscillator 
quanta.  One  has 

H  |N£AZ>  =  (N  +  |  N£A£>  .  18 

U  A  U 


Nilsson  further  parameterised  the  total  Hamiltonian  eq. 
15: 

H  =  H  +  R  19 

p  o  o 

where 

1  C 
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and 


R=r)U-2£*s-|jZ 


20 


with 


y  = 


2D 

C 


and 


6  “o<6) 

n  =  K  o 

0) 


U  is  obtained  from 


H~  =  -6fico 
6  o 


j  y|-  p2  y°(0c|>)  =  KfiwQnu 


21 


The  diagonalisation  of  the  dimensionless  matrix  R  is 
then  carried  out  in  the  chosen  representation.  R  is 
treated  as  a  function  of  n  and  is  diagonalised  for  a 
sequence  of  n_values  giving  eigenvalues  rN^(ri)  where 
ft  is  the  eigenvalue  of  j ^  =  A  +  £ .  The 

corresponding  energy  eigenvalues  of  the  total  Nilsson 
Hamiltonian  are  thus 


Em  .  ,  =  (N  +  4)  ftco  (6) 

Nilsson  2  o 


xfia) 


Nft 
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The  eigenvectors  are  then  an  expansion  in  terns  of  the 
basis  vectors 


xfta  ~  lfia>  £A  aMftcJNME>  *  23 

Associated  with  each  level  are  a  much  used  set  of  num¬ 
bers,  the  asymptotic  quantum  numbers  [Nn^Alft  and  the 
parity  where  n^  is  the  number  of  oscillator  quanta 
along  the  3-axis.  Although  n^  is  not  a  constant  of 
motion  it  becomes  very  nearly  so  for  large  deformations. 
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Incorporating  these  eigenvalues  with  those  of 
the  rotational  Hamiltonian,  the  diagonal  terms  of  the 
total  Hamiltonian  H,  eg.  5,  are: 


<IMRj|H|  IMKct>=  [1(1+1)  -  2k  +  £|C_; 


2 J 


jfia 


j  (j+D 


-  6K,3S(-> 


I+ig 


( 1+^) a]  +  E 


Nilsson 


24 


and  the  off  diagonal  terms  are 


<IMK'a.|H|lMKa>=-|3-  I  CjJJla,C  jfia{  [  6R,  >K.1+  (-^  «K,  (K_1}  ] 


[  (I+K)  (I-K+l)  ( j+K)  ( j-K+1)  ]  "2  + 


6K'  ,K+1[  (I~K)  (I+K+1)  (j“K)  (j+K+D^-fij^,  ^j(j+l)  } 

25 

where  use  has  been  made  of  the  following  matrix  relations 

<IMK| I+ | IMK±1>  =  [ (ITK) (I±K+1) ] ^ 

<jfi|j_|jfi±l>  =  [ ( j+fi)  ( j  ±£H1 ) ] ^ 

The  last  term  in  the  diagonal  matrix  elements  arises 
from  the  Coriolis  coupling  term  for  K  =  This  serves 

to  define  the  "decoupling  parameter"  a 

a  =  Z  (~)^+1"2  (j+k)  |C  .2  |  2  26 

J  1  jka 1 
J  J 

The  decoupling  parameter  can  have  a  very  strong  in¬ 
fluence  on  the  appearance  of  a  K  =  ^  rotational  band. 

If  it  is  small  the  band  will  display  the  usual  1(1+1) 


247 


spacing;  however,  for  larger  values  the  band  will  have 
an  anomalous  behaviour. 


Electromagnetic  Transitions 

Brockmeier  et  al  (Br  65)  and  Nilsson  (Ni  55)  have 
discussed  electromagnetic  transition  probabilities  with 
and  without  the  band  mixing  Coriolis  coupling  term, 
respectively.  The  formalism  of  Brockmeier  is  now  pre¬ 
sented  for  the  reduced  transition  probabilities  of  (a,L) 
radiation  between  the  band  mixed  states  (u  I)  and  (u'l1), 
where  u  is  an  additional  index  to  label  the  state 


B  (GL;pI->U  '  I '  )  =  C  (CJL) 


Z 

a  K 

a  'K' 


F  (a 'K' ;u ' I ' I | aK;uI) 

0  JL 


27 


where  C(aL)  is  a  constant  and 


F  (a'K'ju'I*  aK ; ul )  =  A  (aK;I)A  ,  (a ' K ' ; I ' ) G_  (a ' K ' ; aK) 

OJ_i  U  U  GjL 

T  1  +V  1 

[ (IL,K  K'-K  I'K' )+(-)  (IL,K  -K ' -K | I ' -K ' ) b  (a'K'jaK)] 

Q  ij 


28 


A  (aK;I)  are  the  coefficients  of  expansion  of  the  band 
mixed  states  in  terms  of  the  basis  states  eg.  11. 


i)  Magnetic  dipole  transitions: 


In  this  case  the  constant  C(aL)  is 


C  (Ml )  = 


efi 
2MC  74¥ 
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and  G^l  and  b^L  are  calculated  to  be 


GMl(a'K';aK)  SKK.[  6aa'gRK  +  £A  aMftaaMfta '  (gsE  +  g£A)] 


+  ^  I  vl 


K-K' 


^  a 


K  '  -K  |  ,  1  yj  £AA,  M'ft'a' 


^AA' (gs  gR)  +  ^ZE' (g£  9r) ^ (^~A<) U+A>) ]  29 


(_) 1+N 

bMl(a'^;a^)GMl(a-^;a^)  =  — —  ^  ~£A^a“M^a 


^  aCAUa0Aii^i 


[6Ao6A'o(gs  gR)  +  (6Ao6A'1+6A16A'o)  (gil  9R)  ^  ( £+1 )  ] 


As  before  K  =  ft  and  K'  =  5!'. 


ii)  Electric  quadrupole  transitions: 


In  the  case  the  constant  C(aL)  is 


C(E2)  =  S2(e  +  eeff) 


where 


30 


3  = 


IT 


rMcu 


The  calculated  relations  for  G^_  and  b^T  are 

<JL  CJL 


GE2(a'K';aK)  =  G^(a'K';aK)  +  G™~  (aK)  6aa  ,  6 


coll 


KK 
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coll  1 

where  G„0  (aK)  =  — ^  Q  and  Q  is  the  intrinsic  quadru- 
ia z  2  3Z  ^  o 

pole  moment  given  by 


Q  =  0.8  z  R  6(1+0.56) 
o  o 
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where  Rq  =  1.2  A  fm,  6  is  the  deformation  parameter 
and  z  is  the  atomic  number. 


G®|  (a'K* ;aK) 


Z 

£  A 

£'  A' 


<5ZZ,aA,A,fi,a,aAA^a 


(i«'/p2/nju 


(£2AA’-A  U’A' )  (£200  U'O) 
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and 


sp 

bE2(a'K';  K)GE2 (a'K' ;aK)  =  (-) 


K ' +J5+N ' 


'2' 1  V  ^  r  I  ^  n  I  A  ■  .  i  U 


£  Z-.-Z'  £'A'ft 


£  A 
£ '  A 


'1 


'a' 


(N£'/P2/:n£) 


/  2  A+l 

J2  A  1 +1 


( £2 A-A ' -A | £ ' -A ' )  (  £200  |  £ ' 0 ) 
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where  p  =  r/$  and  a  supplementary  condition  that  b^^ 

ij  Z 

may  be  non-zero  only  if  one  of  the  following  conditions 
is  fulfilled: 

1  3 

1)  K  =  j  ,  K '  =  j  or  vice  versa 

2)  K  =  K'  =  ~  and  a  /  a' 


£Afia 


